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The  light-scattering  problem  of  a  sphere  on  or  near  a  plane  surface  is  solved  by  using  an  extension  of  the  Mie  theory. 
The  approach  taken  is  to  solve  the  boundary  conditions  at  the  sphere  and  at  the  surface  simultaneously  and  to 
develop  the  scattering  amplitude  and  Mueller  scattering  matrices.  This  is  performed  by  projecting  the  fields  m  the 
half-space  region  not  including  the  sphere  multiplied  by  an  appropriate  Fresnel  reflection  coefficient  onto  the  half¬ 
space  region  including  the  sphere.  An  assumption  is  that  the  scattered  fields  from  the  sphere,  reflecting  off  the 
surface  and  interacting  with  the  sphere,  are  incident  upon  the  surface  at  near-normal  incidence.  The  exact  solution 
is  asymptotically  approached  when  either  the  sphere  is  a  large  distance  from  the  surface  or  the  refractive  index  of 
the  surface  approaches  infinity. 


INTRODUCTION 

Light  scattering  from  a  cylinder  wee  solved  independently 
by  Lord  Rayleigh1  in  1881  and  by  von  Ignatowaky2  in  1905. 
A  few  years  later.  Mie1  solved  the  scattering  from  a  sphere. 
Liang  and  Lo4  extended  this  theory  for  two  spheres. 
Kattawar  and  Dean1  confirmed  the  existence  of  resonances 
discovered  by  Wang  et  al who  measured  microwave  scat¬ 
tering  from  two  dielectric  spheres.  Bobbert  and  Vlieger7 
used  operators  to  solve  for  the  scattering  of  a  small  sphere 
near  a  surface.  Yousif  solved  for  the  scattering  from  two 
parallel  cylinders  and  recently  solved  for  the  scattering  from 
a  cylinder  on  a  surface.  Youaifs  results  have  been  verified 
experimentally  by  comparing  his  theory  to  the  light-scatter¬ 
ing  Mueller  matrix  elements  measured  for  a  cylinder  on  a 
surface.3  His  approach  to  the  problem  of  a  cylinder  on  a 
surface  has  been  adopted  to  solve  the  scattering  of  a  sphere 
on  or  near  a  surface.  The  approach  is  similar  to  that  taken 
by  Rao  and  Barakat,10  who  calculated  the  scatter  by  a  con¬ 
ducting  cylinder  partially  buried  in  a  conducting  medium. 

The  method  used  to  calculate  the  electromagnetic  radia¬ 
tion  scattered  by  the  system  is  to  address  how  the  radiation 
may  interact  with  the  sphere.  The  incident  plane  wave 
strikes  the  sphere  either  directly  or  after  interacting  with  the 
medium  at  the  aurfhes,  in  which  case  it  is  an  image  plane 
wave.  Far-field  rartintlnai  is  a  raauit  of  a  superposition  of 
fields  either  directfpfrem  the  sphere  or  from  the  image 
sphere.  The  fields  emanating  from  the  sphere  may  also 
reflect  off  the  lurftna  and  internet  with  the  sphere  again. 
The  magnitude  and  phaaa  a i  the  fields  about  the  image 
sphere  are  altered  by  the  reflections  that  taka  place  at  the 
surface  by  the  Fresnel  coefficients.  One  assumption  is  that 
the  scattered  fields  from  the  sphere,  reflecting  off  the  sur¬ 
face  and  interacting  with  the  sphere,  are  incident  on  the 
surface  at  near-normal  incidence.  The  exact  solution  is 
asymptotically  approached  when  either  the  sphere  is  a  large 
distance  from  the  surface  or  the  refractive  index  of  the  sur¬ 
face  approaches  infinity. 

Figure  1  shows  the  geometry  of  ihe  scattering  system.  A 
sphere  of  radius  a  is  located  on  the  r  axis  a  distanced  below  a 
plane  surface  bounding  two  media  of  different  refractive 
indicaa  oriented  perpendicular  to  the  r  axis.  The  incident 


radiation  is  a  plane  wave  traveling  in  the  x-t  plane,  oriented 
at  angle  a  with  respect  to  the  z  axis.  The  wavelength  and 
wave  vector  for  the  plane  wave  in  the  nonabsorbing,  non¬ 
magnetic  incident  medium  below  the  substrata  are  A  and  k, 
respectively.  The  complex  wave  vectors  for  a  plane  wave  of 
the  same  frequency  in  the  3phere  and  the  medium  above  the 
surface  are  h**  and  k — ■-  respectively. 

SOLUTION 

The  solution  to  the  light  scattering  from  a  sphere  on  or  near 
a  surface  is  treated  here  in  several  subsections.  In  Subjec¬ 
tion  1  the  scalar  wave  equation  is  solved,  and  the  vector  wave 
functions  and  Debye  potentials  Sic  developed.  In  Subsec¬ 
tion  2  the  scattering  coefficients  are  solved  for  a  general  case 
In  Subsection  3  the  specific  case  of  a  sphere-surface  system 
illuminated  with  plane-wave  incident  radiation  is  solved. 
Finally,  in  Subsection  4  the  Mueller  matrix  elements  are 
derived  from  the  scattering  coefficients. 

L  Wave  Equations 

The  starting  point  for  this  scattering  problem  is  the  vector 
wave  equation,  which  is  derived  in  almost  any  text  on  elec¬ 
tromagnetic  theory11: 

7>H  -  „  ^  -  0.  11.11 

at*  <* 

where  E  and  B  are  the  electric  and  magnetic  field  vectors, 
respectively.  With  the  assumption  of  a  time  dependence  of 
the  form  exp <-»«f )  these  equations  can  be  used  for  deriving 
the  scalar  wave  equation,  which  in  spherical  polar  cooHi- 

nates  may  be  expressed  as 


r*dr\  r*sim 


.  1  .j 

♦  — - ; - r  +  A  u  ■  '3 

r 2  sin*  A*2 
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Fig.  1.  Goomotry  of  featuring  lytUm.  A  »pbar*  of  radius  a  is 
located  a  distance  d  from  a  surface.  A  plane  wave  travels  in  the  x-t 
plane  at  ancle  a  with  respect  to  the  *  axis. 


A  separable  solution  of  the  following  form  is  assumed: 

u(r.d.*)-R(r)e(d)<N„).  (1.3) 

The  particular  functions  satisfying  Eqa.  (1.2)  and  (1.3)  are 


*(*>)  -  exp (imw). 


R(r)  m  tn(kr) 


Zn+in^r)i 


where  P?(cox  t»  are  aaeociated  Legendre  polynomials  and 
zn(kr)  are  either  spherical  Beeael  functions  of  the  first  kind, 
jn(kr),  or  spherical  Henkel  functions  of  the  first  kind, 
h  ,"(*r).  The  electric  and  magnetic  fields  may  be  expanded 
in  terms  of  the  vector  wave  functions  Mil  and  Nil  (i  *  1 
corresponds  to  ;„(Ar)’s  being  used  and »  •  3  corresponds  to 
h (*r)’s  being  used): 


E-y^M^+g.jissL, 

njn. 

H--Vp,X  +  ^  d-5) 

1&H 

njn 

Th#  vector  wave  function!  M!^  and  W,*  aaauxne  the  follow¬ 
ing  fonn: 

Mlfi, •  *[  r„(*r)P?(  cos  d)exp(tmw)] 

LSIflCF  J 

-  *  |r,(hr)  —  [P?(coe  d))exp(imd)J, 


Nlii,  ■  f  ^  tnikr)n(n  +  l)£j(coe  t»exp(im*)j 

+  *  Injkr)]  ~  [P?(cot  d)jexp(tm*)J 

+  *  {£  £  [n'ikr)]  slfs  ^(C0#  U-6) 

Rather  than  expanding  the  electric  and  magnetic  fields  in 
terms  of  the  vector  wave  functions  as  in  Eq.  (1.2),  it  is  often 
convenient  to  expand  the  fields  in  terms  of  the  Debye  poten¬ 
tials: 

E  -  V  7  X  (rp^uJi)  +  i  V  X  7  x  (i wO. 

a 

njn 

H  -  ~  Y  v  X  (r^O  +  Jvxvx  (rplu^). 
tuu  a 

(1.7) 

The  Debye  potentials  are  pwuSl  and  The  coefr.- 

cienta  pM  and  9m  are  common  to  both  potentials  and  vector 
wave  functions,  making  it  convenient  to  consider  the  scattar¬ 
ing  problem  from  either  perspective. 

2.  Scattering  Coefficients 

The  general  case  for  scattering  from  a  sphere  on  or  near  a 
surface  can  be  treated  by  repreeenting  the  incident  field  on 
the  sphere  se  a  superpoeition  of  terms  in  the  expansion  given 
by  Eq.  (1.5): 

E““«  j  £  + 

»• 0 

+  (11) 

*•0  »•-» 

The  coefficients  aOT  and  will  be  solved  for  the  special 
case  of  an  incident  plane  wave  in  Subsection  3.  The  radial 
functions  used  in  expending  the  vector  (unctions  for  the 
incident  radiation  are  the  spherical  Beasel  functions  ;„(*/■>. 
In  this  section  the  scattering  coefficients  will  be  solved  ir. 
terms  of  these  coefficients.  The  scattered  field  is  also  ex¬ 
panded: 

E“-y  y 

iH)  fit*—* 

H-«  — y  y  UM'X  +  r^fC.  (2.2) 

UUfi  4-f  ^ 

I)o0  ggawg 

The  radial  functions  are  the  spherical  Henkel  functions  of 
the  first  h  '."(hr).  The  fields  inside  the  sphere  are  also 
expanded: 

i»0  *«•- * 
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The  redial  functions  for  the  internal  fields  ere  the  spherical 
Bessel  functions  j*(k^*r)-  To  this  point,  the  analysis  is 
identical  to  those  given  by  Raft.  11-14  for  the  case  of  a  single 
sphere  (Mie  theory),  except  that  here  the  coefficients  a„„ 
and  bnm  will  include  a  reflection  from  the  surface.  In  addi¬ 
tion  to  the  three  fields  described  by  Eqe.  (2.1M2.3),  a  fourth 
field  is  incident  upon  the  sphere.  The  field  is  a  result  of 
the  scattered  field  of  the  sphere’s  reflecting  off  the  surface 
and  striking  the  sphere.  This  interacting  field  is  also  ex¬ 
panded: 

•  <1 

cmoc  ST  V  g  +  j.  m<ii 

n  *0  tt  »  -n 

H,nI «  —  V  V  h^yiz  +  (2-4) 

iWM  —  — 

',*0  <n*—  n 

The  radial  functions  used  for  the  interacting  field  are  the 
spherical  Bessel  functions  jn(kr). 

The  next  step  is  to  apply  the  boundary  conditions  at  the 
surface  of  the  sphere.  The  coefficients  for  the  internal 
fields,  c,,m  and  <&»,,  can  be  eliminated,  and  the  scattering 
coefficients  can  be  solved: 

^  Snm) 

x  ~  *>bpt 

m  (<*nm  +  tnmW> 

fnm  *  +  km) 

k^j*i«(ko)+n{k,^a)  ~  k^n(k^a)^ka) 

-  (b«n  +  ™ 

where  Hr)  and  f(r)  are  the  Riccati-Beaeel  functions  defined 
|  by 

¥«<r>  *  ri«(r)'  f«<r)  "  rk^(r),  (2.6) 

:  and  the  primee  denote  derivativea  with  reepect  to  the  argu¬ 
ment. 

The  interaction  coefficients  can  be  solved  by  using  the 
relatiooahip  between  the  scattered  field  and  the  interaction 
field.  The  interaction  field  ie  due  to  the  scattered  field’s 
reflecting  off  the  eurihee  and  striking  the  sphere.  The  in  ter - 
action  field  is  the  ftMpef  the  scattered  field.  Identically, 
the  interaction  potantial  la  the  image  of  the  scattered  poten¬ 
tial. 

Figure  2  shows  the  tiro  coordinate  systems,  one  centered 
about  the  real  sphere  and  one  at  an  image  location  of  the 
coordinate  system,  a  distance  of  2d  along  the  positive  r  axis. 
It  is  easiest  to  consider  potentials  when  solving  for  the  inter¬ 
action  coefficients.  The  interaction  potentials  are  the  scat¬ 
tering  potentials  centered  about  the  image  coordinate.  The 
interaction  coefficients,  however,  are  inverted  and  undergo  a 
reflection  at  the  surface.  The  inversion  ess  be  explained  by 
taking  advantage  of  the  nature  of  the  associated  Legendre 
polynomials  that  make  up  the  potentials;  ie., 

flrt«oe<c  -  d)|  -  (-i)'***,/Jy[coe(d)|.  (2.7) 


Fig.  2.  Image  coordinate  surftce  is  located  a  distance  2d  from  the 
sphere  coordinate  system  aloof  the  positive  t  axis.  Fields  in  the 
image  coordir.ate  system  are  inverted;  eg.,  the  image  of  the  incident 
plane  wave  travels  in  the  x~t  plane  at  angle  »  -  a  with  reepect  to  the 
t  axis. 


An  approximation  is  mads  in  order  to  account  for  the  inter¬ 
action.  The  interacting  radiation  is  assumed  to  strike  the 
surface  at  normal  incidence.  A  Fresnel  reflection  at  normal 
incidence  is  than  used  to  account  for  the  reflection  loss  at  the 
surface.  This  same  approximation  was  used  by  Yousif  to 
solve  for  the  scattering  from  a  cylinder  on  e  surface  that  has 
been  verified  by  experiment.  Justification  for  this  approxi¬ 
mation  can  be  seen  by  tracing  rays  from  the  image  coordi¬ 
nate  system  to  the  edges  of  the  sphere,  as  in  Fig.  3.  The 
mMimim  angle  at  which  e  ray  emanating  from  the  center  of 
the  imag-  coordinate  system  strikes  the  surface  is  30*.  and 
the  Fresnel  coeffideate  are  fairly  constant  from  normal  inci¬ 
dence  out  to  this  angle  for  most  optical  materials.  For 
highly  conducting  optical  surfaces,  such  as  mirrors,  the  Fres¬ 
nel  coefficients  are  nearly  constant  from  0*  to  90*.  and.  as 
the  refractive  index  approaches  •,  the  Fresnel  coefficients 
are  constants  and  there  is  no  approximation.  The  Fresnel 
reflection  and  transmission  coeffidantt  for  this  system  can 
be  written  as  follows: 


7W*,) 


Rrn(d,) 


TthW 


oo0  ^  “  ZJl  * 

(»,/n,)2  sin2  djl/2 

Z,  cos  d,  +  Z,[l  - 

(»,/«,)*  sin2  dj 1/2  ’ 

2Z,  coed, 

Zt  cae  tf,  +  Z,(  1  - 

(n,/n,)2  sin2  dju2‘ 

Z\  cos  -  Z,[l 

-  (rij/n,)2  sin2  d,]12 

Z,  cos  dj  +  Z7(  1 

—  (nl/»1)1  sin2  d,)12 

2Z,  cos  d, 

Zt  cos  d,  +  Zjll  ~ 

(n,/*,)2  sin2  d,)1 

where 
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Maximum  infli  of  incidence  on  the  surface  y  for  an  intar- 
{ ray  occurs  for  a  ray  traced  from  the  image  coordinate  system 
«d  at  the  center  of  the  image  sphere)  to  the  edge  of  the  sphere 
•  greatest  when  the  image  sphere  touches  the  real  sphere  (the 
e  touches  the  surface).  This  angle  can  be  no  greater  than  30*. 
Preanel  equations  are  fairly  constant  near  normal  incidences. 


TE  and  TM  subscripts  pertain  to  the  transverse  electric 
transverse  magnetic  cases,  and  the  subscripts  1  and  2 
ir  to  the  optical  parameters  of  the  media  below  and  above 
surface,  respectively.  The  interaction  terms  can  be  ex¬ 
ited  in  the  image  coordinate  system  as 

ere  the  primes  denote  quantities  in  the  image  coordinate 
tern.  Note  that  /!te(0*)  *  Rtm(0*)  ”  R( 0*). 
tt  this  point  the  scattering  coefficients  are  essentially 
ved.  but  they  are  coefficients  of  functions  in  a  different 
irdinate  system.  The  next  step  is  to  find  a  representation 
the  functions  in  the  unprimed  coordinate  system.  Each 
the  primed  functions  can  be  expanded  in  terms  of  func- 
ns  in  the  unprimed  system;  i.e., 

m 

-  ft(0*)(-l)*~#„, 

m 

-  RWK-nn.  £e<rKi  (2-iD 

ie  interaction  coefficients  can  be  solved  by  extending  sum- 
itions  over  n  and  m.  Collecting  like  terms  yields 

m 

f«,-«(o*)  y  (“D">"*^cj"'. 


n  •ifni 

ibstituting  these  expressions  into  Eqs.  (2.5)  yields 


L  -X)  J  \ 


(2.13) 


The  only  unknowns  are  the  translation  coefficients 
These  were  solved  by  Bobbert  and  Vlieger  by  way  of  recur¬ 
sion  relations.  In  this  paper  the  translation  is  along  the 
positive  r  axis  rather  than  the  negative  r  axis.  This  creates 
some  minor  differences.  The  coefficients  can  be  deter¬ 
mined  by  using  the  following  equations: 

cl,00'  -  (2n  +  l)l/Xl,(2M).  (2.14) 


while 

cJT10'  -  -(2»  +  l)‘X1'(2*d). 

[(«'  -  m  +  1  )(n'  +  m)(2n  +  l)Jl/lc!,"  -m'  ■ 

-  {(n  -  m  +  l)(n  +  m)(2»  + 

JT(n  -  m  -f  2)(n  -  m  +  l)'|l/* 

^TTaj  J  c-» 


-2*dj^ 


(2.15) 


(w  +  m)(n  ♦  m  -  l)!1* 

(2n  - 1)  I  C'-1  ’ 


(2.18) 


An'jn)  m 

cn  c*  » 


(2.17) 


(2.18) 


With  the  use  of  these  equations,  the  translation  coefficients 
can  be  determined,  and  the  only  remaining  unknowns  in 
Eqs.  (2.12)  are  the  coefficients  for  the  incident  field.  In 
Subeectioo  3  these  coefficients  will  be  solved  for  the  particu¬ 
lar  case  of  an  arbitrarily  incident  plane  wave. 


3.  The  Cam  at  Plane-Wave  Incident  Radiation 
The  incident  radiation  on  the  sphere  can  be  separated  into 
two  parts.  One  part  strikes  the  sphere  directly,  and  the 
other  part  reflects  off  ths  surface  before  striking  the  sphere. 
The  pert  that  reflects  off  the  surface  will  undergo  a  Fresnel 
reflection  on  striking  the  surface,  and  it  will  also  be  out  of 
phase  by  an  amount  exp(t2Adcoe  a).  The  Fresnel  reflection 
term  depends  mi  whether  the  plane  wave  is  TE  or  TM  with 
respect  to  the  plana  of  incidence  (tha  x-r  plane).  Both  cases 
must  be  considered  separately.  Equations  (2.1)  must  be 
separated  into  two  parts:  one  for  the  TE  case  and  one  for 
the  TM  case.  The  fields  can  be  expressed  by  including  both 
parts  of  ths  incident  radiation  as  follows: 

Eft  » (1  +  RTI(o)exp(2ikd  cos«)(-l)'**'"I 

x  i  i 

iH)  "l«-8 


(2.12) 
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£  ^  unm‘nvn  T 


**•<?  m*-n 


EfS,  *  U  +  RjMiaitipi^d  <-*  a)(— l)"*'"] 

*t  t  a™M- + oc. 

a«0  m«—  * 

H™  -  —  (1  +  /?TM<a)exp(2:/W  cos«)(-lV'*m] 


-  ^ 

xV  V 


6™MU'  +  c™N<n 

uim  kT1im  T  unm  a/n- 


(3.1) 


■he  factor  of  (-l)"4'"*  it  a  result  of  Eq.  (2.7),  since  the 
(fleeted  plane  wave  is  traveling  in  the  x-t  plane  at  an  angle 
f  r- a  with  respect  to  the  x  axis.  Coefficients  a™  a™,  6™, 
nd  61“  are  now  the  coefficients  for  a  plane  wave  traveling  in 
he  x-t  plane  at  an  angle  or  with  respect  to  the  t  The 
«th  to  the  solution  is  outlined  by  Stratton10  (the  case  for  a 
•  0*  is  solved).  The  result,  after  much  algebra,  is 

™  ~  +  t)  IK»  -  m)(n  +  m  +  l)ll/5/sT+,(coa  a) 

-  ((n  -  m  +  l)(n  +  mjJl/J/*7~1(coe«)| 

2in*-  dP^icot  a) 

’  n(n  +  1)  8^  ‘  (3‘2) 


lT«  _  in*7(2n  +  1) 
""  o(n  +  1) 


(coa  a) 


[~(w  -  m  +  l)(n  -oi  +  2) 
[  (2n  +  l)(2n  +  3) 

f(w  +  oi  +  l)(n  +  oi  +  2) 
L  (2.1  +  l)(2n  +  3) 


'«/*  _ _ _ 

;  +  Pjr,‘(eoB  «) 


2j*+i  cos  q) 

n(n  +1)  sin  a 


From  these  expressions,  Iqs. 

form: 


(3.3) 

(3.4) 

(3.5) 

(Z5)  take  on  the  following 


/S  *  j[l  +  *TT<«X-ir  "  exp(2 ikd  cos  a)]6™ 

+^(0*)  y 

n  J 

-  |[1  +  e*p(2x ikd  cos  a)]6™ 

+  )?TM(os)  y  (-lr^/T^c’''  "jo;.  (3.6> 

n  •i/nl  J 

In  practice  the  infinite  number  of  terms  in  these  equations 
are  not  needed.  Experimental  precision  determines  the 
number  of  terms  needed.  In  any  case  the  scattering  coeffi¬ 
cients  may  be  solved,  since  there  always  remains  an  equiva¬ 
lent  number  of  equations  and  unknowns.  The  next  step  in 
this  analysis  is  to  develop  the  Mueller  scattering  matrix. 


4.  The  Mueller  Scattering  Matrix 

Two  axae  are  used  for  determining  the  Mueller  scattering 

matrix,  ona  perpendicular  to  the  plana  of  incidence  (w)  end 

one  parallel  to  the  plane  of  incidence  (£).  The  metrix  will 

also  be  derived  for  the  far-field  solution,  where  hr  »  n*.  In 

this  limit  the  spherical  Henkel  functions  reduce  to  spherical 

waves: 

h'"(hr)  ~  1*.  (4.1) 

i hr 

The  field  brio*  the  surface  is  composed  of  two  parts.  One 
part  results  from  the  sphere  directly.  The  other  pert  results 
from  the  Held  reflecting  off  the  eurface.  This  part  is  cen¬ 
tered  about  the  image  coordinate  system  and  is.  in  fact,  the 
interaction  term.  In  the  far  field  the  contributions  of  spher¬ 
ical  wavea  can  be  added  without  translating  the  image  coor¬ 
dinate  system.  This  image  field  will  contain  a  Fresnel  re¬ 
flection  term,  s  phase  term  equivalent  to  exp<-2 ikd  cos  d>. 
and  «n  inversion  term  equivalent  to  (-1)"*".  The  field 
above  the  surface  results  only  from  the  scattering  directly 
from  the  sphere  that  is  reduced  by  the  Fresnel  transmission 
factor  whan  paesing  through  the  surface. 

The  starting  point  for  calculating  the  Mueller  scattering 
matrix  is  the  amplitude  scattering  matrix,  which  is  defined 
by 


Sj  s,1  f*  £tm1 

[e?\  -ihr  {_ 

St  SiJ[£?lj 

•2  -  |(1  +  »nt(«><-ir-  exp<2iAdcwa)j«£2 

+  ftn<0*)  jr 

r 

•2f  -  |(i  +  *m<«><-l)"*"  «P<2*kd  cos«»e™ 

+  Rtu(0*)  ^  (-l)',+"e™c!lB^)J<?:. 


The  scattering  amplitude  matrix  elements  ere  solved  by 
expanding  the  scattered  electric  fields  in  terms  of  the  vector 
wave  functions  and  then  expanding  the  vector  wave  func¬ 
tions  in  terms  of  the  polarization  directions  by  using  Eqs. 
(1.8).  The  scattering  amplitude  matrix  elements  for  Idl  >  r/ 
2  assume  the  following  form: 

•  * 

S,"Z  2 

n*0  m*-» 

x  [l  +  Rjti*  -  exp  (-2 ikd  coe  i>>) 
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Price.#) +  •!?-#:(«»  4  (4.3) 
|_sin  <>  dv  J 

rhs  elsmsnts  behind  the  surface  (Itfl  <  r/2)  take  on  the  form 


s,  -  Y  Y  (-dvTtiW) 

•  n 

s2»-«Y  V  (-i)v~rTM(i» 

n»0  m--* 

x  [St  ^  ^  +  ^  £  ^(c“  w} 

s,--iY  V  (-i>vm-rTM(d) 

i»0  m«-i 

S,-V  y  (-OVTraW 

n«0  *»«••« 

x  [S  Pr(c0" + *™  ^ />:(co* J)}  (4-4) 

For  the  transmission  case,  i)  is  the  ancle  incident  upon  the 
surface.  The  actual  field  inside  *Se  surface  is  determined 
from  Snell’s  law.  Most  measurements  are  taken  by  piscine 
a  detector  in  a:?  on  the  opposite  side  of  a  thin  slab  of  the  bulk 
material.  If  the  rear  surface  of  the  bulk  material  is  antire¬ 
flection  coated,  these  ere  the  elements  that  would  result 
The  Mueller  scatterin«  matrix  can  be  determined  directly 
from  these  elements  by  usinc  relations  givsn  by  Bohern  and 
Huffman 

A  few  symmetries  that  exist  simplify  the  calculations.  A 
careful  examination  of  Eqs.  (3.6)  reveals  that 


*f\rh  a  *'  e*m* 

/2  -  (-i 

-  (-i  r/Ef.  (4.5) 

where  rft  *  -m.  If  the  Mueller  scattering  matrix  is  mea¬ 
sured  in  the  plane  of  incidence  (*  *  0®),  then  from  these 
symmetries  and  Eqs.  (4.3)  and  (4.4),  S3  and  S,  reduce  to  zero, 
since  for  the  normalized  associated  Legendre  polynomials 

^"(x)  -  (-l)m/*(x).  i4.6) 

Of  the  Mueller  scattering  matrix  elements,  only  four  remain 
that  are  not  zero  or  simple  multiples  of  another  element. 
These  may  be  determined  with  the  foilowingyrelations: 

Su  -  l/j(IS,lJ  +  ISA 
S„  -  VjdSjl*  -  IS, I*). 

S33  ■  RefSjSj*), 

Su  -  ImiSyS,*).  (4.7) 

In  the  plane  of  incidence  these  four  matrix  elements  are 
sufficient  to  categorize  the  light  scattering  from  a  sphere  on 
or  near  a  plane  surface. 

CONCLUSION 

Relatione  for  calculating  the  light-scattering  Mueller  matrix 
for  a  sphere  of  arbitrary  siza  and  optical  parameters  on  or 
near  a  surface  of  arbitrary  complex  refractive  index,  illumi¬ 
nated  by  light  of  arbitrary  wavelength  and  incident  angle, 
are  provided.  When  the  matrix  is  determined  in  the  plane 
of  incidence,  the  Mueller  matrix  can  be  characterized  by  four 
elements. 

The  particular  functions  used  in  the  derivation  facilitate 
the  «-l«ttlatin»i  of  the  matrix  elements.  Incorporating  Eqs. 
(2.7),  (4.5),  and  (4.6)  into  a  scattering  program  will  signifi¬ 
cantly  reduce  the  number  of  calculations  and  the  run  time. 
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Light  scattering  from  a  sphere  on  or  near  a  surface:  errata 
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A  few  errors  in  the  far-fie!d  expressions  for  the  light  scat¬ 
tered  from  a  sphere  on  or  near  a  surface1  are  corrected  in 
these  errata.  The  equations  for  the  interference  coeffi¬ 
cients  are 


-  Rm<-irv*„u£.  <2.10) 

The  interference  coefficients  must  be  translated  by 
means  of  the  vector  wave  functions.2  3  This  interference 
field  may  be  expressed  as 
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where 


c: 


■  C.  ”  - 


2kd  n  +  m  +  1 
2/»  +  3  n  +  1 


C/,.1 


2W  n-m 
- - - - C„-i  , 

2n  -  1  n 


~2tkd 


■me. 


n(n  +  1) 

are  given  in  Refs.  2  and  3.  Equations  (2.11)  lead  to 


(2.11b) 
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The  expressions  for  the  incident  fields  are 
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The  scattering  coefficients  are 

eJi  -  j[l  ~  RTe(a)l-l)"""  exp(2i*d  cos  a)]a,T| 

»  ] 

+  RtE(0°)  2  <-1  -  elicc  "\q;. 

e™  ”  |[1  -  exp(2ikd  cos  a>!a~* 

•  1 

+  EtmIO  2  t-iv'-l/ISD.*  "•  -  -  ■  }<?;. 

fll  “  j(l  +  Rn(a)‘- 1)*""  exp(2i*d  cos  a'ib.’f 

+  RteIO*)  2  -  e;*D,’  -  q\ 
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S]  —  -t  2  2  *  —*)'*  exp(im^) 
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-  ^rwlO4)  2  (-1)**"[AT?C.-  -  -  el *DS  x  |[1  -  RTM(n  -  d)(-l)"*"’  exp<-2 ikd  cos  0>] 


The  scattering  amplitude  matrix  elements  below  the  sur¬ 
face  i  d[  >  ir/2)  are 

w  a 

Si  “  2  2  1  ~i>"  expiimtp) 


{(1  -  Ht 

rTC  — 

'  •""m  O'"! 


TE(r,  -  exp(-2iAd  cos  d)] 


x  -^-rPTlcos  d)  *  [l  -  /?TE(ir  -  d)(-l>" 
sin 

x  exp(-2i*d  cos  d)l«^-^P?(cos  d)\< 
dv  1 


Sj  ■  -i  2  2  <-*)"  exp(tmy) 


x  1(1  *  R-ttiiir  -  &){-!)'"*  exp(-2iAd  cos 


X  PTicos  tf)  +  (1  +  Rrniir  -  d)(— 1)*** 

sin  0 


x  exp(-2tAd  cos  tf)!/-! 


■^PTlcosd)}. 


eTEm 

x  — — —  P^lcos  +  (1  +•  PTv(ir  - 

sin  d 

x  exp(-2iAd  cos  di]/',T|-4p"(cos  d)  [  ■ 
&o  ! 


S«  »  2)  2  <_*)*exp (imifi) 


(1  +  /?TE(fl-  “  d)(-Vn'm  expt -2ikd  cos  &i] 


x P^lcosd)  +  (1  -  RTt(n  -  tJH-lf 
sin  d  » 


x  exp(-2tAd  cos  d)]e™  — P?(cos  #>  >  ■  '4.3i 
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Coherent  fluorescent  emission  and  scattering  from  a  uniform  cylinder 
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The  coherent  light  inelasticaily  scattered  from  an  ideal  fluorescent  cylinder  illuminated  by  an  ar¬ 
bitrary  electromagnetic  wave  is  analyzed  using  an  extension  of  Rayleigh's  theory.  Scattering 
coefficients  are  solved  specifically  for  a  plane  wave  illuminating  the  cylinder  normal  to  the  axis  of 
the  cylinder.  The  resulting  inelasticaily  scattered  radiation  is  not  entirely  unlike  that  of  elastically 
scattered  radiation  from  a  perfect  cylinder.  These  results  are  illustrated  by  numerical  calculations 
for  several  specific,  but  arbitrary,  cases. 


I.  INTRODUCTION 

Inelasticaily  scattered  light  is  subject  to  the  morpholo¬ 
gy  and  optical  properties  of  the  particle  that  interacts 
with  the  incident  electromagnetic  radiation.  Experimen¬ 
tally,  morphology-dependent  resonances  have  been  ob¬ 
served  in  the  elastically  scattered  light  emitted  from  glass 
fibers  coated  with  fluorescent  dyes  by  Owen  et  al.1 
Fluorescence  enhancement  by  a  uniform  cylinder  was 
studied  by  Sekerak2  and  Abromson.3  Uniform  spherical 
particles  were  examined  by  Biswas  et  al.4  and  Benner 
et  al. 5  Theoretically,  Raman  and  fluorescent  emission  by 
a  molecule  embedded  in  a  dielectric  cylinder  has  been  de¬ 
rived  by  Chew,  Cooke,  and  Kerker.*  Kerker  and 
Druger7  theoretically  studied  fluorescent  molecules  em¬ 
bedded  in  spheres.  Das  and  Metiu*  took  a  quantum- 
mechanical  approach  in  examining  fluorescence  enhance¬ 
ment.  Ching,  Lai,  and  Young  examined  microspheres  as 
optical  cavities.9  For  a  uniform  distribution  of  molecules 
within  a  cylinder,  the  coherent  process  can  be  treated  by 
adding  the  resulting  amplitudes,  molecule  by  molecule, 
over  the  cross  section  of  the  cylinder.  A  shortcut  to  this 
summation  may  be  achieved  by  examining  the  polariza¬ 
tion  induced  at  the  new  frequency  within  the  cylinder. 
This  is  the  method  used  in  this  paper. 

In  this  paper  we  derive  and  examine  the  effects  of 
coherent,  inelastic  scattering  due  to  the  geometry  im¬ 
posed  by  an  infinitely-long,  circular  cylinder  composed  of 
fluorescing  material.  The  approach  taken  in  this  paper  is 
an  extension  of  Rayleigh’s  solution  to  the  elastically  scat¬ 
tered  light  from  aa  inflnitely-long,  circular  cylinder.10 
The  derivation  closely  follows  Bohren  and  Huffman's 
derivation,  and  we  represent  the  final  scattered  and  radi¬ 
ated  fields  in  terms  of  Moeller  matrix  elements. 1 1 

Inelastic  scattering  occurring  in  particles  is  a  process 
which  may  be  solved  in  three  parts.  First,  the  boundary 
conditions  of  the  particle  are  used  to  determine  the  inter¬ 
nal  electromagnetic  fields  within  the  particle.  Second, 
these  internal  fields  are  used  to  determine  the  fields  pro¬ 
duced  by  the  fluorescing  medium.  And  third,  the  bound¬ 
ary  conditions  of  the  particle  are  used  to  determine  the 
inelastic  fields  emitted. 

The  fluorescence  process  outlined  here  is  very  similar 
to  coherent,  Raman  processes.  For  incoherent  Raman 
processes,  the  internal  fields  are  not  intense  enough  to 
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stimulate  the  individual  molecules  to  emit  as  they  are  for 
coherent  Raman  processes  where  more  intense  internal 
fields  are  generated,  often  by  employing  a  second  laser. 
Although  our  theory  is  limited  as  a  model  for  Raman 
scatter  because  it  incorporates  only  an  incident  field  of 
one  frequency,  an  extension  of  it  to  multiple-field  interac¬ 
tions  would  follow  the  same  basic  procedure.  However, 
we  expect  the  main  qualitative  features  of  the  scatter  gen¬ 
erated  with  this  simplified  model  to  remain  unchanged. 

Knowledge  of  the  interaction  of  the  incident  electric 
fields  and  the  fluorescing  medium  which  produces  the 
emitted  fields  is  necessary  to  accurately  determine  the 
inelasticaily  scattered  radiation.  We  assume  that  the  po¬ 
larization  within  the  fluorescing  medium  is  proportional 
to  the  internal  electric  fields  due  to  the  incident  elec¬ 
tromagnetic  field.  We  also  assume  all  media  are  isotro¬ 
pic,  homogeneous.  With  these  assumptions,  the  deriva¬ 
tion  is  straightforward. 

II.  THEORY 

The  starting  point  is  to  derive  the  internal  fields  within 
the  uniform,  fluorescing  cylinder  having  the  same  fre- 


FIG.  I.  The  scattering  geometry  showing  the  fluore-ung 
cylinder  of  radius  a,  centered  on  the  origin,  oriemed  p j  r  j  . ; c i  to 
the  z  axis. 
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quency  of  the  incident  radiation.  This  is  done  by  apply-  ing  fields  of  the  incident  radiation  may  be  expanded  as 
ing  the  boundary  conditions  of  the  cylinder  on  the  elec-  follows  (see  Bohren  and  Huffman,  for  instance): 
tromagnetic  field  and  follows  the  same  procedure  used  to  x 

derive  the  scattered  field.  E,nc-1  =  2  .M'Ji  +  6„  .N'J’  ,  Ta» 


A.  The  fields  for  the  incident  wavelength 


H.n(- 1 


la) 

* 

2 

'lb) 

Figure  l  shows  the  scattering  geometry  of  the  system.  wiMi  ,  =  1  1 

A  cylinder  of  radius  a,  centered  on  the  origin  is  parallel 

to  the  z  axis.  The  first  step  is  to  derive  the  internal  elec-  where  k ,  =2 ir/A.,  is  the  wave  number,  cu,  is  the  angular 
tnc  fields  of  the  cylinder.  For  wavelength  A.,,  the  com-  frequency,  and  fx,  is  the  permeability  of  the  incident 
plex  refractive  index  of  the  cylinder  is  mt.  The  illuminat-  medium  at  A.,.  The  vectors  M',‘p  andN,''p  are  given  by 


2u\q  ) 

M "=(k]-hh'n  in — - — ~^~e, — Z^"'(p  )e.  *“«*’**'  , 

'  Pj 

ZU\Q  ) 

N^=<*/-A»2)l/2/lc,  ihZ\"'{p,  fc,-hn -  ^  L  e„  -H  fc,2  -  2 1 1 /2Z“ \p,  ler  . 


The  superscript  identifies  the  Bessel  functions  used  to 
represent  the  fields.  For  example,  i  *  1  corresponds  with 
the  use  of  the  functions  i—2  corre¬ 

sponds  with  the  use  of  the  functions  Z'2l(p)  =  Y„(p),  and 
i  =3  corresponds  with  the  use  of  the  functions  Zh3’(p) 
=  Hl,u(p):=‘J„{p)+iY„(p).  The  functions  representing 
the  incident  radiation  are  of  the  form  u„ine 
=/,(p,  )«'"  V  11  The  internal  fields  are  also  generated: 

E,n,.i=  i  <3a) 

ft  »  —  so  1  1 

X  .  13b) 

wiMi  ft*-*  1 

where  k\—kymy  is  the  wave  number  inside  the  cylinder 


medium  for  wavelength  A.,  and  p\  is  the  permeability  of 
the  cylinder  medium  at  wavelength  X,.  The  functions 
representing  the  internal  fields  are  of  the  form 
u„  m=JK(p\)e‘n*e  The  scattered  fields  are 

E^,=  i  ,  I4a> 

*  *  -  * 

2  •  '4b» 

a\f*\  b  -  -  » 

The  functions  representing  the  scattered  fields  are  of  the 
form  u„  ia~Hl"(p,)e‘”'e‘ 

Setting p, ,  =m’, //u„  f,*x,  sin0o,  17, -x , ( m \  - cos:0o >. 
and  x,  =  fc,a,  four  boundary  conditions  on  E  and  H  must 
be  satisfied: 


=  db.i"«i’?i«f<;(i?i)/Mr.i+c,,lAr,/»y,(ij,)A,pr ,  . 

and 

a,.l£l-M£l>  +  <b.  ll  W  ■ 


B.  Plane-wave  illaaiaadoa 

When  the  coefficients  a,  1  and  b„ ,  are  known,  then  the 
other  field  components  of  Eq.  (5)  can  be  solved.  For  the 
case  of  plane-wave  illumination,  Ar,  *  cos0o, 
p,=*k,i' sin^o,  and  p\»k\rtm\ -cos20o)l/2  where  0O  is 
the  angle  of  incidence  measured  from  the  2  axis.  We  con¬ 
sider  first  the  expansion  of  a  plane  wave  of  the  form 
E**3t£o«  '*  v  traveling  along  the  negative  x  axis. 

The  expansion  coefficients  for  the  TE  case  are 


a„.,=0,  6.  ,«£0<-i)"/*,  sin0o  . 

where  E0  is  the  magnitude  of  the  electric  field.  Next  •<.< 
consider  the  expansion  of  a  plane  wave  of  the  t  rn 
E=»er£0#  traveling  along  the  negative  1  1 3  ' 

The  expansion  coefficients  for  the  TM  case  are 

6., ,  -  0,  a,. ,  *  -  i£0(  ~ ‘  >V*  1  *in0o  • 

For  the  case  of  a  normally  incident  plane  -  •.  ■  • 
(S0“ir/2),  many  of  the  terms  in  Eqs.  (5a)-'5d'  ar.- 
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'since  h,  =Q).  The  internal  and  scattered  field  coefficients 
are  simple  expressions.  For  the  TE  case  (incident  electric 
field  parallel  to  the  z  axis), 

mifMiWvO-tySiV'lrix) 

f  =—  b  i  ~~ - - - - — - — •  .  ( oo) 

lnA  nA  7|) 

c,.,=e,,=0.  (8c) 


medium  (in  the  case  of  Raman  scatter,  a  corresponds  to  a 
gam  term).  We  consider  the  coherent  field  produced  hav¬ 
ing  the  specific  characteristic  wavelength,  angular  fre¬ 
quency,  permeability,  and  wave  number  tn  the  incident 
medium,  of  k2,  ai2,  p2,  and  k2,  respectively.  The  refrac¬ 
tive  index,  wave  number,  and  permeability  of  the  cylinder 
at  A.>  are  m,,  k2,  and  p\,  respectively.  We  first  expand 
the  field  oscillating  at  <u:  into  its  components. 


e„.2=  2 


a„  ,MlJl.  +  ’N'3’. 

”.py  t.p. 


1  1  Id) 


For  the  TM  case  (incident  magnetic  field  parallel  to  the 


z  axis), 

C"A  miy,(Tj,)«;(f,)-m,y;(Tj, >«.(§,) 

(  74J 

y,(^l)y;(»7l)-m,^(^)j„(?l) 

(9b) 

© 

II 

II 

r 

(9c) 

C.  The  fluorescent  aediu 

The  absorption  of  light  in  the  fluorescent  medium  is 
proportional  to  the  magnitude  of  the  electromagnetic 
field.  We  agree  with  the  assumption  of  Chew,  Cooke, 
and  Kerker  that  the  transmitted  electric  field  within  the 
cylinder  induces  a  polarization  within  the  medium  oscil¬ 
lating  at  a  shifted  frequency,  <u2;  i-e.,  the  distribution  of 
excited  fluorescent  molecules  emitting  a  particular  wave¬ 
length  k2  is  given  as  follows: 

?(r,<p,z,t)=saEintA(r,<p,z)e  '“J‘ ,  (10) 

where  a  is  an  excitation  characteristic  of  the  fluorescent 


h*.2= 


-'k'. 


T  a.  ,N'3’.  +6,  -M’3 


1 1  lb) 


tn  order  for  the  fields  to  behave  as  expanding  cylindrical 
waves,  the  functions  of  the  fields  created  by  fluorescence 
are  of  the  form  u„  n~H'n,\pi)e‘'"fe  V.  The  source  of 
these  fields  is  the  electric  polarization  expressed  tn  Eq. 
( 10).  The  coefficients  <j„2  and  j  can  be  solved  for  by  ex¬ 
panding  the  electric  polarization  in  terms  of  the  vector 
harmonics  of  generating  function  u„  „.  Since  the  electnc 
polarization  is  composed  of  vector  harmonics  of  generat¬ 
ing  function  u„  we  can  start  by  expanding  the  internal 
electric  field  vector  harmonics  in  terms  of  the  electric 
field  vector  harmonics  oscillating  at  cu2: 


Mm.=  T  .  ,  (12a) 

"■fit  ~  1  Pj 

-ia» 


As  outlined  in  Stratton,12  we  first  develop  orthogonal  re¬ 
lationships  for  the  vector  harmonics  by  introducing  two 
new  sets  of  vector  harmonics: 


(M  ^)*=(k/-h2)1/2 


Pi 


IN'1’  )*  =( A ? — A 2 )l/2//c. 


- ihZ'"'(p, - hn -  - -  e, +{kf-h2 )xnZ^(p, It, 


i(  n  f  + hi) 


1 1 3a) 


13b> 


The  generating  functions  for  these  vectors  are  of  the  form  u„., =/,(pi)e'"^e'*3Z.  The  orthogonal  relationships  are 

/o2  "MV.;-  • '  ^  a  -  a  ?  -  A 1 ) 1  ,  (p,  _  ,(pi )  -»-yn  ^ ,( p',  ,  ( Pi )  1 .  ua> 

r2’,N'1',-(N1V  .rdfl,  -h])uHk?  -h\)wl 

j0  *.p,  -»,pj  ^  *•*  k\k\  1  * 

~[j„ .  ,(p;  )/„ .  mi +/„  «.  m  )J„ +1<pi>  j 


+(k'12-h|)1/2(ki2-hj)1/2/,(p'1)y,(pi) 


/” 2>r M'2'. *( N*  ’  .  )*dg>=  /2’,Nai.-(Mu.1  .  )*d? 

J  0  *-*l  *  f>l  J0  *  n  f*l 


ith 

= s*.*'  *  i2  -  * !  _  ,<P'2  >//,  _  ,(p;  > -y,  *  ,<p'z)«,  ^  ,(pi >1  • 

K  2 


Ubi 


14C I 
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By  scalar  multiplying  (12a)  and  (12b)  by  ( M^.  )*  and  integrating  both  equations  over  the  fluorescing  cylinder,  then  re¬ 
peating  the  process  with  ( N'^1  \  )*  we  can  solve  for  the  following  coefficients: 

raM,1'.-(M"’,)Vdr  fV!l,-(N"’,)V<fr-  fV'.-IN11'.  )Vdr  faN,21.<M'".)Vdr 

J0  it.p,  » -Pj  J  o  «.p:  "P:  Jq  <1.0,  »P; _ Pp  "P;  l'.P; 

~6‘': -*i  rJM'J,.-iM“'.)*rdr  fV:’.-(N,".lVdr-  f  aN,2.-iM,".)Vdr faM‘2l,-(N',,.)V dr  ’ 

J0  ».p,  /t.pj  Pq  "-Pi  ".p:  J0  ’’•P:  ’-P2  ■'0  "Pj  "Pj 


fV2’.- 

'V". 

•(  N* 1 1 .  )*r 

dr- 

- /V". 

)*rdr  I 

rv2-. 

•(Nni. 

iVdr 

p  0  "P; 

d.p,  j 

0 

J 

fl.ps 

p  0  ,,pi 

H.pv 

j 

0  "-pj 

n.p. 

rv21. 

•(M'".)Vdr  I 

fV2’. 

(N'".)V 

dr- 

-  TV2’, 

V  dr  / 

’aM'2'.- 

(N'1',) 

*rdr 

J  Q  "P; 

n.p2  J 

0  "-Pi 

n,px 

•A  0  '••Pi 

n.p. 

J 

0  "-Pi 

'l.p. 

{'■N"UM"\)-rdr  faNiV.-(Nll,.),rdr -  m,)Vdr /V’.-fN  u\)Vdr 

•'0  "-Pi  "-Pi  «  0  "  Pi  "•pj  J  0  "‘Pi  "’P?  *  0  "-Pi  *’Pj 

C’’~’  *>•*'  f‘kti,.-(M,t>.)*rdrf'N'i,.-(Ntt,.)mrdr- /V2'  •(M‘l,.)Vdp/'M'I,.-(N,,’.)Vdr  ’ 

0  ff,P2  •  0  '••Pj  *0  '••Pj  *0  *•/*}  n,^2 


raM,2’,-(M"',)Vdr  f,Nu,.*(N‘".)Vdr-  raN'",-(M"),)*r dr  f‘,Mt2'.-(N"’,)V dr 

J0  >i.pj  *.p2  Pq  <i.p,  ».Pj _ P  Q  ’’•Pi  "P; _ Pp  11 -P;  "  P; _ 

raM'2'.-(M",,)*rdr  raN'i’.  (N",,)*rdr-/aN'2’,  (M'".)Vdr /aM'2l,-(N'1,,)*rdr  ' 

J  0  .V.pj  n.p*  Jq  n.pj  Jq  q  "*P2  "‘Pj 


The  result  that  /i,=/i2  has'  some  interesting  consequences.  For  plane-wave  illumination,  the  result 
Ji,  =/i2  =  fc,  cosdo=fc2cos02  means  that  the  scattered  light  of  wavelength  X2  will  not  travel  along  the  path  of  the  light 
of  wavelength  A.,  unless  the  light  is  normally  incident  on  the  cylinder.  As  shown  in  Fig.  2,  a  cylinder  illuminated  by  a 
plane  wave  scatters  light  elastically  into  a  cone  at  angle  0O  with  respect  to  the  z  axis.  The  inelastically  scattered  light 
travels  along  a  cone  at  angle  d2  with  respect  to  the  z  axis. 

When  A,  =0,  which  is  the  case  of  a  plane  wave  at  normal  incidence,  the  coefficients  are  greatly  reduced,  since  the 
orthogonality  condition  forces  ( 14c)  to  zero: 


_  *1  fna[J^l(k’lr)J,_l<k'2r)+J^l(klr)J^,(k’2r)}rdr 
fc'2/J’[»1,_1(k'2r)/._1(ltir)+/f.  +  l(fcir)/(l  +  1(*'2r)]rdr  ’ 

f0X‘}-(K'}/rdr  ^  k\  f°J„(k\r  )Jn(k\r)r  dr 
°'~6 f*N,2,.-(N<,’.)Vdr*  kj  f*JI,(*Jr)/.(jkJr)rdr  ’ 

p  o  "*Pi  ,-Pi  •  o 

B,=C,=  0. 

The  integrals  in  Eqa  (16m)  and  (16b)  are  Lommel's  integrals  and  have  the  following  solutions:11 
rxJJk’jX  \J,lk\x  /■  ■■  ^[kflUk'id  )Jn(k\a )—k\j;(k\a )]  , 

0  l  AC  j  J  — M  AC  j  / 

/%y,(*ix)/f,(*ix)dx  =  {[fl2j;(k2a)/f;(k^)- L2--^  hm'pWi*)— ^777  • 


The  coefficients  of  the  emitted  electromagnetic  fields  may 
be  written  as 


ai,lSM  ~  •  ciA'^n  ' 
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(18a)  where  e\  is  the  cylinder  permittivity  at  the  emitted  fre¬ 
quency.  An  induced  Held  within  the  cylinder  is  needed  to 
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satisfy  the  boundary  conditions.  This  field  corresponds 
to  the  field  created  by  fluorescence  reflecting  off  the  inner 
surface  of  the  cylinder.  Using  functions  of  the  form 
2  *  induced  field  may  be  ex¬ 
pressed  as 


e ,n«u=  2  . 


Hind.;  = 


Finally,  the  scattered  fields  may  be  expressed  as 


•  <20b) 

The  generating  functions  of  the  scattered  fields  are  of  the 

form  un^1=H'n"{p2)e,n*e'''i‘. 


FIG.  2.  The  scattered  radiation  from  a  cylinder  illuminated 
D.  Boundary  conditions  at  the  wavelength  of  flnorescence  by  a  plane  wave  traveling  in  the  x-z  plane  at  an  angle  t)0  with 

respect  to  the  z  axis.  Elastically  scattered  light  travels  in  a  cone 
Setting §2~*2  sinfl2,  Tf2-Xj(m{-cos%),  at  angle  9q  with  respect  to  the  z  axis.  Inelastically  scattered 

and  x^-kid,  four  boundary  conditions  on  E  and  H  must  light  travels  in  a  cone  at  angle  Bt  with  respect  to  the  z  axis  deter- 

be  satisfied:  mined  by  Ic,  cos0o=**2  cosOj. 


Si )  ]'  +  /„.2*  inH;  (§  ,)/*:. 

bnWiHn'{Vi^m1+d%1r}]Jn[r)1)/m1=fn  j\Hl,u(^1)  , 


For  the  case  of  a  normally  incident  plane  wave 
(d0— ir/2),  many  of  the  terms  in  Eqs.  (21a)-(21d)  are 
zero,  and  the  scattering  amplitude  coefficients  are  simple 
expressions.  For  the  TE  case. 


fn.l~bi.l 


*,.2  =  0  • 


For  the  TM  case, 


?,.2~“,.2 


/,.2=0.  (23b) 

it  can  be  shown  that  the  following  relationship  holds  for 
these  coefficients: 


•  24a» 

/_„.2=*(-l )"/„.*  .  '24b» 

With  analytical  expressions  for  these  coefficients,  the 
solution  to  the  scattered  field  from  a  coherently  fluoresc¬ 
ing  cylinder  is  complete.  Analytical  expressions  may  be 
derived  for  any  illuminating  radiation  expressed  by  the 
coefficients  a„  (  and  t  by  solving  the  simultaneous 
equations  (Sa)-(Sd)  and  <21a)-(2ld). 

For  large  arguments,  the  first-order  Hankel  function  is 
given  asymptotically  by  the  following: 

jl/2 

H'„"lp)~  —  «'*(  —/■)*€ \p\  »n:  .  25) 

irp 

At  large  distances  from  the  cylinder  \k:r  sine).  !  .the 
scattered  electric  field  may  be  expressed  as  follows 
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l(k\—h\) 


m 


1/2 


—  flr/4) 


2e"'*{~t)n 


h,  ^  (k2-h\ )]/2 

-  ie„  .e,  T-fn.iCp -*■ - - - /„  :c: 

*2  k2 


2b> 


E.  Scattered  light  for  normally  incident 
plane-wave  illumination 


We  now  look  at  the  specific  case  of  light  scattered  in 
the  far  field  from  a  fluorescing  cylinder  illuminated  by  a 
normally  incident  plane  wave.  We  examine  the  light  in 
two  separate  polarization  states:  the  TE  state  corresponds 
with  light  polarized  along  the  vector  e,  and  the  TM  state 
corresponds  with  light  polarized  along  the  vector  — e,. 
The  amplitude  scattering  matrix  is  given  in  the  far  field 
by 


2-te.»c* 

—  /? 

2 

*0 

irk2r 

Tx 

T< 

X 

Ty 

Ty 

i  i4c,p  — 

e  1 


ETE,inc 

ETM.inc 


ir/4) 


where 


Tt 


an 

2 


hi. 

Eo 


( —  i)V"*/*,2 


(27) 


/o.2  +  2  2  (~i/AiCos(n<p) 

n  m  I 


7*2-  2 


(28a) 


k  j  » 

~<T~  *02+2  2  ('')\jCO*(«p)  , 

*0  „-( 


(28b) 


r3  =  r,=0.  I’8c> 

The  four  independent,  nonzero  Mueller  matrix  elements 
may  be  expressed  by 


Tu  =  {UTx‘r+  Tv2)  , 

<29a) 

Tn  =  \(  ]T x\2 -  •  T2\2 ) /Tu  . 

>29b) 

Tjj-Reir.rp/r,,  , 

'29c) 

rn-imir.rji/r,,  . 

29d) 

0- 

IIL  RESULTS 


The  numerical  results  presented  here  are  for  plane- 
wave  incident  radiation  with  the  waves  traveling  along 
the  negative  *  axis.  The  detector  plane  is  the  x-y  plane. 
The  scattering  angle  0=ir— <p,  so  that  0=0*  is  the 
forward-scattered  light  and  0=180*  is  the  backscattered 
light.  We  study  in  detail  a  few  specific,  but  arbitrary, 
cases.  Since  photons  of  wavelength  A.,  are  absorbed  to 
create  photons  of  wavelength  Az,  it  is  unrealistic  to  con¬ 
sider  the  cylinder  to  be  composed  of  a  dielectric  medium. 
The  imaginary  part  of  the  complex  refractive  index 
(which  corresponds  to  the  Einstein  B  coefficient)  for  the 
cylinders  we  calculated  is  set  at  the  arbitrary  value  of 
Im(m,)  =  Im(m2)s*0.001. 

The  program  used  to  calculate  the  inelastically  scat¬ 
tered  electromagnetic  fields  requires  the  calculation  of 
the  cylindrical  Bessel  and  Hankel  functions  of  complex 
arguments.  Some  of  the  subroutines  used  to  calculate 


FIG.  3.  The  elastic  Mueller  scattering  matrix  for  a  small  liber:  a  s»0.0l5A,  m  *  l.SO+O.OOli. 
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0  Seattaring  Angle  (deg)  180  0  Scattaring  Angia  (dag) 


FIG.  4.  The  inelastic  Mueller  scattering  matrix  for  a  small  fiber:  a  =0.015A.,,  kj  —  1 . 2A,.  m,  =m2  =  1.50+0.001/. 


these  functions  are  modifications  of  those  provided  in  necessarily  apply  for  the  inelastically  scattering  case.  In 

Ref.  14,  and  all  of  them  were  checked  using  relationships  most  cases  examined,  these  criteria  were  adequate, 

and  values  given  in  Ref.  15.  Furthermore,  the  subrou¬ 
tines  were  first  used  to  calculate  the  elastically  scattered  ^  A  n  cylinder 

radiation,  which  was  compared  with  the  matrix  elements 

calculated  using  the  program  listed  in  Ref.  11.  This  pro-  It  is  instructive  to  consider  the  scattering  from  a 
gram  was  then  modified  to  calculate  the  inelastically  scat-  cylinder  of  small  kxa  (one  whose  fields  do  not  change  ap- 

tered  radiation.  preciably  throughout  its  cross  section)  since  the  number 

Calculations  of  the  scattered  electromagnetic  fields  are  of  cylindrical  harmonics  needed  to  describe  the  fields  are 

quite  involved.  The  exact  solution  [Eqs.  (20a)  and  (20b)]  greatly  reduced.  Figure  3  shows  the  elastic  light- 

requires  the  summation  of  an  infinite  number  of  terms.  scattering  Mueller  matrix  for  an  a  =0.015/.. 

Criteria  for  terminating  the  series  are  discussed  in  Ref.  11  m— 1.50  +  0.001/  cylinder.  We  note  that  the  Mueller  ma- 

for  the  elastically  scattering  case,  but  these  criteria  do  not  trix  of  this  small  cylinder  is  only  slightly  different  than 


°  Scattering  Angle  (deg)  180  0  Scattering  Angle  (dag)  180 

FIG.  S.  The  inelastic  Mueller  scattering  matrix  as  a  function  of  scattering  angle  for  three  m,  mm j~l.5-i-0.00i.  ju 

a=0.3A.,  (a),  a  =  1.0X,  (O ),  and  a  “3. 0A,  (X).  The  wavelength  A.,  =  1.2A.,. 
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tat  of  a  small,  Rayleigh  particle.  For  a  Rayleigh  parti- 
le,  matrix  element  S,2  is  zero  at  the  end  points,  and  ele¬ 
ment  Sjj  starts  at  + 100%  and  ends  at  — 100%  polanza- 
ion.  These  elements  are  useful  references  of  comparison 
or  the  inelastically  scattered  Mueller  matrix  elements. 

A  specific,  but  arbitrary,  case  of  inelastic  light  scat¬ 
tered  from  a  small  fluorescing  cylinder  is  shown  in  Fig.  4. 
For  this  case.  a=0.015A(,  A.2  =  1 . 2A|(  and  m,=m2 
=  1.50  +  0.001/.  What  is  most  apparent  in  this  set  of  ma¬ 
trix  elements  is  that  they  are  almost  constant,  with  the 
p<.  anzation  matrix  elements  S)2»  100%,  S3J=0%,  and 
Sj4  =0%.  Matrix  element  S12  being  nearly  100%  means 
that  nearly  all  the  light  is  polarized  parallel  to  the 
cylinder  axis  (TE).  For  the  elastic  scattering  of  Fig.  3, 
the  light  from  this  polarization  is  also  constant  as  a  func¬ 
tion  of  scattering  angle.  If  we  were  to  examine  just  the 
TM-polarized  light,  we  would  see  a  dip  in  intensity  near 
90*  similar  to  Rayleigh  scattering  (it  is  actually  shifted  to¬ 
ward  the  forward  scatter  slightly  in  the  inelastic  case), 
but  the  TE  mode  is  so  much  more  dominant  in  the  inelas¬ 
tic  case  that  we  are  unable  to  see  any  effects  of  the  TM 
mode  in  the  Mueller  matrix  of  Fig.  4. 

B.  Cylinder  size  parameter 

Figure  5  shows  the  light-scattering  Mueller  matrix  of 
the  m,  =  m2  =  1.5+0.001/  cylinder  having  three  different 
radii:  a=0.3A„  a  =  1.0*,,  and  a  =3.0*,.  The  fluores¬ 
cent  wavelength  is  constant  at  A2  =  1.2A,.  As  the  size  of 
the  cylinder  increases,  the  frequency  of  oscillations  in  the 
matrix  elements  increases.  This  is  similar  to  what  occurs 
for  elastic  scattering.  The  most  significant  difference  be¬ 
tween  the  elastic  and  inelastic  scattering  matrix  elements 
occurs  in  the  total  intensity:  it  does  not  necessarily  reach 
the  maximum  value  in  the  forward  scatter  (0=0*).  This 


phenomenon  is  present  in  the  scatter  from  all  three 
cylinders  of  Fig.  5. 

Of  more  interest,  perhaps,  are  the  absolute  magnitudes 
of  the  total  intensity  of  the  fluorescent  light.  For  elastic 
scattering,  intensity  increases  as  cylinder  radius  increases, 
but  for  coherent,  inelastic  scattering,  the  intensity  does 
not  increase,  even  though  the  size  parameter  has  been  in¬ 
creased  by  an  order  of  magnitude.  This  is  a  result  of  the 
coupling  that  occurs  between  the  wavelengths.  As  the 
cylinder  radius  approaches  infinity,  Eq.  (17a)  approaches 
zero,  and  A„  and  Z>„  approach  zero.  This  is  analogous  to 
taking  the  Fourier  transform  of  a  portion  of  a  sine  wave. 
For  a  very  small  portion  (consider  a  5  function),  all  fre¬ 
quencies  are  represented  approximately  equally,  but  as 
the  portion  increases,  we  see  more  of  the  fundamental 
and  less  of  the  other  frequencies,  until  only  the  funda¬ 
mental  frequency  remains  as  the  wave  traiS  extends  to 
infinity. 

C.  Flaoreacett  wavelength 

Another  parameter  which  may  be  directly  measured  is 
the  wavelength  *2,  which  is  just  one  of  the  wavelengths  in 
the  fluorescent  continuum.  Figure  6  shows  the  light- 
scattering  Mueller  matrix  of  the  m,  =m2  =  1.5 +0.001/', 
a=  A.  |  cylinder  emitting  at  three  different  wavelengths: 
A2=1.2A„  *2=1.8*.,,  and  *2=2.4*,.  It  should  be  noted 
that  the  frequency  of  oscillations  in  the  matrix  elements 
decreases  as  the  wavelength  increases.  This  is  similar  to 
what  occurs  in  elastic  scattering.  We  compare  these  ma¬ 
trix  elements  with  those  of  Fig.  7,  for  the  same  cylinder 
illuminated  by  a  *,=a/2  plane  wave.  We  note  that  the 
shapes  of  the  matrix  elements,  that  is,  the  position  and 
amplitudes  of  the  maxima  or  minima,  are  very  similar  to 


FIG.  6.  The  inelastic  Mueller  scattering  matrix  as  a  function  of  scattering  angle  from  an  mi  =  m,  =  1.3+0.001/',  a  =  A,  radius  nher 
is  measured  at  wavelength  A, =  l .  2A,  (  •  ),  A,  =  1 . 8A,  ( O  ),  and  A,  =  2. 4A,  (  X  ). 
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FIG.  7.  The  inelastic  Mueller  scattering  matrix  as  a  function  of  scattering  angle  from  an  mx  =tm2=11.5+O.OOl»,  a  —  2k,  radius 
fiber  is  measured  at  wavelength  k2  *  1 . 2a  ( * ),  k2  =  1 . 8a  ( Q  ),  and  X.2  *  2. 4c  (  X  ). 


those  of  Fig.  6.  From  this  we  conclude  that  the  shapes  of 
the  matrix  elements  are  much  more  dependent  on  param¬ 
eters  other  than  (he  incident  wavelength. 

We  can  also  study  the  scattered  light  as  a  function  of 
the  wavelength  k2  by  fixing  the  scattering  angle  and  cal¬ 
culating  the  Mueller  matrices  as  a  function  of  wavelength 
k2.  Figure  8  shows  the  light-scattering  Mueller  matrix 
elements  for  an  m  =  1.5+0.001/,  o=A.(  cylinder  as  a 
function  of  k2  at  three  scattering  angles:  0*.  90*,  and  180*. 


The  horizontal  axis  is  now  a  wavelength  scale.  The  fre¬ 
quency  of  oscillations  in  these  matrix  elements  is  greatest 
when  the  difference  between  X,  and  k2  is  small.  We  com¬ 
pare  these  elements  to  those  of  Fig.  9  for  the  same 
cylinder  system  of  Fig.  8  illuminated  by  k,-a/ 2  light. 
The  locations  and  amplitudes  of  maxima  or  minima  in 
the  matrix  elements  are  very  similar  to  those  of  Fig.  8 
even  though  the  frequency  of  the  incident  light  has 
changed  greatly. 


FIG.  8.  The  inelastic  Mueller  scattering  matrix  as  a  function  of  wavelength  kt.  The  light  scattered  Tim  an 
m,  ■,ml“l.S+O.OOW.  a  *A.,  radius  fiber  is  detected  at  angles  0*0'  (*),  $m9Cf(o ),  and  6mlMTlx )  as  wavelength  increase*  from 
12  to  2.4  times  the  fiber  radius  a. 
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FIG-  9.  The  inelastic  Mueller  scattering  matrix  as  a  function  of  wavelength  A.j.  The  light  scattered  from  an 
«!=!«:  =  1.5+0.0014  <»  “2A,  radius  fiber  is  detected  at  angles  0*(T(e),  <J=s90'(O),  and  B~  18CT  (  X  )  as  wavelength  increases  from 
1 .2  to  2.4  times  the  fiber  radius  a. 


IV.  CONCLUSION 

These  theoretical  results  demonstrate  that  the 
coherent,  inelastically  scattered  light  is  a  complicated 
function  of  many  variables.  In  addition  to  the  parame¬ 
ters  involved  in  elastic  scattering,  we  also  have  to  take 
into  account  the  wavelengths  emitted  and  the 
wavelength-dependent  optical  properties  of  the  fluoresc¬ 
ing  material.  Direct  experimental  verification  of  these 
conclusions  will  need  careful  measurements  of  the  low  in¬ 


tensity  of  radiation  because  the  total  intensity  does  not 
increase  with  cylinder  radius  as  docs  the  intensity  of  the 
elastically  scattered  light. 
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The  coherent  light  inelastically  scattered  from  an  ideal  fluorescent  sphere  illuminated  by  an  elec¬ 
tromagnetic  wave  is  analyzed  using  an  extension  of  Mie  theory.  Scattering  coefficients  are  solved 
specifically  for  plane-wave  illumination.  The  resulting  inelastically  scattered  radiation  is  similar  to  that 
of  elastically  scattered  radiation  from  a  sphere.  These  results  are  illustrated  by  numerical  calculations 
for  several  specific,  but  arbitrary,  cases. 


I.  INTRODUCTION 

Inelastically  scattered  light  is  subject  to  the  morpholo¬ 
gy  and  optical  properties  of  the  particle  that  interacts 
with  the  incident  electromagnetic  radiation.  Experimen¬ 
tally,  fluorescent  emission  from  uniform  spherical  parti¬ 
cles  was  examined  by  Biswas  et  ai.  [1]  and  Benner  et  al. 
[2].  Theoretically,  Raman  and  fluorescent  emission  by  a 
molecule  embedded  in  a  dielectric  sphere  has  been  de¬ 
rived  and  examined  [3-5].  For  a  uniform  distribution  of 
molecules  within  a  sphere,  the  coherent  process  can  be 
treated  by  adding  the  resulting  amplitudes,  molecule  by 
molecule,  over  the  cross  section  of  the  sphere.  In  another 
paper,  we  analyzed  plane-wave  illumination  of  a  uniform 
fluorescent  cylinder  by  examining  the  polarization  in¬ 
duced  at  the  new  frequency  within  the  cylinder  [6],  This 
is  the  method  used  in  this  paper. 

In  this  paper,  we  derive  and  examine  the  effects  of 
coherent,  inelastic  scattering  due  to  the  geometry  im¬ 
posed  by  a  uniform  sphere  composed  of  fluorescing  ma¬ 
terial.  The  approach  taken  in  this  paper  is  an  extension 
of  Mie's  solution  to  the  elastically  scattering  light  from  a 
uniform  sphere  [7].  The  derivation  closely  follows  the 
derivation  of  Bohren  and  Huffman  [8],  and  we  represent 
the  final  scattered  and  radiated  fields  in  terms  of  the 
Mueller  matrix  elements  [9]. 

Inelastic  scattering  occurring  in  particles  is  a  process 
that  may  be  solved  in  three  parts.  First,  the  boundary 
conditions  of  the  particle  are  used  to  determine  the  inter¬ 
nal  electromagnetic  Helds  within  the  particle.  Second, 
these  internal  Mdl  an  used  to  determine  the  fields  pro¬ 
duced  by  the  Inarming  medium.  And  third,  the  bound¬ 
ary  conditions  of  the  particle  are  used  to  determine  the 
inelastic  fields  emitted. 

The  fluorescence  process  outlined  here  is  very  similar 
to  coherent  Raman  processes.  For  incoherent  Raman 
processes,  the  internal  fields  are  not  intense  enough  to 
stimulate  the  individual  molecules  to  emit  as  they  are  for 
coherent  Raman  processes,  where  more  intense  internal 
fields  are  generated,  often  by  employing  a  second  laser. 
Although  our  theory  is  limited  as  a  model  for  Raman 
scatter  because  it  incorporates  only  an  incident  field  of 
one  frequency,  an  extension  of  it  to  multifield  interactions 
would  follow  the  same  basic  procedure.  However,  we  ex¬ 


pect  the  main  qualitative  features  of  the  scatter  generated 
with  this  simplified  model  to  remain  jinchanged. 

Knowledge  of  the  interaction  of  the  incident  electnc 
fields  and  the  fluorescing  medium  that  produces  the  emit¬ 
ted  fields  is  necessary  to  accurately  determine  the  inelasti¬ 
cally  scattered  radiation.  We  assume  that  the  polariza¬ 
tion  within  the  fluorescing  medium  is  proportional  to  the 
internal  electric  fields  due  to  the  incident  electromagnetic 
field.  We  also  assume  all  media  are  isotropic  and  homo¬ 
geneous.  With  these  assumptions,  the  derivation  is 
straightforward. 

IL  THEORY 

The  starting  point  is  to  derive  the  internal  fields  within 
the  uniform,  fluorescing  sphere  having  the  same  frequen¬ 
cy  of  the  incident  radiation.  This  is  done  by  applying  the 
boundary  conditions  of  the  sphere  on  the  electromagnetic 
field  following  the  same  procedure  used  to  derive  the 
scattered  field. 

A.  The  fields  tat  the  incident  wavelength 

Figure  1  shows  the  scattering  geometry  of  the  system. 
A  sphere  of  radius  a  is  centered  on  the  origin.  For  wave¬ 
length  A.,,  the  complex  refractive  index  of  the  sphere  is 
n , .  The  first  step  is  to  derive  the  internal  electric  fields  of 
the  sphere.  The  illuminating  fields  of  the  incident  radia- 


FIG.  1.  The  scattering  geometry  showing  the  resemg 
sphere  of  radius  a  centered  on  the  origin. 
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may  be  expanded  as  follows  (see,  for  instance,  Ref. 


E,nc,i“  2  2  +  * 


nm.  1 


N' 


n  ^0  m  *  —n 


(It 

*m,p, 


(la) 
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•1 


2  2 


8-i m,  1 


N"’ 

'’-im.p. 


(lb) 

iere  kt  =2rr/kl  is  the  wave  number,  is  the  angular 
quency,  and  p,  is  the  permeability  of  the  incident 
sdium  at  A.,.  The  vectors  ^  and  are  given 


nm,p 


=  9 


im 


-r^zX'ik.r )P  ;(cos d)e‘m* 
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nd 
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(2b) 


The  superscript  identifies  the  Bessel  functions  used  to 
represent  the  fields.  For  example,  i»  1  corresponds  with 
the  use  of  the  functions  z{„l\p)~jH{p),  i*2  corresponds 
with  the  use  of  the  functions  and  i*3  cor¬ 

responds  with  the  use  of  the  functions  zj3l(p ) 
=sH'l\p)=*j,{p)+iyl,{p).  The  functions  Pficotd)  are 
the  normalized  associated  Legendre  polynomials  defined 
by 


P  ”(CQid)=*P?(C<M$) 


(2 m  tlH»  —m) 1 
2fn  + mM 


(3) 


The  functions  representing  the  incident  radiation  are  of 
the  form  u„miOKa‘j„{ktr)P^(co»6)tmp.  The  internal 
fields  are  also  generated: 
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4b) 

where  k  \  ~kxnx  is  the  wave  number  inside  the  sphere 
medium  for  wavelength  A.,,  and  p\  is  the  permeability  of 
the  sphere  medium  at  wavelength  A.,.  The  functions 
representing  the  internal  fields  are  of  the  form 
unm  M  —j„(k  \r)P  *lcosd)e,mip.  The  scattered  fields  are 


Ek..!=  2  2  ViMi, +/«,. ,N£„.  ,  1 5a) 


n  *0  m  ^  * 


i.P\  * 


and 


=  2  2  fnm. l^C.p.+^.lN',; 


'•Pi  ' 


1 5b ) 

The  functions  representing  the  scattered  fields  are  of  the 
form  u^mtea~h[nu{kxr)P  "(cos0)e‘T 
Four  boundary  conditions  on  E  and  H  must  be 
satisfied: 

k\>t>Jkia)a^ml  +  k]^nu{kla)en„A  , 

(6a) 

P i ( * \a Km. i  “Mi < * i a Km. i  *  i «  >«„m,  i  .  (6b) 

M^-.(*i«WMI,,  =  ^ii,(*la)6,(,,1+/ii{,(l:1a)/,m,1,  (6c) 


and 


*  i  * i  <  *  'i«  i  *  *  i  K  < k  i  a  >b„w, ,  +  * £ !  k ,  a  )/,„. ,  . 


where 
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|(P>  “d  • 


B.  Ptsas-wsve  ilhunlaaftea 


When  the  coefficients  arm  l  and  b^,.,  are  known,  the 
other  field  components  of  Eq.  (6)  can  be  solved.  For  the 
case  of  plane-wave  illumination  along  the  positive  z  axis, 
polarized  along  the  x  axis,  the  coefficients  can  be  derived 
as 


®nm.l 
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The  elastic-scattering  coefficients  and  the  internal  ri.-  d 
coefficients  may  be  solved  in  terms  of  the  coefficient-  t 
the  incident  field: 


I 


Mnm,  I 


*iMi^<*'l«)jf,(*fa)-k,p'l|;(k,fl)d,(k'1a) 
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=  b  i i <3  (  i a  >  ~  fe  i Mt ^ ( fc  i  a  ^  i a  1 

k\^n(k\a)diH(k{a)~kiti\di'n(k2a)d;^k\a) 
fc  ( *  ',a  ^  ( /c  t  a )  -  A:  t  Ar  t  a  l  A:  ',a )  ’ 
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C.  Tbe  fluorescent  medium 


and 


The  absorption  of  light  in  the  fluorescent  medium  is 
proportional  to  the  magnitude  of  the  electromagnetic 
field.  We  agree  with  the  assumption  of  Ref.  [3],  which 
states  that  the  transmitted  electric  field  within  the  sphere 
induces  a  polarization  within  the  medium  oscillating  at  a 
shifted  frequency,  w,;  i.e„  the  distribution  of  excited 
fluorescent  molecules  emitting  a  particular  wavelength  k2 
is  given  as  follows: 

P(r,0,<p,r)=saEim.,(/\0,^)e  '“I<  ,  (10) 

where  a  is  an  excitation  characteristic  of  the  fluorescent 
medium  (in  the  case  of  Raman  scatter,  a  corresponds  to  a 
gain  term).  We  consider  the  coherent  field  produced  hav- 
ing  the  specific  characteristic  wavelength,  angular  fre¬ 
quency,  permeability,  and  wave  number  in  the  incident 
medium  of  A.2,  <o2,  n2,  and  k2,  respectively.  The  refractive 
index,  wave  number,  and  permeability  of  the  sphere  at  A.2 
are  n2,  k  j,  and  n'z,  respectively.  We  first  expand  the  field 
oscillating  at  a2  into  its  components: 
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In  order  for  the  fields  to  behave  as  expanding  spherical 
waves,  the  functions  of  the  fields  created  by  fluorescence 
are  of  the  form  u„m  a=hj,l,(k'2r)P  ”(cosff)e‘mv.  The 
source  of  these  fields  is  the  electric  polarization  expressed 
in  Eq.  (10).  The  coefficients  anm  2  and  bnm  2  can  be  deter¬ 
mined  by  expanding  the  electric  polarization  in  terms  of 
the  vector  harmonics  of  generating  function  u„m „.  Since 
the  electric  polarization  is  composed  of  vector  harmonics 
of  the  generating  function  unm l  jnt,  we  can  start  by  expand¬ 
ing  the  internal  electric-field  vector  harmonics  in  terms  of 
the  electric-field  vector  harmonics  oscillating  at  a>2: 


*C,r  2 


»  m  ,j»j  "  m 


N( 


O) 


1 1 2a) 


and 


nm.p, 


A 

n‘,m’ 


n'm\ 


12b) 


The  following  orthogonal  relationships  exist  for  the  vec¬ 
tor  harmonics  [10]-. 
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By  scalar  multiplying  (12a)  and  (12b)  by  (Af^’p-  )*  and  integrating  both  equat  ons  over  the  fluorescing  sphere,  then  re¬ 
peating  the  process  with  ( N'J '  .  )*,  we  can  solve  for  the  following  coefficients: 
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/whN'^V,'(N^.p;>*</r  _  f‘[(n  +\'ljn-\(k\r)jn-i{k'lr)  +  ni„  +  l{,k\r)jn  +  l(k'ir)}r1dT 
Dnm~  fj*~*{*^*dT  /0°C( "  1  - » ( *  ^  i ( ^ ,( /c^)]rJ4/r  * 

fl„=C„=  0. 

The  integrals  in  equations  (l 4a)  and  (14b)  are  Lommel’s  integrals  and  have  the  following  solutions  [11]: 

fajn[k:r)jn^k\r]r1dr  =  ——r—-^[kj'H{k’1a)j,,{k\a)-k\j'n(k\a)j„{k'1a)) 
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The  coefficients  of  the  emitted  electromagnetic  fields  may 
be  written  as 


— ^  ^  A 

*2 
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.  (16b) 

e2 

where  ej  is  the  cylinder  permittivity  at  the  emitted  fre¬ 
quency.  An  induced  field  within  the  sphere  is  needed  to 
satisfy  the  boundary  conditions.  This  field  corresponds 
to  the  field  created  by  fluorescence  reflecting  off  the  linear 
surface  of  the  sphere.  Using  functions  of  the  form 
u„m  ind  -jH(k\r)P ”(coidU‘m*,  the  induced  field  may  be 
expre  *  -*s 

E,*U-  i  i  .  (17.) 

*  -o  m  m  -n  *  1 


The  generating  functions  of  the  scattered  fields  are  of  the 
form  hin(fc2r)^  "(coi8)t,m*. 

D.  Boaadary  coodhhws 
at  the  wavelength  of  flaoraeceace 

Four  boundary  conditions  on  E  and  H  must  be 
satisfied: 
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Finally,  the  scattered  laUa  may  be  expressed  as 
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The  scattering  amplitude  coefficienu  are 
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With  analytical  expressions  for  these  coefficients,  the 
solution  to  the  scattered  field  from  a  coherently  fluoresc¬ 
ing  sphere  is  complete. 

E.  Scattered  light  for  plane-wave  illumination 

We  now  look  at  the  specific  case  of  light  scattered  into 
the  far  field  by  a  fluorescing  sphere  illuminated  by  a  plane 
wave.  For  large  arguments,  the  first-order  Hankel  func¬ 
tion  is  given  asymptotically  by  the  following: 

i  —i\ngiP 

K"ipy~  .  — ■  1pi»  i-  <2i) 

We  examine  the  light  in  two  separate  polarization  states: 
the  TE  state  corresponds  with  light  polarized  perpendicu¬ 
lar  to  the  scattering  plane  and  the  TM  state  corresponds 
with  light  polarized  parallel  to  the  scattering  plane.  The 
amplitude  scattering  matrix  is  given  in  the  far  field  by 

£tE.k«  eiklr  S,  S*  ^rtinc 

£tM.«c«  —I  k2r  1^3  (^TM.inc 

where 

*  P '  a 
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(23a) 

-  Pl  a 

S:=  2-2 !(-.)"  ~^enU^-^P\^)U :J  > 

(23b) 

and 
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The  four  independent,  nonzero 
may  be  expressed  by 

Mueller  matrix  dements 

sll=;(!sl!2+!s2i2) , 

2-lai 

Su  =  t(  S2i2-'S,i2)/Sn  - 

24b- 

S33  =  Re(S1Sj  )/Sn  , 

-24c) 

and 

S34  =  Im(S2Sp/S,,  . 

-24d) 

III.  RESULTS 

The  numerical  results  presented  here  are  for  plane- 
wave  incident  radiation.  We  study  in  detail  a  few 
specific,  but  arbitrary,  cases.  Sincf  photons  of  wave¬ 
length  A.,  are  absorbed  to  create  photons  of  wavelength. 
A.2,  it  is  unrealistic  to  consider  the  sphere  to  be  composed 
of  a  dielectric  medium.  The  imaginary  pan  of  the  com¬ 
plex  refractive  index  (which  corresponds  to  the  Einstein 
B  coefficient)  for  the  spheres  is  set  at  the  arbitrary  value 
of  Imlfi)  )  =  Im(/i2  )=0.001. 

A.  A  laall  sphere 

It  is  instructive  to  consider  the  scattering  from  a 
sphere  of  small  It, a  (one  whose  fields  do  not  change  ap¬ 
preciably  throughout  its  cross  section)  since  the  number 
of  spherical  harmonics  needed  to  describe  the  fields  is 
greatly  reduced.  Figure  2  shows  the  elastic  light¬ 
scattering  Mueller  matrix  is  an  a  =0.00  Urn 
n  =  1.50+0.  OOli  sphere.  For  such  a  small  sphere  <a  Ray¬ 
leigh  panicle),  matrix  element  Sn  is  proponional  to 
l+cos20,  matrix  element  S,2  is  equal  to 

— sin20/(  1  +cos20),  element  S33  is  equal  to 


(22) 


FIG.  2.  The  elastic  Mueller  scattering  matrix  for  a  small  sphere:  a  *0.0011.,  n  m  1. 50+0.001/.  /Ia)  represents  total  mten-.ity 
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0  Scstcsnng  anglt  Uag)  ISO  0  Scattering  angle  (dag)  ISO 

FIG.  3.  The  inelastic  Mueller  scattering  nutria  for  a  small  sphere:  a  =0.001A.,,  1.2A.,,  n,  =*n2  =  1. 50+0.0011. 

total  intensity. 


represents 


2  cos<?/(  l  +cos20),  and  element  SH  is  zero  everywhere. 
These  elements  are  useful  references  of  comparison  for 
the  inelastically  scattered  Mueller  matrix  elemenu. 

A  specific,  but  arbitrary,  case  of  inelastic  light  scat* 
tered  from  a  small  fluorescing  sphere  is  shown  in  Fig.  3. 
For  this  case,  a  =0.001A.,,  L2*1.2&|,  and 

n,  =n2  —  1.50+0.001i.  The  shapes  of  the  matrix  ele¬ 
ments  of  Fig.  3  are  identical  to  the  shapes  of  the  matrix 
elements  for  the  small,  inelastic  sphere  of  Fig.  2.  The 


only  difference  between  these  two  sets  of  curves  is  the 
magnitude  of  matrix  element  5,,  (note  the  total  scatter¬ 
ing  has  decreased  by  nine  orders  of  magnitude). 

Figure  4  shows  the  light-scattering  Mueller  matnx  of 
three  different  «,  =cn2:=  15+0.001i  spheres  having  radii 
a  S0.3L|,  l.0A.|,  and  3.0A.,.  The  fluorescent  wavelength 


a 


a 


-3 


FIG.  4.  The  inelastic  Mueller  scattering  matrix  as  a  function  of  sphere  radius.  Curves  are  shown  for  three  (t,”n2“15-0'WMi 
sphere  radii:  a  m0.3K,  (  •  —  •),  a -l.OL,  (O—  O ),  and  a -3.0A.,  <  X  —  x  ).  The  wavelength  1.2A.,.  /,„  represents  total  intensi¬ 


ty- 
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2 


O 


-1 


FIG.  5.  The  inelastic  Mueller  scattering  matrix  as  a  function  of  wavelength  4  The  light  scattered  from  an  n,  =n3  - 1.5—0. 00 1 
a=A.,  radius  sphere  is  measured  at  wavelength  A.j~1.2A,  (•  —  «),  A2=  1.8A,  <o— O ),  and  A.2»2.4A,  (X  — X).  represents  the 
total  intensity. 


is  set  at  Aj*  I.2A,.  As  the  size  of  the  sphere  increases, 
the  frequency  of  oscillations  in  the  matrix  elements  in* 
creases.  This  is  similar  to  what  occurs  for  elastic  scatter* 
ing.  One  major  difference  between  the  elastic  and  inelas¬ 
tic  cases  is  illustrated  in  matrix  element  Sn.  For  elastic 
scattering,  the  total  intensity  of  the  scattered  radiation 
increases  significantly  as  the  sphere  radius  increases. 
Such  an  increase  does  not  necessarily  occur  for  the  case 
for  coherent,  inelastic  scattering.  Figure  4  shows  that  as 
the  sphere  radius  is  increased  from  0.3A,  to  1.0A,,  the  to¬ 
tal  intensity  has  increased  at  all  scattering  angles,  but  as 


the  sphere  radius  is  increased  from  1.0A,  to  3.0A,,  the  to¬ 
tal  intensity  has  decreased  at  all  scattering  angles. 

C  Flaorasceat  wavelength 

Another  —Ac  that  may  be  directly  measured  is 
the  wavelength  AjrWhich  is  a  single  wavelength  selected 
from  the  fludrtscflk  continuum.  Figure  S  shows  the 
light-scattering  ’^Mueller  matrix  of  the 
1. 3  a -A,  sphere  emitting  at  three 

different  wayrjlht^:  A2®1.2A„  A2  =  1.8A,,  and 


-tool 


12  t  aavtlenetfl  2 A 


FIG.  6.  The  inelastic  Mueller  scattering  matrix  u  a  function  of  wavelength  A2.  The  light  scattered  from  an  n ,  =  n  -  '00!/, 
a-A,  radius  sphere  is  detected  it  angles  fl-CT  (•  —  •),  S-W  (O— O  I,  and  6- 180*  <  X  —  X )  as  wavelength  A.,  mere  • .. .  •  ,m  I  2 
to  2.4  times  the  sphere  radius  a. 
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FIG.  7.  The  Mueller  matrix  tor  elastically  ( • — * )  and  inelasticaily  (O—O  )  scattered  light  from  an  nt  “rij  =  1.3+0.001/,  a  =  A, 
radius  sphere.  The  inelasticaily  scattered  light  is  detected  at  A.j*  I.OOIA.,.  /,M  represents  the  total  intensity. 


X2  —  2. 4A.|.  It  should  be  noted  that  the  frequency  of  oscil¬ 
lations  in  the  matrix  elements  decreases  as  the  wave¬ 
length  increases.  This  is  similar  to  what  occurs  in  elastic 
scattering. 

We  can  also  study  the  scattered  light  as  a  function  of 
the  wavelength  A.z  by  fixing  the  scattering  angle  and  cal¬ 
culating  the  Mueller  matrices  as  a  function  of  wavelength 
k:  Figure  6  shows  the  light-scattering  Mueller  matrix 
elements  for  an  n,  =n2  =  1.5+0. OOli,  a  sphere  as  a 
function  of  k2  at  three  scattering  angles:  O’,  90*.  and  180*. 
The  horizontal  axis  is  now  a  wavelength  scale  varying 
from  1.2  to  2.4  times  the  sphere  radius.  Some  features 
occur  in  these  curves  that  occur  in  elastic  scattering. 
Matrix  elements  SI2  and  5M  are  zero  at  0*  and  180*, 
while  matrix  element  S}1  is  100%  at  O’  and  — 100%  at 


180*. 


IV.  CONCL1 


To  conclude,  we  show  how  iha's^^dBPT scattered  ra¬ 
diation  differs  fro*  the  ineleifinalMirsfrernil  radiation. 
Figure  7  compares  the  Moeller  ma^p  fhre&tsricaJly  and 
inelasticaily  scattered  light  from  an  n,  =*n2 *  1.5 
+•0.001/,  a  “A.,  radios  sphere.  The  inelasticaily  scat¬ 
tered  light  is  detected  at  1.001A.,.  To  the  experimen¬ 


talist,  this  is  equivalent  to  illuminating  the  sphere  at  400 
nm  and  comparing  the  scatter  at  400  nm  (which  is  almost 
completely  elastic)  with  the  scatterer  at  400.4  nm  (which 
is  completely  inelastic).  We  note  that  the  shapes  of  the 
two  sets  of  elements  are  very  similar,  with  the  main 
difference  being  the  amplitudes  of  the  peaks.  This 
difference  is  due  to  the  difference  between  elastic-  and 
inelastic-scattered  light,  not  the  small  wavelength 
difference. 

These  theoretical  results  demonstrate  that  the 
coherent,  inelasticaily  scattered  light  is  a  complicated 
function  of  many  variables.  In  addition  to  the  parame¬ 
ters  involved  in  elastic  scattering,  the  wavelengths  emit¬ 
ted  and  the  wavelength-dependent  optical  properties  of 
the  fluorescing  material  must  be  considered.  Direct  ex¬ 
perimental  verification  of  these  conclusions  will  need 
careful  measurements  of  the  tow  intensity  of  inelastic  ra¬ 
diation  because  the  total  intensity  does  not  increase  with 
sphere  radius  as  does  the  intensity  of  the  elastically  scat¬ 
tered  light. 
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Light-scattering  Mueller  matrix  from  a  fiber  as  a 
function  of  MgO  contamination 

Gordon  VkJeen  and  William  S.  Bickei 


The  li<ht-»caltermf  Mueiier  matrix  for  in  r  «  0.345-wm-radiux  quarts  fiber,  iiiuminated  at  A  ■  0.44 16  «m.  is 
examined  aa  a  function  of  contamination  with  MfO  crystals.  When  the  M|0  contamination  is  low.  the  matrix 
elements  resemble  thoee  of  a  fiber  of  slightly  larger  radius.  The  MgO  contamination  creates  higher- 
frequency.  smaller-ampiitude  oscillations  in  the  matrix  elements  that  mask  the  lower-frequency  oscillations 
indicative  of  a  perfect  cylinder.  The  contamination  also  causes  scatter  outside  the  plane  of  incidence. 

Key  words:  Mueller  matrix,  light  scattering,  fiber,  contamination. 


If  iwOOUvwflfl 

The  problem  of  light  scattering  from  an  infinitely  long 
circular  cylinder  was  solved  independently  by  Lord 
Rayleigh1  and  by  von  Ignatowsky.2  Theory  was  ex¬ 
perimentally  verified  by  Bell  and  Bickei,3  who  mea- 
sured  the  light-scattaring  Muellar  matrix  from  a 
quartz  fiber  whoee  radius  was  approximately  equal  to 
the  wavelength  of  the  illuminating  radiation.  Certain 
modifications  to  the  basic  cylindrical  system,  such  as 
cladding  with  another  index  material  or  giving  the 
system  an  elliptical  cross  section,  can  be  treated  theo¬ 
retically.  These  solutions  have  been  discussed  in  sev¬ 
eral  texts.4-6  When  the  scattering  system  geometry 
becomes  irregular,  that  is,  whan  structures  or  surfaces 
cannot  be  made  to  conform  to  a  single  orthogonal 
coordinate  system,  the  theoretical  solution  becomes 
extremely  difficult,  if  not  impoaaibie,  to  obtain.  Theo¬ 
retical  progress  toward  solving  for  the  light  scattered 
from  such  systems  naada  experimental  data  so  the 
theorist  can  be  guided  toward  the  (implications  that 
may  be  made  while  still  yielding  acceptable  results. 

In  studying  light  scattering  it  is  important  to  consid¬ 
er  that  no  scattarfaf  system  is  perfect  Impurities  in 
the  material  andmtortteaa  in  geometry  cause  experi¬ 
mental  measurements  of  scattered  light  to  deviate 
from  theory.  Many  systems,  by  their  nature,  are  im¬ 
perfect.  The  addition  of  contaminants  to  the  scatter¬ 
ing  system  (dust  particles  in  air,  minerals  in  water. 
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etc.)  change  the  systems  and  their  resulting  scatter.  fc 
Contamination  monitoring  is  an  especially  difficult  * 
problem  when  the  experimenter  is  not  able  to  separata  ■- 
the  sample  to  be  monitored  from  other  scatterere  that  i 
might  be  present.  For  example,  scattering  signals  W 
from  asbestos  impurities  could  be  dependent  on  the 
dust  present  or  on  humidity  levels.  The  monitor  can¬ 
not  remove  or  ignore  these  other  factors  and  so  must 
deal  with  their  contributions  to  the  scattering  signals. 

In  the  laboratory,  greater  control  over  the  scattering 
system  can  be  achieved. 

The  light-scattering  Mueller  matrix  for  a  circular- 
cross-section,  quartz  (n  ■  1.466  -  O.Oi  at  X  -  0.4416 
Mm)  fiber  of  known  optical  and  geometrical  parameters 
is  predicted  exactly  by  theory.  This  same  fiber,  coated 
by  cubic  MgO  crystals  (n  ■  1.74  -  O.Oi  at  A  ■  0.4416 
Mm),  has  a  different  light-scattering  Mueller  matrix. 

To  see  how  the  contamination  affected  the  scattering, 
we  studied  the  light  scattering  as  a  function  of  MgO 
contamination  of  the  fiber.  By  studying  how  a  perfect 
quartz  fiber’s  scatter  changes  when  we  contaminate  it 
with  particles  of  known  characteristics,  we  may  better 
understand  how  contaminants  affect  other,  irregular, 
systems. 

R.  UgM-Scattering  Technique 
The  polar  nephelometer  used  in  this  study  employs  the 
polarization  modulation  technique  developed  by  Hunt 
and  Huffinan.7  A  complete  discussion  of  the  nephe- 
^lometer  design  and  operation  is  given  by  them  and  also 
by  Perry  et  of.*  and  Bickei  et  of.*  This  technique 
involves  periodically  modulating  the  incident  beam's 
polarization  stats  at  wo  »  50  kHz,  using  a  photoelastic 
modulator,  while  observing  the  signals  carried  out  by 
the  fundamental  frequency  (wo)  and  the  second  har¬ 
monic  (2wo)  of  the  scattered  light  A  lock-in  amplifier 
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fig.  1.  Experimental  apparatus  used  to  mtaaure  the  Mueller  scat¬ 
tering  matrix  'f  a  fiber.  Laser  light  paaaea  through  the  entrance 
optics  and  strikes  the  fiber  at  normal  incidence.  The  scattered 
radiation  pa:set  through  the  exit  optics,  where  it  can  be  intercepted 
by  a  detector  or  a  photographic  plate  to  make  a  photograph. 


is  used  to  demodulate  the  signals.  The  experimental 
scattering  matrix  StJ{9)  can  be  measured  with  the  prop¬ 
er  choice  of  exit  optics.10  Normalization  of  the  matrix 
elements  is  performed  by  servoing  the  photomultiplier 
tube  (RCA  1P21)  gain  by  means  of  a  constant-current 
servo,  which  creates  a  constant  dc  output  signal  over 
the  entire  scan  of  the  detector.  The  normalized  out¬ 
put  S,,(0)  analog  signals  are  collected  and  sent  to  a 
computer. 

Figure  1  shows  the  experimental  apparatus  used  to 
measure  the  S,,(0)  of  a  fiber.  The  beam  from  the  laser 
passes  through  the  entrance  optics  and  illuminates  the 
fiber  at  normal  incidence.  The  light  scattered  by  the 
fiber  may  then  either  pass  through  various  exit  optics 
to  a  detector  or  pass  directly  to  a  photographic  plate. 
A  photographic  plate  can  record  the  out-of-plane  scat¬ 
tering,  which  does  not  occur  for  a  perfect,  uncontami¬ 
nated  fiber  and  which  is  usually  missed  by  the  polar 
nephelometer. 

M.  Initial  System 

The  quartz  fiber  was  made  with  the  method  developed 
by  Bell  and  BickeL3  The  midsection  of  a  quartz  rod 
was  heated  by  an  oxygen-acetylene  flame  until  it  was 
molten.  The  rod  eras  stretched  slightly  and  slowly  to 
make  the  midsection  thinner.  Just  before  the  rod 
separated  into  two  pieces,  the  ends  were  rapidly  pulled 


Fig.  2.  Four  experimentally  mtaaurtd  MuaUtr  Maturing  matrix  tlemtnu  for  tha  initial  fibtf  ftotud  curvas)  and  tht  theoretical  matrix 

tlemanu  for  an  r  •  0.34S-»m-radiua  quarts  fiber  (xoiid  curvet). 
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apart,  This  maneuver  often  produce#  fibers  stretch¬ 
ing  from  one  or  both  enda  of  the  rod.  If  no  fibers  are 
created,  the  process  can  be  repeated. 

Figure  2  compares  the  four  unique,  nonzero  light- 
scattering  Mueller  matrix  elements  for  the  experimen¬ 
tal  fiber  (dotted  carves)  with  the  theoretical  matrix 
elements  of  an  r  •  0.345-vm-radius  quartz  fiber  (solid 
curves).  The  elements  are  nearly  identical  with  differ¬ 
ences  most  likely  resulting  from  a  slight  nonuniformity 
of  the  fiber  cross  section.  These  matrix  elements  are 
the  scattering  signatures  of  the  perfect  fiber  system. 
These  elements  will  change  as  contaminants  are  added 
to  the  system. 

IV.  Contaminated  System 

Figure  3  shows  the  apparatus  used  for  contaminating 
the  -quartz  fiber  with  MgO  crystals.  The  fiber  is  sus¬ 
pended  vertically  from  a  tube,  22  cm  long  and  8  cm  in 
diameter,  the  top  of  which  is  blocked  off  except  for  a  2 
cm  x  2  cm  hole  on  one  end  opposite  the  fiber.  The 
fiber  is  coated  by  burning  a  2.3-cm-long  strip  (~24  mg) 
of  magnesium  ribbon  ( J.  T.  Baker  Chemical  Company) 
placed  2  lui  below  the  lower  end  of  the  tube.  The 
coating  occurs  as  the  smoke  passes  through  the  tube 
and  out  through  the  hole  at  the  top.  The  ribbon  burns 
for  *-10  s,  and  after  the  smoke  hss  cleared  (~2  min),  the 
contaminated  fiber  is  removed.  The  fiber  is  placed 
under  a  microscope  and  photographed.  Then  it  is 
placed  in  the  nephelometer,  where  its  scattering  ma- 


MgO  Smoto 


[.3.  Apparatuaueed  to  contaminate  the  quartz  fibers  with  MgO 
crystals. 


Scattering  Angle  Scattering  Angie 

Fig.  4.  Experimental  light -*cettefinf  Mueller  matrix  of  an  r  »  0.345-mn-radiua  quarts  fiber  aa  a  function  of  MgO  contamination  The 
uncoated  fiber  a  la  menu  are  depicted  by  dote  on  aolid  curves  The  other  curvaa  are  for  the  fiber  coated  once  end  twice  with  .MgO  crvstale. 
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Fi*.  5.  Tbaorttical  lifht-tcaturin*  MudUr  maths  •lamanta  for  a  quarts  fibar  at  a  function  of  fibtr  radiua.  Tha  alcmanta  for  the  initial  r  • 
0.34S-*m-radiua  quarts  fiber  are  thoem  by  dou  on  solid  curves.  Tha  other  curve*  are  for  fiber  radii  of  0.350,  0.355.  and  0.360  «a 


trix  elements  are  measured.  A  photographic  plate  ia 
used  to  measure  the  total  intensity  of  the  out-of-plane 
scatter.  The  process  is  then  repeated,  each  time  with 
additional  coatings  of  contaminants  added,  until  the 
light  scattered  is  no  longer  significantly  affected  by  the 
additional  contaminants. 

Figure  4  shows  the  light-scattering  Mueller  matrix 
of  the  riber  fur  three  different  levels  of  MgO  contami¬ 
nation.  The  light-scattering  total  intensity  matrix 
element  (Su)  increases  with  increasing  fiber  contami¬ 
nation.  while  the  magnitudes  of  the  maxima/ minima  of 
the  polarization  matrix  elements  (Si2,  Sjj,  and  Sm) 
tend  to  decrease  as  the  fiber  becomes  more  contami¬ 
nated.  The  frequency  of  tha  oscillations  increases  as 
the  contamination  increaaaa  The  matrix  elements 
also  display  high-frequency,  small-amplitude  ostilla- 
tiona  These  small  (only  a  few  percent)  signal  fluctua¬ 
tions  are  not  noise  einse  they  are  reproduced  exactly  by 
repetitive  meeaureaMMs  of  tha  fiber  in  exactly  the 
same  orientatioo.  They  change,  however,  if  the  laser 
beam  strikes  tha  fiber  at  a  slightly  different  location  or 
orientation. 

Figure  $  shows  tha  theoretical  matrix  elements  for  a 
quarts  fiber  as  a  function  of  fibar  radiua.  A0.345-um- 
radiua  fiber  (dote  on  solid  curves)  is  shown;  additional 
sets  of  fitments  art  shown  for  fibers  whose  radii  have 
increased  by  0.006  wm.  Figures  4  and  5  show  many 
similarities.  The  low-frequency  oecillation  of  the  con- 


Fif.  6.  Photomicrograph  of  tbs  fiber  stiffen  system  after  lix  coat- 
in**  of  M*0  cryaula. 
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crystals  (dotted  currw). 


Scmmi  instil  (<M«> 


Fig.  a  Photograph  of  th*  light  scatund  by  the  fiber  with  tis 
coetinfs  of  MgO  ciyilsh 


Laminated  fiber  coincidta  with  that  of  a  liber  of  slightly 
larger  radius.  Not  only  have  the  frequencies  increased 
but  the  amplitudes  of  the  maximr /minima  of  the  con¬ 
taminated  fiber  matrix  elements  coincide  with  the  ele¬ 
ments  of  the  larger  fibers.  Note  especially  the  shapes 
of  ths  elements  near  9  “  130*. 

Figure  8,  s  photomicrograph  of  ths  fiber  system, 


shows  that  after  six  coatings  the  MgO  crystals  extend 
many  fiber  diameters  beyond  ths  initial  fiber  bound¬ 
ary.  Figure  7  shows  the  experimental  light-scattering 
Mueller  matrix  of  ths  uncoated  fiber  and  of  the  fiber 
after  six  costings  of  MgO  crystals.  The  phase  infor¬ 
mation  that  characterizes  the  initial  uncoated  fiber  is 
no  longer  present  for  the  coeted  fiber.  Instead,  high- 
frequency  oscillations  dominate  ths  matrix  elements. 
As  ths  number  of  contaminants  increases,  matrix  ele¬ 
ments  Six  and  Sj4  tend  toward  zero.  Matrix  element 
S33  also  tends  toward  zero,  except  in  the  forward- 
scattering  region  (near  9  ■  0*),  where  it  is  positively 
polarized. 

The  Mueller  scattering  matrix  elements  for  a  single 
sphere  or  cylinder  obey  the  following  relationship: 

S|/  +  Sjj*  ♦  S],J  •  1 

(note  that  we  are  using  a  normalized  Mueller  matrix 
representation).  This  relationship  is  obeyed  by  the 
uncoated  fiber  at  all  scattering  angles.  However,  this 
relationship  doss  not  hold  for  ths  contaminated  fiber, 
which  has  lost  its  efficiency  in  producing  polarized 
light 

F;gure  8  is  a  photograph  of  ths  scattering  patterns 
created  by  ths  contaminated  fiber  shown  in  Fig.  6.  A 
sheet  of  Kodak  Panalure  H  photographic  paper  was 
placed  lm  sway  from  ths  fiber.  Both  the  fiber  and  the 
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photographic  paper  were  normal  to  the  incident  beam. 
It  is  intereating  that  the  pattern  ia  especially  intense  in 
the  forward  scatter  very  near  specular  (0  -  0*.  8  »  0#) 
Fiber  asymmetry  along  the  z  axis  causes  out-of-plane 
scatter.  This  out-of-plane  scatter  is  especially  inter¬ 
esting  because  it  is  not  a  random  speckle  pattern.  The 
speckles,  elongated  in  the  9  direction,  betray  the  un¬ 
derlying  cylindrical  geometry.  Although  this  is  not 
evident  in  the  matrix  elements  measured  in  the  scat¬ 
tering  plane  (Fig.  7),  it  is  evident  in  out-of-piane  scat¬ 
ter. 


V.  otscueekm 

Adding  contaminants  to  a  perfect  fiber  changes  its 
light-scattering  Mueller  matrix.  Small  amounts  of 
contaminants  tend  to  increase  fiber  size  slightly,  and 
this  information  appears  in  the  matrix  elements.  Ad¬ 
ditional  contamination  drives  the  polarization  ele¬ 
ments  (with  the  exception  of  S33  in  the  forward -scatter 
region)  toward  zero.  Only  the  out-of-plane  scatter 
contains  a  hint  of  a  fiber  system. 

This  research  was  supported  in  part  by  the  U.S.  Air 
Force  Office  of  Scientific  Research. 
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The  light-scattering  Mueller  matrix  i*  experimentally  determined  for  a  0.26-um-radiua  quartz  fiber  mounted 
on  an  aluminum  surface  at  five  different  incident  angles.  The  results  are  compared  with  those  for  the  experi¬ 
mental  scattering  elements  of  the  lone  fiber  and  the  tone  surface  and  with  theoretical  results  derived  from  a 
simple  fiber-surface  model.  The  experimental  matrix  elements  of  the  fiber-surface  system  do  nA resemble 
the  matrix  elements  from  any  of  these  other  systems. 


INTRODUCTION 

A  perfect  fiber  (or  sphere)  on  a  perfect  surface  is  the  sim¬ 
plest  solvable  surface  defect.  We  measured  the  Mueller 
scattering  matrix  for  a  fiber-surface  system  with  a  polar 
nephelometer.  This  scattering  system  incorporates  two 
separate  perfect  systems  whose  scatter  can  be  determined 
exactly  by  theory.  Fresnel  derived  the  scatter  from  a  per¬ 
fect  plane  surface  before  Maxwell  developed  his  famous 
equations.  Later  in  the  19th  century,  Rayleigh  derived 
the  scatter  from  a  perfect  cylinder.1  Although  the  combi¬ 
nation  of  a  fiber-surface  system  is  more  difficult  to  solve 
theoretically,  we  can  expen  mentally  measure  the  scatter 
of  the  combined  system  just  aa  easily  as  we  can  measure 
the  scatter  of  the  individual  systems.  The  main  experi¬ 
mental  problem  is  to  ensure  that  the  surface,  the  fiber, 
and  the  fiber-surface  system  are  of  sufficiently  high  qual¬ 
ity  to  give  lata  accurate  enough  to  warrant  careful  theo¬ 
retical  attention. 


INSTRUMENTATION 

The  light-scattering  apparatus  used  to  msaeure  the  light- 
scattering  Mueller  matrix  eiemeata  ia  shown  schemati¬ 
cally  in  Fig.  1.  Light  emitted  from  a  He-Cd  laser  (A  » 
0  4416  uni)  passes  through  various  entrance  optics  before 
striking  the  sample  (either  the  fiber,  the  surface,  or  the 
fiber-surface  system).  Elements  were  measured  for  the 
lone  fiber  oriented  along  the  t  axis,  perpendicular  to 
the  incident  beam,  and  for  the  lone  eurfacs.  with  the  plane 
of  the  surface  defined  by  the  z  axis  and  a  line  at  angle  a 
measured  from  the  incident  beam.  The  angle  a  is  the 
complement  to  the  incident  angle  measured  from  the  sur¬ 
face  normal  (see  Fig.  1).  Elements  were  then  measured 
for  the  fiber-surface  system  with  the  fiber  resting  on  the 
surface.  The  fiber  is  oriented  parallel  to  the  z  axis,  which 
is  perpendicular  to  both  the  incident  beam  and  the  normal 
to  the  surface.  Light  from  the  scattering  system  pasees 
.h rough  various  exit  optics  before  reaching  the  detector, 
which  rotates  through  angle  0  about  the  z  axis. 

The  polar  nephelometer  used  in  this  study  tmploys  the 
polarization  modulation  technique  developed  by  Hunt  and 


Huffman.2  Complete  discussions  of  the  nephelometer  de¬ 
sign  and  operation  are  given  in  Ref.  2  and  also  by  Perry 
et  of.2  and  Bickal  ft  ai. **  The  experimental  scattering 
matrix  measured  by  the  nephelometer  is  signified  by  S„*. . 
These  elements  are  actually  combinations  of  the  general 
scattering  matrix  elements  S„.  Iafelice  end  Bickei*  de- 1 
termined  that  five  elements  are  neceesary  to  characterize^ 
the  Mueller  matrix  for  a  near-perfect  surface:  j 

Su*  -  S,i  S||*  ■  Su/Su, 

Sts*  m  (Su  +  Sa)/(S„  +  S«), 

Sjj*  “  (S„  +  Sjj)/(Su  +  Sji), 

Sm*  m  (Si*  +■  SM)/(Su  +  Sn).  (1) 

For  symmetric  scatterers  such  as  spheres  and  fibers, 
Sn  *  Sn  »  Si*  *  0  and  Sn  »  S«.  For  these  scatterers, 
the  experimental  scattaring  matrix  reduces  to  the  normal¬ 
ized  scattering  matrix: 

Sn*  m  Sn ,  Su*  »  Su/Su,  Sjj*  *  1 , 

Sjj*  *  Sjj/Su,  S,**  M  Su/Su  ■  1 2) 

SCATTERING  SYSTEMS 

The  quartz  fiber  waa  made  by  using  the  method  developed 
by  Bell  and  Bickal.7  The  radius  of  the  fiber  was  deter¬ 
mined  by  comparing  the  values  of  the  experimentally  mea¬ 
sured  matrix  elements  with  theoretical  values.9'10  A  beet 
fit  of  the  experimental  and  theoretical  curves  occurred  for 
a  fiber  radius  of  0.260  £  0.006  Mm.  The  optical  con¬ 
stants  of  the  fiber  were  held  constant  (for  quartz,  the  re¬ 
fractive  index  a  ■  1.466  ♦  0.01  end  permeability  m  3  mo). 
The  experimental  matrix  elements  for  the  quartz  fiber 
(dotted  curves)  end  the  theoretical  matrix  elements  for 
a  0.26-Mm-radiua  quartz  fiber  (solid  curves)  are  shown 
in  Fig.  2.  Repeated  measurements  of  the  fiber  matrix 
elements  coincide.  Deviations  occur  because  the  fiber  is 
not  perfect. 

The  surface  studied  in  this  experiment  is  the  same  sur¬ 
face  as  that  studied  in  great  detail  by  Iafelice  and  Bickei.* 
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Fig.  7.  The  four  unique  Mueller  matrix  element*  for  •  fiber - 
surface  system  illuminated  at  a  »  78.75*  from  an  aluminum 
surface. 


0.26-m~  radius  quart*  fiber:  Su*.  ^u*.  Sb*.  and  S*m. 
Meaaurements  of  the  other  matrix  element*  show  that  the 
relation*  expressed  by  Eq».  (2)  are  valid  for  the  fiber- 
surface  system  and  that  only  these  four  measurement*  are 
needed  to  characterize  the  scatter  completely.  Figures  4, 
5, 6.  and  7  show  'he  fiber-surface  matrix  element*  for  a  ■ 
22.5*.  45*.  67.5’,  and  78.73*.  respectively.  The  inter** 
specular  peak  that  occurs  when  the  scattering  angle  3  ■ 
2a  causes  sharp  spikes  on  the  polarization  matrix  ele¬ 
ments.  W*  note  that  as  the  illumination  angle  is  varied, 
the  matrix  elements  varv  greatly  and  do  not  appear  to  ap¬ 
proach  a  limiting  value. 

The  experimental  matrix  elements  of  the  fiber-surface 
system  are  dissimilar  to  the  elements  of  both  tht  lone 
fiber  and  the  lone  surface.  The  oscillation  frequency  of 
the  elements  of  the  fiber-surface  system  is  much  greater 
than  that  of  the  lone  fiber  or  the  tone  surface.  This  is  not 
surprising  since  the  fiber-surface  system  approximates  a 
double-fiber  system:  fiber  and  an  image  fiber  whoa*  di¬ 
mensions  are  twice  those  of  the  kme  fiber.  In  general, 
increasing  the  linear  dimensions  of  a  scattaring  system  in¬ 
creases  the  oscillatory  frequency  of  its  scattering  matrix 
elements.  The  large  difference  among  the  matrix  ele¬ 
ments  of  the  fiber,  the  surface,  and  the  fiber-surface  sys¬ 
tem  weals  that  tha  three  systems  art  fundamentally 
different.  Therefore  one  cannot  use  tha  matrix  elements 
of  a  lone  fiber  to  predict  tha  matrix  elements  of  a  fiber- 
su’-iace  system. 


MODEL 

The  difference  in  oscillatory  frequencies  between  the  two 
systems  cannot  be  accounted  for  simply  by  treating  (ha 
fiber  as  a  perturbation  of  the  surface.  A  different  model 
must  he  used  to  predict  the  scattering.  Three  models 
v>«re  proposed  by  Nahm  and  Wolfs  to  predict  the  Matter¬ 
ing  from  spheres  on  a  conducting  plans."  Nahm  and 
Wolfs  concentrated  on  the  total  intensity  matrix  element 
Su.  Our  polarization  mea*u*»msnts  contain  additional 
information  that  doe*  not  play  a  part  in  bidirectional  re- 
f  iectanc;  distribution  function  theory.  Our  treatment  of 
this  system  uses  a  mods!  similar  to  the  double-interaction 
model.  In  this  model,  light  reaches  the  detector  by  travel¬ 


ing  one  of  four  different  paths  ( Fig.  8).  The  light  may 
either  (I)  strike  the  fiber  and  scatter  directly  to  the  detec¬ 
tor;  (II)  reflect  off  the  surface,  strike  the  fiber,  and  scat¬ 
ter  to  the  detector;  (III)  strike  the  fiber  and  scatter  to  the 
surface,  where  it  is  reflected  to  the  detector;  or  i  IV)  re¬ 
flect  off  the  surface,  strike  the  fiber,  and  scatter  to  the 
surface,  where  it  is  reflected  to  the  detector.  This  model 
includes  no  interactions  between  the  fiber  and  the  sur¬ 
face,  for  instance,  light  scattered  by  the  fiber  that  reflecta 
off  the  surface  before  interacting  with  the  fiber  again. 
The  resulting  amplitude  elements  for  these  four  rays  can 
be  derived  a*  follows: 


r.“( 9)  =■  7\(d)  +  R,[. t/2  -  atexp  i5<a)r„(d  -  2al 

+■  RJiv/2  -  (d  -  aijexp  i6(d  -  air,nJ  -  2al 
+  RJtw/2  -  (3  -  a)]ff,(n72  -  alexp  i[5(ai 
+  5<d  -  al]r„(d).  *  '3) 

Here  r„(d)  is  the  scattering  amplitude  matrix  element 
T[(3)  or  TjUH  (n  *  1  is  for  the  TE  mode,  and  n  »  2  is  for 
tha  TM  mods;  see  Ref.  8)  for  the  cylinder  with  no  surface 
present,  and 


5(d) 


4i rr  sin  d 
A 


(4) 


is  a  phase  difference  resulting  from  the  paths  that  the 
light  may  travel  before  reaching  tha  detector.  The 
Fresnel  reflectance  coefficient*  R.  are  given  by 


Mi»i  cos  d,  -  -  (ni/ni)*  sin2  d.]12 

Mani  coe  d,  +  Mi«atl  -  (rti/rta)2  sin2  d,  j1 2 


*a(d.) 


Mirta  co*  d,  -  Marti[l  -  (rt1/rtt)i  sin8  d  ]‘ 2 
Mi  rtf  coa  d,  +  Mart|[l  -  («i/rtj)*  sin*  d,  ]“  2 


1 5) 

16) 


where  the  subscript*  on  m>  and  n,  are  the  permeability  and 
the  complex  refractive  index  for  the  medium  on  the  inci¬ 
dent  side  of  the  surface  (i  *  1)  and  the  transmission  side 
of  the  surface  U  *  2). 

Figure  9  shows  the  experimental  Mueller  scattering 
matrix  and  tha  theoretical  results  for  a  0  26-Mm-radms 
quartz  fiber  resting  on  an  aluminum  surface  in  =  0 .3  + 
5.0*)1*  calculated  from  the  scattering  amplitudes  de¬ 
termined  from  Eq.  (3).  Although  the  experimental  and 


1 

i 

l 


Fig.  8.  Paths  that  an  incident  beam  may  follow  to  :  -e  ■  twr  »nd 
after  interaction  before  being  detected  on  the  inc.ar-i  -  d«  of 
the  surface. 
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Fig.  2  Th*  four  unique  Mueller  matrix  element*  calculated  for  a 
fiber-surface  system  (quart*  fiber  of  radius  0  26  nm)  illuminated 
at  ncar-grann,  incidence  (a  »  11.25*1  from  an  aluminum  surface 
i  (slid  curves)  are  compared  with  the  experimental  Mueller  me* 
}  tnx  element,  (dotted  curves). 

jP  theoretic  J  curves  do  not  coincide,  it  is  interesting  to  note 
it  number  of  similarities.  Most  important  is  that  the 
^■unbars  of  minima  and  maxima  in  the  corresponding  ele- 
mment'  of  both  sets  of  curves  are  approximately  equal, 
■fcne  of  the  crudest  ways  to  estimate  fiber  (or  sphere)  di- 
fwneter  is  to  match  the  number  of  maxima  or  minima  in  a 
^particular  interval  with  theory  (this  is  similar  to  deter- 
■Billing  slit  width  in  a  diffraction  exparimant).  Theas 
Dds  fit  the  fiber  on  the  surface  much  better  than  they  fit 
■i  lone  fiber  The  theory  for  total  intensity  matrix  eie- 
^Brat  6  ,,  fits  the  experimental  data  quite  well.  Both  show 
‘  intensity  falloff  of  approximately  3  orders  of  magni- 
from  the  specular  peak  to  -180*.  Although  this 
may  predict  the  total  intensity  of  the  scattered  light 
^Mrty  wed.  it  doee  not  accurately  predict  the  polarization 
Bitt«  of  the  scattered  light,  especially  with  regard  to  po- 
B  isation  magnitudes.  However,  the  numbers  of  maxima 
vd  minima  compart  favorably.  Comparing  polarization 
'*  is  important  from  an  experimental  point  of  view 
^Keaua*  polarization  curves  from  two  scattarara  that  give 
y  identical  Su  can  be  quite  different  It  ia  impor- 

F;  from  a  theoretical  point  of  view  because  theories  that 
predict  Su  cannot  bo  taken  far  granted  in  predicting 
correct  polarization.  That  ia  the  caee  here. 


INCLUSIONS 

te  light-scattering  ■aaaarsmants  and  analysis  of  the 
tar-surface  syeteas  lead  to  the  following  conclua tons: 
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1.  The  S,J  for  a  fiber-surface  system  in  no  way  re¬ 
sembles  the  S„  for  the  isolated  fiber.  This  means  that  the 
matrix  elements  measured  for  a  fiber  (or  sphere)  on  a  sur¬ 
face  are  not  useful  in  determining  the  properties  of  the 
defect  without  an  appropriate  theory. 

2.  S,j  are  strongly  dependent  on  the  angle  of  inci¬ 
dence  a.  Studies  of  complete  scattering  as  a  function  of 
incident  angle  may  be  necessary  to  characterize  surface 
defects. 

3.  Polarization  elements  lother  than  Su)  carry  suffi¬ 
cient  extra  information  to  warrant  their  measurement  for 
a  complete  analysis  of  the  surface. 

4.  The  most  stringent  test  of  any  theory  that  predicts 
scatter  is  whether  it  accurately  predicts  the  polarizations 
Often  surfaces  with  very  similar  Su  will  produce  very  dif¬ 
ferent  polarization  curves. 
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Two  limiting  expressions  occur  for  scattering  from  very  small  spheres.  One  occurs  when  the  refrac¬ 
tive  index  becomes  small  (Rayleigh  scattering),  and  the  other  occurs  when  the  sphere  becomes  perfectly 
conducting  (Thomson  scattering).  We  explore  the  scatter  from  small  spheres  having  real  refractive  in¬ 
dices.  For  such  spheres,  resonance  conditions  occur,  and  the  resulting  scattering  coefficients  are  no 
longer  proportional  to  the  volume  of  a  sphere. 

PACS  numberis):  42.25.Fx 


INTRODUCTION 

Electromagnetic  scattering  from  small  spheres  was  ex¬ 
plored  by  Rayleigh  [1]  and  Thomson  [2],  who  calculated 
light-scattering  expressions  for  two  very  different  special 
cases:  where  the  spheres  have  small  refractive  indices, 
and  where  they  are  perfect  conductors,  respectively. 
They  did  this  before  Mie  [3]  and  Lorenz  [4]  derived  a  for¬ 
malism  for  the  scatter  from  arbitrary  spheres.  The  range 
of  validity  of  these  special  limits  has  been  explored  by 
Kerker,  Scheiner,  and  Cooke  [3],  They  found  that  as  the 
sphere  size  becomes  smaller,  Rayleigh  theory  is  valid 
over  a  larger  range  of  refractive  index,  and  Thomson 
theory  is  valid  over  a  smaller  range  of  refractive  index. 

We  might  expect  the  scattering  that  occurs  from  small 
spheres  that  lie  in  the  region  between  the  Rayleigh  and 
the  Thomson  limits  to  be  composed  of  some  combination 
of  the  modes  present  at  these  two  limits.  However,  the 
extinction  efficiencies  of  small,  dielectric  spheres  as  a 
function  of  refractive  index  (Fig.  1)  are  not  smooth,  but 
complicated  by  sharp  resonances.  Resonances  in  the 
light  scattering  from  spheres  have  attracted  a  great  deal 
of  attention  recently  [6-16].  Resonances  appear  as 
strong,  narrow  enhancements  in  the  scattering  of  a  parti¬ 
cle.  The  large  internal  fields  that  are  built  up  within  such 
particles  can  create  interesting  effects  and  have  been  used 
to  investigate  various  phenomena  such  as  fluorescent  and 
Raman  scattering  [17-24]. 

It  is  well  known  that  as  the  sphere  size  becomes  small 


FIG.  1.  Extinction  efficiency  for  a  small  (r  «0.0U)  sphere  as 
a  function  of  real  refractive  index. 


with  respect  to  the  illuminating  wavelength,  the  equa¬ 
tions  predicting  the  scatter  are  greatly  simplified,  since 
only  a  few  sets  of  coefficients  arC  necessary  to  character¬ 
ize  the  scatter.  We  take  advantage  of  these  simplified  ex¬ 
pressions  and  derive  the  resonance  conditions  directly 
from  the  a„  and  6,  coefficients  [23-27]  rather  than  from 
the  A „,  5„,  C„,  D„  coefficients  [28],  which  is  the  stan¬ 
dard  method.  We  then  explore  the  scattering  behavior 
and  cross  sections  on  and  near  resonance.  Studying  reso¬ 
nances  in  small  spheres  gives  insight  into  the  resonances 
that  occur  in  larger  spheres  and  even  in  more  complicat¬ 
ed  particles. 

We  note  that  the  resonance  conditions  in  these  small 
spheres  are  met  when  the  sphere  refractive  index  is  large 
(m  >  ir/x ).  Therefore,  it  would  seem  that  this  work 
would  constitute  only  a  theoretical  exercise,  which  could 
only  provide  insight  into  other  resonance  situations. 
However,  in  a  recent  paper,  Scully  [29]  has  shown  via 
quantum  coherence  that  when  operating  near  an  atomic 
resonance  between  an  excited  state  and  a  coherently 
prepared  ground-state  doublet,  a  large  enhancement  of 
the  refractive  index  (by  many  orders  of  magnitude)  may 
be  achieved  with  zero  absorption.  In  this  case,  the  light- 
scattering  resonances  that  we  examine  cannot  only  be 
realized,  but  may  prove  to  be  a  useful  tool  in  characteriz¬ 
ing  the  optical  propenies  of  such  materials. 


SCATTERING  COEFFICIENTS 


The  electromagnetic  scattering  a  large  distance  from  a 
sphere  of  radius  r,  illuminated  by  a  unit-normalized  plane 
wave  traveling  in  the  positive  z  direction,  and  polarized 
in  the  2  direction,  may  be  expressed  by  two  scattering 
amplitude  functions  given  by 
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where  5|  is  measured  in  they-z  plane  and  S:  is  measured 
in  the  x  - z  plane.  The  Mie  scattering  coefficients  are  given 
by 

m  0',  ( x  )0,  ( mx )  -  0',  ( mx  )0,  ( x ) 


a.= 


b.= 


m|f,(x)0„(mx)-0;,(mx)f,(x) 
m0',(mx)0„(xl  — 0',(x)0,(mx) 
m  0',  ( m x  )|,  ( x )  -  &  ( x  )0,  ( mx  ) 


(2a) 


where  m  is  the  complex  refractive  index  of  the  sphere, 
x-litr/K  and  i b„  and  5,  are  the  Riccati-Bessel  func¬ 
tions. 

First  we  examine  what  happens  for  the  special  limiting 


case  where  the  sphere  size  is  small  ir  «X).  In  this  limit¬ 
ing  case  the  scatter  is  determined  primarily  by  the 
lowest-order  terms  of  the  series  given  by  Eq.  (1).  The 
Riccati-Bessel  functions  for  the  n  -  I  case  are  given  by 


0,  (p> 


sinp  _ 


ccsp  ,  £,(p)=exp(ip)(  —  ip~'  —  1 )  .  (3) 


(2b)  For  small  arguments,  these  functions  are  approximately 


(4) 


For  small  x,  the  scattering  coefficients  given  by  Eq.  (2) 
are  approximately 
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Similarly,  the  second-order  coefficients  may  be  written  as  in  Eq.  (6). 


+x ' 
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2m 2 
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Two  additional  limiting  conditions  exist  for  small 
spheres.  One  is  the  Rayleigh  limit  valid  when  I  mix  « 1. 
In  this  case  the  scattering  coefficients  further  reduce  to 


2ix 1  m z  —  1 
3  m 2  +  2  ’ 

a,  -bn  -0  for  n  >  1  . 


6,-0  , 


(7) 


The  other  limiting  condition  occurs  when  the  sphere’s  re¬ 
fractive  index  approaches  that  of  a  perfect  conductor 
m  —i<x.  In  this  case  the  scattering  coefficients  reduce  to 


o, 


a,  —6,  -0  for  n  >  1  . 


(8) 


Note  that  in  both  these  limiting  cases  the  scattering 
coefficients  are  proportional  to  x  .  Figure  2  shows  the 
angular  scattering-intensity  distributions  for  a  Rayleigh 
and  a  Thomson  sphere. 


RESONANCES 

We  now  examine  the  6,  mode  in  more  detail.  When 
the  refractive  index  m  is  increased  along  the  real  axis,  the 
sine  terms  in  Eq.  (5)  do  not  contribute  to  the  scatter  when 
mx  =*Nir  where  the  index  N  is  an  integer.  Resonances  in 
the  6,  coefficients  occur  at  approximately  these  locations. 
Two  interesting  results  occur  that  are  worth  pointing  out. 
The  first  is  that  the  scattered  fields  are  no  longer  propor¬ 


tional  to  the  sphere  volume.  As  a  result,  the  scattering 
efficiencies  of  such  particles  are  greatly  increased  corn- 
paired  to  particles  just  off  resonance.  Figure  3  shows  the 
first  by  resonance  (following  a  previous  convention  [8]; 
this  is  written  as  6  J).  For  purposes  of  illustration,  we 
chose  to  examine  small  spheres  with  an  arbitrary  but 
definite  radius  r=0.01X.  The  resulting  resonances  will 
necessarily  occur  only  at  large  values  of  refractive  index 
(m  *  50 N)  due  to  the  small  size  of  these  spheres. 

The  second  interesting  result  is  that  the  mode  of  oscil- 


FIG.  2.  Angular  scattering-intensity  distributions  for  an 
r*0.0lX,  m  ”2.0  Rayleigh  sphere  (O),  and  an  r“0.0lX  Thom¬ 
son  sphere  ( X ). 
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FIG.  3.  Extinction  efficiency  for  a  small  (r=0.0U)  sphere  as 
a  function  of  real  refractive  index  near  the  b  \  resonance.  Also 
shown  are  the  scattering  intensities  for  the  TE  and  TM  modes 
on  the  6 5  resonance  (m  =» 49.98006). 


lation  for  the  resonant  sphere  (bt)  is  completely  different 
than  that  for  the  Rayleigh  sphere  (<fp.  The  resulting  field 
distributions  (also  shown  in  Fig.  3)  for  the  h, -resonant 
sphere  will  necessarily  be  different  from  the  Rayleigh 
sphere  (shown  in  Fig.  2).  The  TE  and  TM  Rayleigh- 
sphere  intensities  are  proportional  to  the  TM  and  TE  6,- 
resonant  sphere  intensities,  respectively.  For  the  6, 
mode,  the  incident  electromagnetic  field  induces  a  dipole 
moment  perpendicular  to  the  incident  electric  field. 

The  resonances  do  not  account  for  all  the  reso¬ 
nances  shown  in  Fig.  1.  Resonances  in  the  a ,  and  b: 
modes  occur  when  tan(mx)  is  approximately  equal  to  mx 
[for  large  mx,  this  occurs  approximately  when 
mx  —iN  +  {hr].  The  a  J  and  b\  resonances  are  shown  in 
Fig.  4.  These  resonances  are  much  narrower  [half-width 
of  the  order  &(mjc)/(mx)~  10-7]  than  the  b  j  resonances 
[half-width  of  the  order  A(mx)/(mx)-~  10“5].  Reso¬ 
nances  in  the  higher-order  modes  are  not  as  prominent, 
since  the  resonant  fields  are  proportional  to  xk  where 
k>  3. 

Equations  (3)  and  (6)  not  only  can  be  used  to  predict 
where  resonances  occur,  but  can  also  provide  information 
on  the  shapes  of  the  resonances.  We  will  now  take  a 
closer  look  at  the  b,  resonance.  If  we  express  the  com¬ 
plex  refractive  index  as  m  mmr+im,,  where  m,  and  m, 
are  both  real  quantities,  we  can  expand  the  sine  and 
cosine  function  about  the  resonance  locations. 


7 1  *S  Retractive  index  71  53 


FIG.  4.  Extinction  efficiency  for  a  small  (r  «0.0U>  sphere  as 
a  function  of  real  refractive  index  near  the  a  \  and  b  l  resonances. 
Also  shown  are  the  scattering  intensities  for  the  TE  and  TM 
modes  on  the  oj  (m  *71.50079)  and  b\  resonances 
(rvs  —71.510 13). 
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light-scattering  resonances  in  small  spheres 
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sini  mx )  =  (  -  1  >'vsin( mx  -Sir) 

-(  - 1  )s[i  sinh(m,x  H-Acoshlm,* )]  , 
cos(  mjt ) * ( —  1  )‘vcos<  mx  - Nv ) 

~( -  1  )-v[ coshlm, x)-/Asinh(m,;r )]  , 


(9) 


where  m,x-Sir,  and  A  =  m, x  -  Sir  , 
Near  resonance,  the  scattering  coefficient  b{  can  be 
simplified: 

1  2  « 

x  m _ 1_ 
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x  cosh(m,x )  +  /' 


sinh(m,x ) 


— sinh( m,x)  —  i  cosh<  m,x ) 
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where  A  = - 


—  im  coshlm.* ) 


(10) 


FIG.  5.  Extinction  efficiency  near  the  6!  resonances:  (a)  for 
an  m “30.0  sphere  where  N“  l;  <b>  for  an  m - 100.0  sphere 
where  N“2;  (cl  for  an  m  “50.0  sphere  where  N  “2:  (d)  for  an 
m  “5O.0+0.0W  sphere  where  (V“l;  (e)  for  an  m  “50.0+O.OW 
sphere  where  .V“l;  and  (0  for  an  m  “50.0+0.  If  sphere  .V“  I. 
Radius  r“0.0lA,. 


Equation  ( 10)  may  be  simplified  if  we  assume  the  absorp¬ 
tion  is  small  lm,jt  «  1 ), 


x  +mm,  t  (A+x  /m  ) 
ix1/m  +  im,x 

Multiplying  Eq.  (11)  by  its  complex  conjugate  yields  a 
Lorentzian  function  centered  at  A —  —x/m,  with  a  half¬ 
width  of  x1/m  +m,x.  We  note  that  as  the  absorption  is 
increased  from  zero,  the  amplitude  of  b,  will  decrease, 
and  the  half-width  will  increase.  When  mm,  »x,  the 
amplitude  will  be  proportional  to  1/m,  and  the  half¬ 
width  will  be  proportional  to  m,.  These  dependencies 
were  reached  empirically  for  larger  spheres  having 
moderate  refractive  indices  [9].  This  type  of  analysis  may 
also  be  performed  for  the  a,  and  b2  resonances. 

The  dependence  of  the  line  shapes  on  the  size  and  re¬ 
fractive  index  is  shown  in  Fig.  S.  In  Fif,  3(a),  the 
m  —  50.0  sphere  passes  through  the  b  |  resonance  as  its 
radius  is  increased.  In  Figs.  5(b)  and  5(c),  we  can  exam¬ 
ine  the  b  f  resonances  as  the  radius  of  an  m  — 100.0 
sphere  and  an  m  =  50.0  sphere  is  increased,  respectively. 
Going  to  a  higher  index  by  increasing  the  refractive  index 
[Fig.  5(b)]  results  in  a  narrower  resonance.  Going  to  a 
higher  index  by  increasing  the  size  parameter  [Fig.  5(c)] 
results  in  a  broader  resonance.  The  latter  result  has  been 
discussed  previously  for  spheres  much  larger  than  the 
wavelength  [7,10,11],  Figure  5  also  shows  the  shape  and 
the  width  of  the  b\  resonance  as  the  refractive  index  is 
changed  from  m  «  50.0+0.  OOIt  [Fig.  5(b)]  to 


Anpald^) 


FIG.  6.  Angular  scattering-intensity  distributions  for  a  small 
<r“ 0.0U)  sphere  having  refractive  indices:  (a)  m  “43.666  10, 
(b)  m  -  59. 226  10,  (c)  m  “  71. 50*  44,  and  (d)  m  “  71 . 5 19  92. 
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m  =50.0+0.01/  (Fig.  5(e)]  to  m  =50.0+0. 1/  [Fig.  5(f)]. 
These  figures  verify  that  increasing  the  imaginary  part  of 
the  refractive  index  results  in  a  reduction  of  the  height 
and  an  increase  in  the  width  of  the  resonant  peak.  This 
conclusion  has  also  been  discussed  when  examining  large 
spheres  [7-9,13]. 

The  angular  intensity  distributions  for  spheres  having 
refractive  indices  between  resonant  values  do  not  neces¬ 
sarily  resemble  the  distributions  of  either  a  Rayleigh  or  a 
Thomson  sphere.  Figure  6  shows  that  cases  exist  when 
the  forward  scatter  id?-0*)  is  down  several  orders  of  mag¬ 
nitude  [in  cases  of  Figs.  6(b)  and  6(c),  a ,  +6 ,  +  5b 2  /3  -  0] 
and  when  the  backscatter  ti?~  180")  is  down  several  or¬ 
ders  of  magnitude  [in  cases  of  Flgs.,6(a)  and  6(d), 
a,  —  h,  +  562/3~0].  These  are  important  points  because 
most  scattering  studies  of  spheres  would’jbad  one  to  be¬ 
lieve  that  the  intensity  in  the  forward-scafering  or  back- 
scattering  directions  would  net  emtnd  liberal  orders  of 
magnitude  below  the  scattered  intensity  at  other  scatter¬ 
ing  angles. 

- 


SUMMARY 

For  the  limiting  case  when  sphere  size  is  much  smaller 
than  the  incident  wavelength,  the  equations  describing 
the  scatter  are  greatly  simplified.  When  two  additional 
limiting  conditions  on  the  refractive  index  are  applied, 
these  equations  are  simplified  further.  When  the  refrac¬ 
tive  index  is  increased  along  the  real  axis,  resonance  con¬ 
ditions  develop  that  complicate  the  resulting  scatter. 
Precisely  because  these  resonances  occur,  no  limiting 
condition  can  be  reached  as  the  complex  refractive  index 
is  increased  along  the  real  axis.  Resonances  have  been 
studied  extensively  for  large  spheres.  Studying  the  reso¬ 
nances  in  smaller  spheres  in  which  the  equations  are 
greatly  simplified  gives  insight  to  the  resonance  beha.ior 
occurring  in  larger  spheres. 
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A  facet  model  is  used  to  predict  the  polarization  states  of  light  scattered  from  rough  surfaces.  The  states  are 
compared  with  experimentally  determined  elements  from  various  types  of  rough  substrates.  The  experimental 
data  are  found  to  be  quite  similar  to  what  is  predicted  by  the  model. 


INTRODUCTION 

Before  the  development  of  Maxwell’s  equations  and  an 
electromagnetic  theory  of  light,  Fresnel  was  able  to  de¬ 
duce  the  reflection  and  transmission  coefficients  for  light 
interacting  with  a  perfectly  smooth  surface.  Since  then 
researchers  have  been  developing  theories  to  obtain  sur¬ 
face  characteristics  (such  as  the  rms  roughness)  from  the 
light  scattered  from  the  surface.  These  models  often  de¬ 
velop  along  a  Fourier- transform  approach.  The  bidirec¬ 
tional  reflectance  distribution  function  is  a  standard 
measurement  of  the  scattered  radiance  divided  by  the  in¬ 
cident  irradiance  on  the  surface.14  It  usually  does  not 
take  the  polarization  states  of  the  source  or  of  the  scat¬ 
tered  light  into  consideration.  Therefore  valuable  infor¬ 
mation  is  missing.  To  characterize  the  light  scattered 
from  a  system  completely,  one  can  measure  the  Mueller 
scattering  matrix.*"1  This  4X4  matrix  contains  all  the 
polarization  scattering  information. 

Recently  there  has  been  considerable  interest  in  light 
scattered  from  rough  surfaces.*"11  These  studies  include 
measurements  of  the  scattered  a-  and  p-polarization  states. 
In  only  a  few  instances  have  authors  characterized  the  po¬ 
larization  states  more  completely  by  expressing  the  result¬ 
ing  scattered  light  in  the  form  of  a  Mueller  matrix.14"* 
The  purpose  of  this  paper  is  to  take  e  closer  look  at  the 
polarization  state  of  the  light  It  is  not  our  interest  to  ex¬ 
amine  the  intensities,  which  have  been  examined  in  detail 
in  the  papers  cited  above.  The  simplest  model  that  can 
predict  the  polarization  state  of  light  scattered  from  a 
rough  surface  uses  a  geometric-optics  approach  in  which 
the  light  is  Fresnel  reflected  from  planar  surface  facets.17 
Although  the  validity  of  this  model  is  questionable  in 
nearly  all  physical  situations  (it  is  valid  only  when  the 
wavelength  approaches  zero),  the  model  does  have  one  dis¬ 
tinct  advantage  in  that  the  calculations  are  easily  made. 
In  this  paper  we  dbcuss  and  develop  this  facet  model  We 
then  compare  the  Mueller  scattering  matrices  predicted 
by  this  facet  model  with  experimental  matrices  for  several 
surfaces  having  various  surface  roughnesses. 

FACET  MODEL 

Figure  1  shows  the  geometry  of  the  scattering  system. 
The  scattering  surface  is  located  perpendicular  to  the 
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z  axis  in  the  x-y  plane.  The  incident  radiation  is  a  plane 
wave  traveling  in  the  x-z  plane,  oriented  at  angle  d,,*  with 
respect  to  the  z  axis.  The  wavelength,  wfove  vector,  and 
refractive  index  for  the  plane  wave  in  the  nonabsorbing, 
nonmagnetic  incident  medium  above  the  surface  are  A,  k,, 
and  nt,  respectively.  The  complex  wave  vector  and  re¬ 
fractive  index  for  a  plane  wove  of  the  same  frequency  in 
the  medium  below  the  surface  are  k>  and  n»,  respectively. 
A  detector  having  area  A* ,  at  a  point  de¬ 

tects  light  scattered  into  the  solid  angle 

faai  An,)* 

The  facet  model  used  in  this  study  is  based  on  a 
geometric-optics  approach  in  which  ray  optica  are  valid. 
Light  is  assumed  to  travel  in  straight  lines.  When  a  ray 
of  light  intersects  a  surface,  it  is  Fresnel  reflected  at  that 
location  in  such  a  way  that  the  incident  angle  at  that  loca¬ 
tion  is  equal  to  the  reflection  angle  where  these  angles  are 
measured  from  the  normal  to  a  tangent  plane  (or  surface 
facet)  at  that  location  (Fig.  2).  Ws  further  simplify  our 
model  by  considering  only  the  light  scattered  in  the  x-z 
plane  (pM  -  0)  and  ignoring  rays  that  reflect  off  of  more 
than  one  surface  facet.  These  limitations  eliminate  the 
ability  of  the  model  to  predict  any  cross  polarization 
(s  -*  p  or  p  —  s)  or  effects  such  as  enhanced  backscatter. 
Of  course  it  is  possible  to  develop  a  physical  theory  that 
has  more  basis  in  reality,  but  it  is  always  interesting  to  see 
how  well  a  zeroth -order  solution  predicts  experimental  re¬ 
sults.  And  perhaps  for  some  applications  the  simplest  so¬ 
lution  is  all  that  is  required. 

From  Fig.  2  we  see  that  only  the  facets  oriented  at  angle 
d  with  the  incident  light  rays  will  reflect  rays  and  contrib¬ 
ute  to  the  amplitude  in  the  direction  d*.'.  Upon  reflec¬ 
tion,  the  light  is  attenuated  by  a  Fresnel  reflection  factor 
that  depends  on  the  state  of  linear  polarization  of  the  inci¬ 
dent  light  and  complex  refractive  indices  of  the  media 
above  and  below  the  surface: 


9  ...  *1  cos  d  -  n»[l  -  (ni/n,)1  sin*  d]14 

cos  d  +  n,(l  -  (n./n,)1  sin1  d]"’ 

R  (A\  n*  609  *  ~  ~  ,in>  d]*” 

1  "  n,  cos  d  +  n,[l  -  («,/«,)» sin1  d]1'1’ 

where 
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Fig.  2.  Geometry  of  the  facet  reflection  of  a  rough  surface. 


d  -  (#*.'  -  »)/2.  (3) 

The  aubacripta  TE  end  TM  correapond  to  the  incident 
electric*  end  magnetic-field  vectora  being  tranaverae  to 
the  scattering  plane  (parallel  to  the  y  axis),  respectively. 
We  can  calculate  the  electric-far-field  amplitude  by  sum¬ 
ming  the  contributions  of  the  facets  reflecting  light  in  the 
direction  [a  time  dependence  of  expftatf)  is  assumed]: 

_  _  Et.'-Rt.W  (  k,  •  A(r) 

Et. - -—J-— 

X  txp(i(kt  —  k|)  ■  rja.^-i^di,  (4) 

where  the  T*  is  need  to  represent  the  TE  or  the  TM  mode 
of  the  scattered  and  the  incident  electric  fields  (£?•*“  and 
Er.im,  respectively),  is  the  cross-sectional  area  of  the 
incident  beam  upon  the  substrate,  and  dir)  is  the  unit  nor¬ 
mal  at  position  r  on  the  surface.  The  dot  product  is  a 
geometric  size  factor,  the  exponential  function  accounts 
for  the  phaae  differences  from  rays  reflecting  at  different 
locations  on  the  surface,  and  the  Kronecksr  delta  function 
ensures  that  only  rays  obeying  Snail’s  law  contribute  to 


the  amplitude.  In  terms  of  the  scattering  amplitude  ma¬ 
trix,  this  may  be  expressed  as 


(note  that  Sj  -  S«  -  0  for  light  scattered  in  the  x-z  plane 
for  the  facet  model,  since  there  is  no  cross  polarization). 
Of  the  16  elements  of  the  Mueller  matrix,  only  4  are 
nonzero  or  simple  multiples  of  other  elements  in  this  case. 
We  use  a  normalized  Mueller  scattering  matrix  given  by 


Su'  ”  S„  V4(|5i|*  +  |Si|’),  (6) 

Sis'  “  Sa/sa  -  Vi  (|S,|*  -  ISaIVS,, ,  (7) 

Sja  ™  Sja/Sn  “  Re(SiS,*)/Su ,  (8) 

Su  ”  Su/Sn  “  Im(SiSj*)/Su .  •"  (9) 


When  we  normalize  the  polarization  matrix  elements  by 
the  total  intensity  matrix  element  Su  (expressing  them  as 
a  percent  polarization),  the  integral  in  Eq.  (4)  cancels 
from  these  elements  and  contributes  only  to  the  total  in¬ 
tensity  matrix  element  (Su).  The  polarization  matrix  ele¬ 
ments  (which  we  want  to  examine)  are  not  affected  by 
this  function.  In  our  calculations  of  matrix  element  Su, 
we  arbitrarily  set  this  integral  to  unity.  For  a  perfectly 
conducting  substrate  (it,  -*  <*>),  the  scatter  would  then  be 
independent  of  the  scattering  angle.  Any  deviation  in  the 
scatter  from  a  horizontal  line  (when  plotted  as  a  function 
of  the  scattering  angle)  is  due  to  the  complex  refractive 
index  of  the  substrate. 

Finally,  by  expressing  the  matrix  elements  as  a  function 
of  the  scattering  angle  dn',  which  is  measured  from  the 
path  of  the  incident  radiation,  we  find  that  the  polarization 
state  [Eqs.  (7)-(9)]  is  independent  of  the  incident  angle 
dta,.  The  light- scattering  Mueller  matrices  predicted  by 
the  facet  model  for  a  copper  substrate  (n,  -  1.1  -  2.5i  at 
A  -  0.4416  (un)a  and  an  aluminum  substrate  (n,  -  0.5  - 
5-0i  at  A  —  0.4416  #tm),“  both  illuminated  at  d,„  —  90*. 
are  shown  as  a  function  of  scattering  angle  in  Fig.  3. 
We  choee  this  incident  angle  (grazing  incidence)  because  it 
provides  the  maximum  amount  of  information  in  the 
graphs.  Choosing  a  nongrazing  incidant  angle  would  re¬ 
move  only  the  leftmost  portions  of  the  graphs  (the  forward 
scatter)  owing  to  vignetting  by  the  substrate.  The  differ¬ 
ences  between  these  two  seta  of  matrix  elements  is  due 
mainly  to  the  differences  in  the  abeorption  coefficients. 
These  sets  of  matrix  elements  ere  compared  with  the 
experimentally  measured  matrix  elements  for  the  rough 
surfaces. 


LIGHT-SCATTERING  MEASUREMENT 

The  polar  nephelometer  used  in  this  study  employs  the  po¬ 
larization  modulation  technique  developed  by  Hunt  and 
Huffman.**  A  complete  discussion  of  the  nephelometer 
design  and  operation  is  given  in  Ref.  20  and  also  by  Perry 
et  a l.,n  Bickel  et  al.,a  Bell,"  and  Iaftlks.*4  The  measure¬ 
ment  requires  modulating  the  incident  beam’s  polarization 
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Fig.  5.  Experimental  light-ecattering  Mueller  matrix  alemanta  of  a  rough  aluminum  eubatrata  illuminated  at  6i»  -  0*.  16*.  30*.  45*.  60*. 
and  75*  U  -0.4416  pm). 


state  at  u),  -  50  kHz  with  a  photoelaatic  modulator  while 
observing  the  signals  carried  by  the  fundamental  modula¬ 
tion  frequency  (a»«)  and  second  harmonic  {2m,)  of  the  scat¬ 
tered  light  (A  -  0.4416  pm).  A  lock-in  amplifier  is  used 
to  demodulate  the  signals.  The  experimental  scattering 
matrix  S„*(0)  can  be  measured  with  the  proper  choice  of 
entrance  and  exit  optica.  Normalization  of  the  matrix  ele¬ 
ments  is  performed  fay  servoing  the  photomultiplier  tube 
(RCA  1P21)  gain  via  a  constant-current  servo,  requiring 
a  constant-dc  output  signal  from  the  detector.  This  is 
performed  over  the  entire  scan  of  the  detector.  The 
normalized-output  S/(d)  analog  signals  are  collected  and 
sent  to  a  computer.  The  experimental  matrix  elements 
measured  are  actually  combinations  of  the  Mueller  matrix: 


Su*-«a,  (10) 

Sis*  ■  Sjs/Sm,  (11) 

S„*  -  <3»  +  S»)/(S„  +  S„) ,  (12) 

S\t*  “  (Si«  +  Sm)/(5u  +  S„).  (13) 


Since  there  is  no  mixing  of  the  s-  and  the  p-polarization 
states  in  the  facet  model  (Sj  —  S«  —  0),  the  matrix  ele¬ 
ment  S«  -  Su  “  Sl4  -  0.  The  normalized  matrix  ele¬ 
ments  predicted  by  the  model  are  equivalent  to  the  - 
measured  combinations  of  matrix  elements  (S</  -  S„*). 


SCATTERING  FROM  RANDOMLY 
SCRATCHED  SUBSTRATES 

The  first  sets  of  experimental  data  are  from  randomly 
scratched  metal  substrates.  Figure  4  shows  the  four  in¬ 
dependent  light-scattering  Mueller  matrix  elements  for  a 
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Fig.  6.  Micrograph  of  randomly  sanded  aluminum  substrate. 
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Micrographs  of  thraa  diffsrant  substrata*  roughen *d  by 


randomly  scratched  copper  substrate  as  a  function  of  scat¬ 
tering  angle  ««,’  at  six  different  incident  angles  -  0*, 
15*.  30*.  45’,  60*.  and  75*).  Note  that  when  IdJ  <  90* 
(when  there  is  no  vignetting),  the  polarization  elements 
shown  in  Fig.  4  are  very  similar  to  the  polarization  ele¬ 
ments  shown  in  Fig.  3  predicted  by  the  facet  model. 
Figure  5  shows  the  four  independent  light-scattering 
Mueller  matrix  elements  lor  a  randomly  scratched  alu¬ 
minum  substrate  as  a  function  of  scattering  angle 
illuminated  at  the  sauna  incident  angles  as  those  of  Fig.  4. 
Note  that  the  peaks  (fat  element  Su),  valleys  (in  element 
SM),  and  zero  crossing  points  (in  element  S„)  Tor  the  alu¬ 
minum  substrata  are  shifted  slightly  toward  smaller 
angles,  as  can  be  seen  in  the  theoretical  elements  pre¬ 
dicted  in  F%.  3.  Moat  notable  in  the  experimental  set  of 
matrix  elements  is  that  element  S*  shows  more  depen¬ 
dence  on  incident  angle  for  the  aluminum  substrate  than 
for  the  copper  substrata.  Note  also  that,  as  the  incident 
angle  ie  increased  (toward  grazing  incidance),  the  percent 
polarization  of  element  S»  rises  toward  zero  in  the 
backscatter  (d«,  -  180*).  A  micrograph  of  tha  aluminum 
substrate  used  in  this  study  is  shown  in  Fig.  8. 


We  note  that  the  experimental  matrix  elements  of  the 
aluminum  substrate  resemble  the  theoretical  matrix  ele¬ 
ments  of  the  copper  substrat  3  of  Fig.  3  better  than  they  do 
the  theoretical  elements  of  the  aluminum  substrate.  One 
possible  reason  for  this  is  that  our  experimental  substrate 
is  not  pure  aluminum.  When  exposed  to  air,  the  sample 
is  coated  by  a  layer  of  aluminum  oxide  (this  occurs  much 
more  rapidly  for  aluminum  than  for  copper).  One  effect 
of  this  is  a  decrease  in  the  absorption  (the  imaginary  part 
of  the  complex  refractive  index).  The  effective  complex 
refractive  index  of  the  oxidized  aluminum  substrate  is 
closer  to  that  of  copper  than  pure  aluminum. 


SCATTERING  AS  A  FUNCTION  OF  SURFACE 
ROUGHNESS 

The  last  set  of  matrix  elements  ie  from  aluminum  sub¬ 
strates  that  were  roughened  to  different  amounts  by  sand¬ 
blasting.  Eleven  nominally  smooth  mirror  substrates 
were  polished  to  high  reflectivity.  Substrate  0  remained 
a  mirror  and  served  as  e  reference,  while  the  other  10  mir¬ 
rors  were  sandblasted  to  increasing  roughness  (substrate  1 
was  blasted  for  1  s,  substrate  2  was  blasted  for  2  s,  etc.). 
This  created  a  family  of  related  surfaces  in  which  the 
roughness  ranged  from  zero  (a  mirror)  to  saturation  (when 
additional  sandblasting  <bd  not  change  the  surface).  Mi¬ 
crographs  of  three  of  these  surfaces  are  shown  in  Fig.  7. 

Figure  8  shows  the  four  independent  light-scattering 
Mueller  matrix  elements  for  the  reference  (mirror)  sub¬ 
strate  as  a  function  of  scattering  angle  ’  at  five  differ¬ 
ent  incident  angles  (dta  “  IS*.  30*,  4S*.  60*,  and  75*).  Su 
is  characterized  by  a  strong  specular  peak  at  the  angle  of 
reflection  and  a  steep  falloff  in  intensity  away  from  the 
specular  peak.  The  polarization  matrix  elements  of  this 
substrate  only  slightly  resemble  those  predicted  by  the 
facet  model  shown  in  Fig.  3.  Note  that  there  is  e  strong 
dependence  on  incident  angle  in  all  three  seta  of  polariza¬ 
tion  matrix  element a  in  Fig.  8,  not  just  5a  as  was  the  case 
for  the  scratched  substrate. 

Figure  9  shows  the  four  independent  light-scattering 
Mueller  matrix  elements  for  substrate  4  (sandblasted  for 
4  s)  as  a  function  of  scattering  angle  at  five  different 

incident  angles  (A*  *  15*.  30*.  4S*.  60*.  and  75*).  This 
set  of  matrix  elements  resembles  tboee  of  the  facet  model 
more  closely  (especially  matrix  element  S|«).  In  the 
backscatter  region,  elements  Sa  and  Sa  still  show  a  strong 
dependence  on  incident  angle,  becoming  more  positive  as 
is  increased.  Element  Su,  however,  is  almost  con¬ 
stant  in  tha  backscatter  region  as  di««  is  increased. 

Figure  10  shows  the  four  independent  light-scattering 
Mueller  matrix  elements  for  substrate  10  (sandblasted  for 
10  a,  almost  to  saturation)  aa  a  function  of  scattering 
angle  d*.'  at  five  different  incident  angles  (di»  —  15*.  30*. 
45*,  60*,  and  75*).  Tbs  scattering  elements  of  this  sub¬ 
strate  ere  very  similar  to  those  predicted  by  the  facet 
model.  Again,  however,  a  definite  dependence  on  incident 
angle  in  the  backscatter  region  of  Sa.  which  is  not  pre¬ 
dicted  by  the  facet  model,  occurs.  As  seen  in  the  matrix 
elements  for  the  other  aluminum  samples,  element  Sa  be¬ 
comes  more  positive  as  d*.  is  increased. 


Scattering  Angle 


Scattering  Angle 


Fig.  &  Experimental  light-featuring  Mueller  matrix  elemenU  of  lubatrau  4  (•lightly  roughened)  illuminated  at  di*  —  15*.  30*,  45’.  60 
and  76*  (A  -  0.4418  urn). 
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Fig.  10.  Experimental  light-scattering  Mueller  matrix  dementi  of  roughened  subatratc  10  illuminated  at  9im  —  15*.  30*.  45*.  60*.  and  75* 
(a  -  0.4416  Mno- 


CONCLUSION 

In  this  study  we  used  a  facet  mode l  to  predict  the  polariza¬ 
tion  states  of  light  scattered  from  a  rough  surface.  This 
model  predicts  that  the  polarization  matrix  elements  de¬ 
pend  solely  on  the  complex  refractive  index  of  the  surface 
and  are  entirely  independent  of  the  angle  of  surface  illu¬ 
mination.  We  emphasize  that  this  is  a  simple  model  that 
can  explain  only  gross  scattering  features.  It  is  inade¬ 
quate  in  explaining  the  c rose-polarization  effects  (*  — »  p 
and  p  —  t)  that  occur  in  light  scattered  from  rough  sur¬ 
faces.  Nevertheless,  we  find  that  the  experimental  ma¬ 
trix  elements  from  substrates  of  two  different  materials 
roughened  in  different  ways  sgree  quits  well  with  this 
simple  model. 
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Abstract 

We  examine  four  different  theories  which  predict  the  scattered 
radiation  from  a  system  composed  of  a  small  sphere  and  a  plane  and  show 
that  in  the  far  field  the  scattering  amplitude  components  predicted  by  these 
theories  are  essentially  the  same.  We  express  the  scattering  intensities  in 
Mueller  matrix  representation  and  examine  the  far-field  Mueller  matrix  as  a 
function  of  various  parameters. 

•  Gorden  Videen  is  now  at  Dalhousie  University,  Physics  Department,  Halifax 
N.  S.  Canada  B3H  3J5 
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Introduction 

Electromagnetic  scattering  from  isolated  systems  has  become  well- 
established  in  the  past  century.1'1  Lord  Rayleigh  did  much  to  initiate  this  re¬ 
search  when  he  solved  for  the  scattered  electromagnetic  field  from  a  small 
sphere  or  dipole  to  explain  the  color  of  skylight.4  The  scattered  electromag¬ 
netic  fields  surrounding  an  isolated  cylinder  were  solved  by  Rayleigh*^  in  1881 
and  by  von  Ignatowsky*  in  1905.  The  scattering  from  a  sphere  has  been  derived 
numerous  times  and  its  complex  history  has  been  analyzed  extensively.*'™ 

More  recently,  researchers  have  concentrated  on  scattering  from 
more  complex  systems.  Scattering  by  multiple  cylinders  by  Twersky,*  using  an 
order-of- scattering  approach;  Row, 10  using  a  Green’s  function  approach;  and 

Olaofe;"  and  Yousif  and  Kohler’*'’1  by  satisfying  the  boundary  conditions  at  the 

1 

cylinder  surfaces.  The  two-sphere  system  has  been  examined  by  Liang  and  Lo“ 
and  Bruning  and  Lo"  by  satisfying  boundary  conditions  at  the  surfaces  of  both 
spheres.  The  order-of-scattering  approach  has  been  examined  by  Fuller  and 
Kattawar. Scattering  systems  composed  of  a  cylinder,  sphere  or  dipole  rest¬ 
ing  on  or  near  a  plane  surface  have  been  explored  by  Rao  and  Barakat,"  Bob- 
bert  and  Vlleger,1*  Wind,  Vlieger  and  Bedeaux,**  Nahm  and  Wolfe,*1  Vldeen,** 
Videen.  Wolfe  and  Bickel,’*  Undell  et  mJ,u  and  Muinonen  et  a/.** 

This  last  class  of  problems  is  closely  related  to  the  Sommerfeid***r 
half-space  problem  in  which  a  radiating  dipole  is  located  near  an  interface.  In 
the  Sommerfeld  problem,  a  dipole  is  oscillating  vertically  or  horizontally  to 
the  interface,  and  is  the  source  of  the  electromagnetic  waves.  In  the  scattering 


3 


problems,  the  source  is  typically  at  infinity  producing  plane  waves  at  the 
interface.  These  interact  with  the  surface  interface  and  the  object  placed  near 
it.  This  scattering  system  has  been  discussed  extensively  by  Banos,21  Wait,2* 
and  Rahmat-Samii,  Mittra,  and  Parhami.10 

In  this  manuscript  we  investigate  four  different  theories  that  predict 
the  scattered  far-field  radiation  from  a  small  sphere  placed  near  a  surface  and 
show  that  they  all  yield  essentially  the  same  results.  Error  analysis  of  the  re¬ 
sults  are  considered  using  an  order-of-scattering  approach.  We  then  examine 
these  results  in  the  form  of  the  Mueller  matrix.11 

The  Solutions 

We  will  be  discussing  four  different  solutions  from  different  refe¬ 
rence  sources  that  have  not  used  the  same  mathematical  and  physical  conven¬ 
tions.  We  therefore  adopted  the  following  description  and  used  it  for  this 
paper.  Differences  that  arise  in  other  manuscripts  will  be  stated.  Figure  1 
shows  the  geometry  of  the  system.  A  small  sphere  of  refractive  index  m.*  and 
radius  a  is  suspended  a  distance  d  above  a  surface  separating  two  media:  free 
space  where  z  >  0  and  a  medium  of  complex  refractive  index  nw  where  z  <  0. 
The  incident  plane  wave  of  amplitude  E’*  travels  in  the  x-z  plane  and  is  elec¬ 
trically  polarized  either  parallel  to  the  y  axis  (TE)  or  parallel  to  the  vector  e 
(TM>  given  by 

l  >  S  coso  -  2  si  not  (1) 


4 


where  a  is  the  angle  of  incidence  of  the  plane  wave  measured  from  the  z  axis. 
To  make  the  dipole  approximation,  it  is  necessary  that  lm.**la  «  1.  The  time 
dependence  used  in  this  paper  is  expl-iut).  The  results  will  be  in  the  form  of 
the  scattering  amplitude  matrix  given  by 

(E5‘  )  (s*  Sj\  ( E't"u\ 

\ Si  Si j  [e z)  j2) 

where  E£*  and  E“*  are  the  scattered  electric  field  components  in  the  d  and  <p 
directions,  respectively,  and  k  is  the  wave  vector  for  the  medium  above  the 
surface. 

L  Sommerfeld- Rayleigh  System 

The  first  theoretical  treatment  we  examine  is  a  modification  of  the 
Sommerfeld  half-space  problem.  Sommerfeld  was  able  to  solve  for  the  electric 
fields  from  a  radiating  dipole  located  above  an  interface.  In  our  scattering 
problem,  the  dipole  is  driven  by  the  incident  plane  wave.  Rayleigh  was  able  to 
show  that  a  plane  wave  induces  a  dipole  moment  in  the  direction  of  the  inci¬ 
dent  electric  field.  We  consider  the  radiation  from  a  dipole  in  the  x-z  plane 
oriented  at  angle  3  from  the  z  axis.  The  dipole  moment  is 

p  »  pit  stn0  ♦  t  cos0)  Six)  Sly)  S(z-d).  11.1) 

Sommerfeld's  far-fieid  vector  potential  has  the  following  form 


n  =  n.t  ♦  n.t 


(1.2) 


s 


where 


ifcr  ikr 

n.  =  -r2—  sin0  [  —  +  Rt ,(W^t  1 
4ns.  L  r  r  J 


n,  = 


{_  'l«r  ihr  *r'  1 

cos(3|[—  -  Rtm(8)^t-  J  ♦  2sin3cosipsin^cosdRvt(&)^r  j 


and  r'  is  measured  from  a  source  point  located  a  distance  d  below  the 


r’  =  r  *  2dcos$ 


RreO.) 


cosfli  -  [  m,J  -  sin1^]171 
cosdi  «■  [mJ-  sinI8<]’/a 


Rtm(^) 


mJcosft  -  [mwa-  sina$t]t/a 
m«.r2cos$>  ♦  [m../—  sin2$i]*/2 


cos&i  -  [m« w-  sin^&i]'^ 
nw2cos8i  *  [m»-  sin3di],/a 


The  electric  fields  can  be  calculated  from 


B  *  k*n*  V  (V  ■  H) 


which  yield  the  following  electric  far-field  components 


(1.3) 

surface 

(1.4) 


(1.5) 


(1.6) 


(1.7) 


(1.8) 
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Er  =  0 


E&  =  -k3^^-  j|[l -Rru(0)exp(2ikdcos$)  J  sin$cos(3 

-  £l+RTu(8)exp(2ikdcos$)  JcosOcosqpsingJ 


(1.9) 


(1.10) 


E<p  =  “k  \t tu 


j~  £l  ♦  Rn($)exp(2ikdcos0)]  sinqpsinp. 


(1.11) 


A  dipole  moment  p  is  induced  in  the  direction  of  the  incident  electric  field 
with  a  magnitude 

p  =  4tc£.  ”***--—  aJE“*.  (1.12) 

m.*  ♦  2 

The  field  incident  on  the  small  sphere  has  two  sources.  One  source  is  the  in¬ 
cident  plane  wave  (including  the  plane  wave  which  reflects  off  the  surface). 
The  second  source  is  the  scattered  field  of  the  sphere-surface  system.  This 
scattered  field  is  proportional  to  the  induced  dipole  moment  which  is  propor¬ 
tional  to  the  sphere  radius  cubed,  and  at  the  sphere  center  is  also  proportio¬ 
nal  to  l/2kd.  This  term  is  very  small  compared  with  the  contribution  of  the 
incident  plane  wave  source.  We  therefore  neglect  this  scattered  term  and 

a 

assume  the  dipole  modes  of  the  sphere  are  driven  by  the  incident  plane  wave 
(including  the  plane  wave  which  reflects  off  the  surface).  In  sections  4  and  5 
we  consider  this  assumption  in  more  detail.  We  consider  two  incident  polarization 
states.  A  TE  plane  wave  will  induce  a  dipole  moment  in  the  sphere  given  by 


prt  -  p£l  ♦  RTt(rt-a)exp(-2ikdcosa)J^. 


(1.13) 
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A  TM  plane  wave  will  induce  a  dipole  moment  in  the  sphere  given  by 


pru  =  p  [l  +  RTu(7t-a)exp(-2ikdcosa)]8cosa 


-  p  Ql  —  Rru(7:-a)exp(r2ikdcosa)]£sina. 


The  scattering  amplitude  matrix  can  now  be  calculated  from  1.10  -  1.14: 


Si  =  Ai£l  ♦  RTt(it-a)exp(-2ikdcosa)J  £l  +  R™($)exp(2ikdcos$)]cos<p 

Sa  =  Ai{[[l  -  Riu(7t-a)exp(-2ikdcosa) J  £l  -  Rru(8)exp(2ikdcos8)Jsinasin8  + 

[]l  ♦  RTM(x-a)exp(-2ikdcosa)J  £l  ♦  Rru($)exp(2ikdcosft)Jcosacos9cos<pJ 
Sj  =  Ai[^l  +  Rrt(TC-a)exp(-2ikdcosa)J  £l  *  Rru($)exp(2ikdcosd)  J  cos&simp 


S«  =  -Ai|~l  ♦  RrM(it-a)exp(-2ikdcosa)J  [l  *  RTt(8)exp(2ikdcos$)Jsin<pcosa 


where 


A-  (k.)‘E- 

ikr  m>m  ♦  2 


(1.14) 

(1.15) 

(1.16) 

(1.17) 

(1.18) 

(1.19) 


These  four  components  represent  the  solution  to  the  scattering  problem  and 
will  be  compared  with  solutions  derived  from  the  other  theories. 
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2.  RavLeigh  Limit  of  sphere-surface  scattering 


The  second  method  we  examine  is  a  limiting  approach  of  the  light 
scattered  from  a  sphere-surface  system  as  the  sphere  becomes  optically  small. 
One  of  us  (Videen)  has  derived  an  analytical  solution  for  the  scattered  far 
fields  from  a  sphere-surface  system  (ref.  22).  The  orientation  of  the  scattering 
system  is  slightly  different  than  the  one  used  in  section  A.  Figure  2  shows 
the  scattering  system  used  to  predict  the  scatter  from  a  sphere-surface  sy¬ 
stem.  A  sphere  of  radius  a  is  suspended  a  distance  d  below  a  surface  separat¬ 
ing  two  media:  free  space  where  z  <  0  and  a  medium  of  complex  refractive 
index  nw  where  z  >  0.  The  results  may  be  expressed  as 


Si  s  22  (-i)"exp(im<p)|j]l  *lfc*(rc-d)(-l)"’"exp(-2ikdcosd)J  P.(cosd) 

«*0  «*■-»»  SlUv 

♦  [i  -  Rre(n-O)  (- 1) exp( -2ikdcosd) J  el-jP.lcosd)}  E~?Xp^kr)'  (2.1) 


00  n  r  TM 

S,*-i22  (-i)"exp(im«p)|[  1  -  RrU rc-9) (- 1) """ exp< -2ikdcosd) ]  2^2- R"(cosd) 


sind 


[l  ♦iW*-d)(-irmexp(-2ikdcos$)]  f^^-P^cosa) } E1" (2.2) 


00  m  r  tw 

Sl*- *22  (-i)*exp(im«p){ri-  Rru(7t-d)(-ir"exp<-2ikdcosd)l  2=S1  P.lcosd) 

lu  sind 


[l  ♦  Rri<(rc-d)(-l)"'mexp(-2ikdcosd)]  f£-~P."(cosd)}  E~ <2.3) 
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00  m  ,  _  r  Tm 

S.  -  22  <_i)"exp(im<p)|[l  +  Rte< Tr-^)(-t>" "'exp(-2ikdcos$)J  — ^  P.<cosd) 

*i«0  m*-« 

♦  [l  - Mk-V(-[) "•"’exp(-2ikdcosd)]  eir^P-lcosd)} E"c-^-k-  (2. 


4) 


where 

elm  =  [[1  -  Rrt(a)(-l)"’"exp(2ikdcosa)  ]a£ 


Rrt(O')  2  ("D""  [f»  »D!' 


Tt  _<n>0 

“  6nmCn 


]]q: 


(2.S) 


elm  =  [[1  -  Rru(a)(-l)"*"exp(2ikdcosa)]aIm 


♦  W0*)2<-1>"'"Cf-“D""'"1  :  e2c.w]]q:  <2.6) 

n 

fim  =  £[1  ♦  Rr€(a)(-l)"*”exp<2ikdcosa)]bII. 

♦  Rri<0*)  2  <-!)"  '”Cf-T*  Ci"-‘  -  evmDl"-”’  ]]q;  (2.7) 

i»  •hwl 


=  jj[l  ♦  Rr*»(a)(-l)"'"exp(2ikdcosa)]blli 


*  Rt»<0*>2  [f-T»Ci"”’ 


TU  _  ("V"1 
“  Sum  lAi 


]]Q" 


(2.8) 
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and  the  Lorenz-Mie  coefficients  are  represented  by 


ki>4t4*«  (k.x.a)iMka)  ~  k4«»»j/..(ka)4'«(k.»>'a) 

k.,i<n^-  (k..i-a)^n(ka)  -  kn.w£»(ka)4j"(k.*a) 


k.,h(i4)Mka)4)i.(k.>»a)  -  kti.^^'^k.^ai^.lka) 
k.^tx^’n  (ka)4i»(k.„a)  -  kp.*.^’"  (k.^.a)^n(ka) 


(2.10) 


2  i"'1  c>P^(cosa) 
n(n+l)  <3a 


(2.11) 


b£ 


2  i1"2  mP.(cosa) 
n(n+l)  sina 


(2.12) 


al“  =  ibl 


(2.13) 


b 


TM 

itm 


ia£ 


(2.14) 


— ,(«  ,m) 

L* 


(»>i  2kd  n  ♦  m»  1  (n'.ml 

Crt  “  C«*( 

2n  ♦  3  n  ♦  l 


2k d  n  -  m  <«>.» 

A  .  0.-1 

2n  -  1  n 


(2.15) 


D 


in', mi 


-2ikd 
n(n  ♦  1) 


m  c. 


(2.16) 


cT  =  -/Zn*i  h.m(2kd) 


(2.17) 


o!"'"  =  --/2n+l  h«'"(2kd) 


(2.18) 
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/(n'-m+l)(n’-*'m)(2n-*-l)  cl"  =  /(n-m+tHn+m)(2n*l)  g.  ' 


2kd  /*n-m'l'2}(n-m+l)  2^(1  / (n*n>)(n*m-t) 

y  <2n+3>  ~  /  (2n-l) 


(«  ,m- 1) 

C-.  (2.19) 


c«  =  c«  (2.20) 


(n‘+l)c. 


/  .v  l"'0' 

(n+1)  c«.» 


(2.21) 


where  and  £.(r)  are  the  Riccati-Bessel  functions  of  the  first  and  third 

kinds,  respectively,  and  the  primes  denote  their  derivatives  with  respect  to  the 

<N. 

argument.  The  functions  Pr(cosd)  are  the  normalized  associated  Legendre  poly- 

a 

nomials.  These  results  may  be  simplified  when  the  radius  of  the  sphere  be¬ 
comes  small  in  comparison  to  the  incident  wavelength.  This  can  be  taken  care 
of  by  expanding  the  Bessel  functions  for  small  arguments: 


™  ~  »'  r.~sl:.(a,.u  <2-22’ 

^‘"(p)  ~  -ip'1"*”  1  •  3  •  5 ...  (2n-l).  (2.23) 


With  the  sphere  much  smaller  than  X  we  now  can  consider  some  specific  ran¬ 
ges  for  the  sphere-surface  separation  distance  d.  First  we  consider  the  case 
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when  the  separation  distance  d  »  X.  In  this  case  the  fields  scattered  by  the 
small  sphere  (which  decay  as  l/2kd)  are  negligible  when  compared  with  the 
fields  of  the  incident  plane-wave.  The  interaction  which  occurs  between  the* 
sphere  and  image  sphere  are  analytically  represented  by  the  summation  terms 
of  equations  2.5-2.8.  We  now  consider  these  terms  in  detail.  The  translation 
coefficients  c,  " are  proportional  to  the  spherical  Hankel  functions  which 
decay  as  t/(2kd).  This  dependence  can  be  inferred  from  the  addition  theorem 
for  spherical  scalar  wave  functions. As  a  result,  for  large  values  of  d  the 
summation  terms  in  equations  2.S-2.8  are  negligible.  This  is  confirmed  numeri¬ 
cally  for  several  cases.  Equations  2.S-2.8  may  be  expressed  as 


elS,  =  [1  -  Rn(  a)  ( -D"  ”exp(2ikdcosa)  ]  aZ  Ql  (2.24) 

=  [1  -  Rru(a)(-l)"'”exp(2ikdcosa)]a^  Q".  (2.25) 

f2  =  [1  ♦  RTt(a)(-l)"  ’"exp(2ikdcosa)]b^  Q "  (2.26) 

flT  =  [1  +  Riu(a)(-l)"‘"exp(2ikdcosa)]  bill  Q",  (2.27) 


For  small  sphere  radius  a,  the  coefficients  Ql  and  Q J  are  proportional  to 
(k.^a)*"'  or  greater  powers  in  k„*a.  Because  of  the  dependence  of  the  coeffi¬ 
cients  on  sphere  radius,  only  the  first  term  (n  =  i)  is  significant.  This  is  equi¬ 
valent  to  making  an  electric  dipole  approximation: 


Ql  ~  0  (2.28) 


Q- 


m,»»  —  1 
-  mu*1  ♦  2  . 


(2.29) 
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From  2.28  and  2.29,  equations  2.1-2. 4  reduce  to 

Si  =  Ai[l  ♦  RTt(a)exp(2ikdcosa)]]  [l  ♦  RTi(jt-0)exp(-2ikdcos$)Jcos<p  (2.30) 

Sj  =  Aij[l  -  Rru(a)exp(2ikdcosa)J  []l  -  RTu(7r-$)exp(-2ikdcos8)^]sinasinO  + 

[j  ♦  Rru(a)exp(2ikdcosa)J  £l  +  RTu(7t-$)exp(-2ikdcos$)JcosacosOcos<p^.  (2.31) 


Sj  =  Ai[]l  ♦  Rrt(a)exp(2ikdcos«)J  [j  ♦  RTu(7t-0)exp(-2ikdcos$)Jcos$sin<p  (2.32) 


S«  =  -Ai£l  *  Rru(a)exp(2ikdcosa)J  £l  ♦  Rrc(7t-$)exp(-2ikdcos8)Jsin(pcosa.  (2.33) 


Except  for  the  system  orientation,  these  results  would  be  identical  to  those 
given  in  section  1  for  the  Sommerf eld- Rayleigh  system.  Recall  that  equations 
2.30-2.33  were  derived  assuming  a  large  sphere-surface  separation,  an  assump¬ 
tion  not  actually  necessary.  For  example,  the  coefficients  eZ  and  fZ  are  pro¬ 
portional  to  Q2  and  Q?,  respectively,  which  are  themselves  proportional  to 
(k.'i.a)2"'1.  For  equations  2.30-2.33  to  be  valid,  the  summation  terms  of  equations 
2.5-2.8  must  be  small  compared  with  the  plane  wave  expansion  coefficients 
and  bll  (given  by  2.11-2.14)  which  are  independent  of  the  sphere  radius.  All 
that  is  required  for  this  to  occur  is  for  the  translation  coefficients  to  be 
well-behaved.  The  translation  coefficients  are  a  function  of  the  spherical  Han- 
kel  functions  h.”'(2kd)  which  become  large  only  when  4rcd  «  X.  Therefore, 
equations  2.30-2.33  are  valid  except  when  the  sphere-surface  separation  d 
becomes  much  smaller  than  the  wavelength. 
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3.  EIT  formulation 

The  third  method  we  examine  is  the  exact-image  theory  (EIT)  deve¬ 
loped  by  Lindell  et  al  and  Muinonen  et  al.  The  time-dependence  used  in  these 
manuscripts  is  exp(iut),  so  some  minor  adjustments  will  be  needed.  Some 

minor  adjustments  were  made  to  their  equations  to  account  for  a  slightly  dif- 

#■  ' 

ferent  system  orientation.  Their  results  for  a  small  sphere  near  a  surface  may 
be  expressed  as 

EPi*  =  cAu  e*P---"--r)  ♦  t-QTfc.-E-'expiakdcosa) 

9  47tr  1 

♦  exp(-2ikdcos(M  +  $-^.*Q*^i*E’,"'exp(-2ikd(cosd-cosa))}  (3.1) 

E‘“  =  mV.  e-xP-~ikr-  f  ♦  f.Q-?t.*B“*  exp(2ikdcosa) 

”  \TZT  L 

♦  ■Q*E'"cexp(-2ikdcosd)  +  •Q*^<'E,",exp(-2ikd(co38-cosa))}  (3.2) 

The  bold-faced  cursive  characters  are  dyadic  functions  given  by 


jjC.  =  ii  Rru<7t-a)  ♦  yjRn(jt-a)  -  il  Rr«(Tt-a)  (3.3) 

»  [#•  Rr-(d)  ♦  $$Rn(7t-a)]*[  XX  ♦  JJ  -  ** ]  (3  4) 

Q  =  [  x*  ♦  jj  ]Q,(2d)  ♦  il  Q«(2d)  (3.5) 
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where 


Q,(2d)  =  - - j 


1  -  u>>BK.(2d) 


Q.(2d)  = 


t  ♦  o4uBIC(2d) 


«  .  nii»*i  -  1  i 

B  =  4  its. - : - —  a 

m.*2  +  2 


Mr,= 


nw2-  t  exp(-ikr)r  1  ♦  ikr  ♦  (ikr)2-i  8m, 


]  ♦  m  J  exp(-ik(z  -  ip/C)) 


2  ♦  1)  1  ♦  ik(z  -  ip/C)  i  v  ( m -  1  \  Ja-(p)  , 

'[ikil  -p/c7T  1 2,  “( JTi )  p  “p  (3-9' 


,  ,  „f  m(u,2  -  1  exp(-ikr)rl  *  ikr  ~|  8m.~2  f  °°  ,  ,  .  ,„Vk 

K-(r>"2l-„..i~-«»r  L~iSrtr~  ]  *  ^7^1  J.  exp<-'ku  -  Ip/ai 


f  1  ♦  ik(z  -  ip/C)  l  V1  /  nw2  -  1  \"  Ja»(p)  , 

t  -[Tkil  ^PcW )  4  "^STTl )  p  dp 


(3.10) 


C  =  k(tws.„  -  1)V2. 


Evaluating  the  dyadic  functions  yields  the  following  far- field  solutions 


A'lQ»(2d) 


[j  *  Rit(7C-a)exp(21kdcosa)J  £t  ♦  Rn($)exp(-2ikdcos&)  Jcosp  (3.12) 


Sa  -  {[[l  -  RT«(x-a)exp(2ikdcosot)J  £l  -  RrH<8)exp(-2ikdcos8)jQt(2d)sinasind 

[]l  ♦  Rru(it-a)exp(2ikdcosa)J£t  ♦  Rrn(8)exp(-2ikdcos8)^Q.(2d)cosacos8cos9j  (3.1 
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Sl=  Ai£|<2d>[i  ♦  Rn(rt-a)exp(2ikdcosa) J£l  ♦  RTu(d)exp(-2ikdcos$)  J  cos$sinq?  (3.14) 
S»  =  —  £l  +  RrM(rt-a)exp(2ikdcosa)JQ  +  Rti(d)exp(-2ikdcos&)  Jsinqpcosa  (3.15) 


_  exp(-ikr)  m„»  -  1 
ikr  m.*2  ♦  2 


(ka)JE" 


Except  for  the  time  dependence  and  factors  of  Qt(2d)/B  and  Q<(2d)/B,  these 
solutions  would  be  identical  to  those  given  in  the  previous  sections.  So  let  us 
examine  these  factors  a  little  more  closely.  These  factors  depend  on  the  func¬ 
tions  Kt  and  K.  (given  by  equation  3.9  and  3.10)  which  are  a  result  of  the 
scattered  fields  interacting  with  the  small  sphere.  Recall  that  in  the  previous 
sections  we  assumed  that  this  contribution  was  negligible  so  that  Qi(2d)  ~ 
Q.(2d)  ->•  B.  We  now  examine  this  approximation  in  more  detail.  The  contri¬ 
butions  of  K«  and  K>  will  be  greatest  when  the  permittivity  of  the  medium  be¬ 
low  the  surface  approaches  infinity.  In  this  case,  3.9  and  3.10  simplify  to  / 


Equations  3.6  and  3.7  may  now  be  expressed  as 

Q,(2d)  *  B  [  1  -  exp(-2ikd)  [l  *  2ikd  ♦  (2ikd)*  ][ ggi-~-]  ]  '•  (3.19) 
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Q.(2d)  =  B 


l  -  2exp(-2ikd)  [l  + 


2ikd][ 


2  4 

TTltph  ~  I 
+  2 


-! 


(3.20) 


The  derivation  following  from  3.8  has  two  requirements:  the  sphere  radius  is 
much  smaller  than  the  incident  wavelength  (a  «  X)  and  the  sphere  radius  is 
much  smaller  than  the  wavelength  inside  the  sphere  (m.*.a  «  X).  For  sphere- 
surface  separations  on  the  order  of  X  or  greater,  we  have  Q<(2d)  ~  Qt(Jd)  ~  B. 
This  is  a  direct  result  of  the  first  requirement,  in  which,  a  «  d.  We  now  con¬ 
sider  a  worst  case  scenario  for  our  approximation.  We  note  that  if  the 
sphere-surface  separation  is  equal  to  the  sphere  radius  these  equations  further 
simplify  to 

<*<*»  - » [  i  - 

q.<2“>  ■  B  [  1  - 

The  second  requirement,  that  m(»a  «  X,  helps  ensure  that  we  can  assume 
Q>(2d)  ~  Q.(2d)  ~  B,  saying  that  the  sphere  refractive  index  cannot  be  large. 
For  a  numerical  example,  this  assumption  is  good  to  within  10%  for  a  lossless 
dielectric  sphere  of  refractive  index  m.*.  =  1.6.  Although  this  accuracy  does  not 
seem  outstanding  we  should  remember  that  this  is  the  worst  case  in  which 
the  medium  below  the  surface  is  a  perfect  conductor  and  the  sphere-surface 
separation  Is  much  smaller  than  the  wavelength. 


(3.21) 

(3.22) 
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4.  Projection  model 

The  three  theoretical  methods  thus  far  discussed  all  yield  the  same 
results.  We  now  discuss  one  more  approach  which  is  particularly  insightful. 
The  method  is  based  on  the  double  interaction  model  of  Nahm  and  Wolfe  and. 
dipole  scattering  discussed  by  van  de  Hulst.  In  the  double  interaction  model, 
the  scattered  light  from  a  sphere  near  a  surface  can  reach  a  detector  *hy  one 
of  four  paths.  It  may  1.)  strike  the  sphere  directly  and  scatter  directly  to  the 
detector;  2.)  strike  the  sphere  directly,  scatter  to  the  surface  which  reflects 
light  to  the  detector;  3.)  reflect  off  the  surface,  strike  the  sphere  and  scatter 
to  the  detector;  or  4.)  reflect  off  the  surface,  strike  the  sphere,  scatter  to  the 
surface  which  reflects  light  to  the  detector  (see  Fig.  3).  The  light  scattered  by 
the  sphere,  reflecting  off  the  surface  and  striking  the  dipole  is  assumed  to  be 
very  small  compared  with  the  incident  fields  and  contribute  only  negligibly  to 
the  scattered  fields.  Assuming  the  small  sphere  is  a  dipole,  the  magnitude  of 
the  scattered  electric  field  is  proportional  to  the  dipole  moment  projected 
onto  a  plane  perpendicular  to  the  radius  vector.  Although  it  is  not  necessary, 
for  the  purpose  of  this  discussion,  we  restrict  our  detector  to  the  x-z  plane 
(<p  -  0*).  We  sum  the  four  components  of  the  TE  mode  that  make  up  matrix 
element  Si  as 

Si  a  Al£l  ♦  Ri«(re-a)exp(-2ikdcosa)  ♦  RnU»exp(2ikdcos&) 

♦  Rn(tt-a)Rn(8)exp(2ikd(cos$-cosa))] 


=  Ai£l  ♦  Rr«(jt-a)exp<-2ikdcosa)]  [l  ♦  Rn($)exp(2ikdcosd)] 


(4.1) 
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Similarly,  we  sum  the  four  components  of  the  TM  mode  that  make  up  element 
Sa  as 

Sa  =  Ai[]cos(d-a)  -  Rru(rc-a)cos(rt-d-a)exp(-2ikdcosa)  -  RTM($)cos(7t-$-a)exp(2ikdcos$) 

♦  cos(D-at)  Rru(x-a)  Rtu(D)exp(2ikd(cos$-cosa) )] 

■  "{['  -  R.Tu(rt-a)exp(-2ikdcosa)J  j^l  -  RruO)exp(2ikdcos^) JsinotsinD  + 

£l  ♦  Rm(jr-a)exp(-2ikdcosa)J  £l  +  RruU»exp(2ikdcos&) J  cosacos$  |  .  (4.2) 

Since  <p  =  0*.  elements  Si  =  &  =  0.  Although  this  solution  is  derived  using  a 
very  simple  model,  the  results  are  identical  to  equations  1.15  -  1.19  (<p  =  0”). 


Multiple  scattering 


interactions 


We  derived  the  first-order  terms  in  the  scattering  amplitude  matrix 
in  the  first  section  in  equations  1.15-1.19  where  order  refers  to  the  number  of 
interactions  that  a  beam  of  light  undergoes  with  the  small  sphere  before  being 
scattered  to  the  far  field.  The  order-of-scattering  approach  can  be  used  to 
add  correcting  terms  to  these  results.  Using  this  approach,  the  scattered  field 
fron  one  portion  of  the  system  is  considered  to  be  an  incident  field  on  an¬ 
other  portion  of  the  system.  Then  the  scattered  field  from  this  portion  of  the 
system  is  calculated  and  the  process  is  repeated  until  the  desired  accuracy  is 
achieved.  To  perform  a  rigorous  order-of-scattering  analysis  requires  perform- 
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ing  translations  of  the  scattered  fields  and  the  resulting  equations  are  similar 
to  those  in  section  2,  containing  the  same  translation  coefficients.  Such  a 
derivation  is  not  our  intent.  Insight  into  the  problem  can  be  gained  by  consi¬ 
dering  the  interaction  which  occurs  between  the  scattered  dipole  radiation 
terms  which  are  spherical  waves  reflecting  off  the  plane  surface  and  interact¬ 
ing  with  the  dipole.  Since  the  sphere  is  much  smaller  than  the  wavelength,  we 
can  assume  the  field  to  be  approximately  constant  over  the  volume  of  the 
sphere  and  avoid  the  cumbersome  translation.  It  is  relatively  straightforward 
to  calculate  the  effect  of  the  second  order  terms.  This  case  includes  the 
fields  scattered  from  the  small  sphere,  striking  the  plane  surface  and  inducing 
a  dipole  moment  in  the  small  sphere.  From  equations  1.12-1.19  we  get 

W  3  4tcs,[  ~  2  JkV  ""^kd"  E1  +  Rji(it-a)exp(-2ikdcosa)]  Rr.(0)  E~  f 

(5.1) 


for  the  TE  case  and 


pruj  =  4 ret.  £  +  2  ]  k1*"  **^kd^  E*  *  Rn*(7l~at)exp(-21kdcosot)]  Rrv.(0)  E*"  #cosa 

(5.2) 


for  the  TM  case.  Since  the  scattered  fields  are  proportional  to  the  dipole  mo¬ 
ment  induced  In  the  small  sphere,  it  is  apparent  that  the  second  order  dipole 
moment  will  have  only  a  negligible  effect  since  it  is  proportional  to  kV/d 
times  the  first  order  dipole  moment. 
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We  can  include  the  effects  of  all  the  higher  order  scattering  inter¬ 
action  terms  since  the  induced  dipole  moment  is  due  to  the  incident  field  and 
the  scattered  field  from  the  small  sphere  (from  equations  1.12-1.19): 

pra  =  p£l  ♦  RTf(rc-a)exp(-2ikdcosa)J^  ♦  p™  ~  kV  eX^jt£j  Rn(0)  (5.3) 


pru  =  p[l  ♦  Rru( ic-<x)exp( -2ikdcosa) ]  ficosa  -  p  [l  -  Rr«<TC-a)exp(-2ikdcosa)]*sina 


pTH* 


-  L  kV  *XP'2Vt-dj  Rm<0) 


m.M.2  ♦  2 


2kd 


(5.4) 


where  piu.x  is  the  x  component  of  prw.  These  two  equations  yield  dipole  mo¬ 
ments  of  the  form 


p™  =  fl  -  rn,f*a  ~  —  kV  Rn(0)  1  p[l  ♦  Rit (ir-a)exp(-2ikdcosa) ] f 

L  m«*»  ♦  2  Zko  J 

(5.5) 

p™  =  -L  kV  Rru(O)  J  p[l  *  Rm(jt-a)exp(-2ikdcoso) J  ficosa 


-  p  [l  -  Rnk(rc-at)exp(-2ikdcosa)Jfsina  (5.6) 


Note  that  the  first  terms  in  the  expansions  of  equations  5.5  and  5.6  are  the 
second  order  correction  terms  of  equations  5.1  and  5.2.  The  scattering  coeffi¬ 
cients  can  be  expressed  as 


Si  ~  LAC^l  ♦  Rr«(rc-<x)exp(2ikdcosa)]  [l  ♦  Rrt(d)exp(-2ikdcosd)Jcos<p 


(5.7) 
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Si  =  iA{c[l  -  Riw(7t-a)exp(2ikdcosa)]  £l  -  Rru(O)exp(-2ikdcos0) J  sinasind  ♦ 

[]l  ♦  RrM(rt-a)exp(2ikdcosa)J[^l  ♦  Riu(d)exp(-2ikdcosd)  Jcosacos-&cos<p|  (S.f 


Sj  =  iAC[l  ♦  Rrt(n-a)exp(2ikdcosa)J[]l  ♦  Rru($)exp(-2ikdcosO)  JcosOsin<p  .  (5.9) 


S«  =  -iAC^l  ♦  Rru(7c-a)exp(2ikdcosa)J  £l  ♦  Rn(8)exp(-2ikdcosO)J  sinqpcosa  (5.10) 


m.*. 1  -  1  .  j  i  exp(21kd) 
m.*J  >2  2kd 


Rn(0) 


(5.11) 


since  Rn(0)  =  Rru(O).  Note  that  the  correction  term  C  varies  only  slightly  from 
1  since  kV/d  «  1  .  This  appears  to  contradict  the  analysis  of  section  3  where 
the  scattered  fields  could  vary  by  as  much  as  10X.  It  should  be  noted  that  the 
analysis  in  section  3  was  for  a  very  small  separation  distance  d;  whereas  the 
analysis  in  this  section  was  derived  from  the  radiation  terms  of  the  dipole 
scatter  and  is  invalid  for  small  separations  d.  In  the  near-fleld  scatter  from  a 
small  sphere,  there  are  scattering  terms  proportional  to  (kr)'J  which  become 
significant  along  with  a  radial  component.  These  terms  were  not  considered  in 
this  analysis. 
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Results 


Up  to  this  point  we  have  derived  the  far-field  electric  fields  for  a 
small  sphere  near  an  interface  using  four  different  formal  solutions.  We  now 
examine  the  results  of  these  derivations  which  are  expressed  by  equations 
1.1S-1.19.  Rather  than  examine  just  the  intensities  of  the  TE  and  TM  compo¬ 
nents,  we  will  put  the  results  in  Mueller  matrix  form,  since  the  16-element 
Mueller  matrix  contains  the  phase  information  of  the  scattered  fields.  We  will 
also  restrict  ourselves  to  the  plane  of  incidence  (9  =  0*).  The  particular  Muel¬ 
ler  matrix  formulation  we  use  is  given  by 

&<  =  J  [  IS, I*  ♦  I  Sal2  ], 

Sa  *  \  [  IS,!2  -  I  Sal2  ]/S„, 

Sia  =  Re(SiSa*)/S„,  and 

Sm  a  -Im(SaS,*)/S«.  (S.l) 

With  the  detector  in  the  plane  of  incidence  (9  -  0*),  the  other  scattering  S.j(0)  are 
either  zero  or  simple  multiples  of  these  elements.  It  should  be  noted  that  the 
results  derived  do  not  include  the  specular  component  of  the  scattered  field. 


The  first  set  of  scattering  S<j(&)  we  consider  are  shown  in  Figure  4. 
The  system  is  illuminated  at  normal  incidence  (a  »  180*)  and  optical  constants 
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of  the  medium  below  the  plane  are  chosen  to  correspond  with  an  aluminum 
surface  <nw  =  0.5  ♦  5.0i  at  X  =  0.4416pm).  Scattering  Sij(9)  are  shown  for  four 
different  sphere- surface  distances  d.  Note  that  when  the  system  is  illuminated 
at  normal  incidence,  equations  4.1  and  4.2  reduce  to  two  periodic  functions  in 
d.  One  of  the  functions  is  independent  of  9  and  determines  the  maximum 
amplitude  of  the  signal.  This  amplitude  cycles  through  maximum  and  minimum 
values  with  a  period  of  d  =  X/2.  If  the  medium  below  the  surface  is  a  lossless 
dielectric,  then  Rr«(a)  is  real  and  amplitude  minima  occur  when  d  =  nX/2  and 
amplitude  maxima  occur  when  d  =  X/4  +  nX/2.  If  the  medium  below  the  sur¬ 
face  is  lossy,  then  R Mat)  is  complex  and  a  phase  difference  will  be  introduced. 
The  parameters  chosen  in  Fig.  4  introduce  only  a  slight  phase  difference. 
Matrix  element  Sh(9)  shows  that  the  highest  signal  intensities  (of  the  four 
curves  shown)  occur  when  d  =  X/4  and  the  lowest  signal  intensities  occur 
when  d  *  X/2. 

The  oscillatory  structure  in  the  matrix  elements  of  Fig.  4  is  also 
periodic  in  d  and  is  very  nearly  periodic  in  9.  From  equations  4.1  and  4.2,  we 
would  expect  the  total  intensity,  Sn(9),  to  be  at  minima  values  approximately 
when  2dcos9  *  nX  (where  n  is  an  integer)  for  a  dielectric  medium  below  the 
surface.  For  the  d  *  X/2  curves  of  Fig.  4,  we  would  expect  a  minimum  to 
occur  at  9  3  0*.  The  minima  have  been  displaced  slightly  because  the  medium 
below  the  substrate  is  slightly  conducting.  Note  that  spikes  and  inflection 
points  occur  in  the  polarization  matrix  elements  at  the  angular  locations  where 
the  intensity  signals  are  at  minima. 
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Another  set  of  matrix  elements  is  shown  in  Figure  5  for  the  same 
systems  as  Fig.  4,  but  illuminated  at  near-grazing  a  =  101*.  In  Fig.  4,  where 
the  illumination  was  at  normal  incidence,  the  sets  of  matrix  elements  were 
very  similar.  Changing  the  incident  angle  away  from  normal  incidence  breaks 
the  symmetry  and  increases  the  complexity  of  the  matrix  elements. 

Conclusion 


Four  theoretical  models  were  examined  and  found  to  predict  the  • 
same  scattering  signals  from  a  small  sphere  suspended  above  a  plane  surface 
interface.  The  approximation  breaks  down  when  the  sphere-surface  separation 
also  becomes  much  smaller  than  the  wavelength.  In  this  case  the  maximum 
error  in  the  approximation  is  on  the  order  of  10X  of  the  theoretical  value. 
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Figure  Captions 

Fig.  1.  The  geometry  of  the  scattering  system.  A  sphere  of  radius  a  is  located 
a  distance  d  above  a  surface.  A  plane  wave  travels  in  the  x-z  plane  at 
angle  a  with  respect  to  the  z  axis. 

Fig.  2.  The  geometry  of  the  scattering  system.  A  sphere  of  radius  a  is  located 
a  distance  d  below  a  surface.  A  plane  wave  wave's  in  the  x-z  plane  at 
angle  a  with  respect  to  the  z  axis. 

Fig.  3.  The  geometry  of  the  scattering  system  showing  the  four  paths  light 
may  travel  to  reach  the  detector  in  the  double  interaction  model. 

Fig.  4.  The  Mueller  matrix  elements  from  a  small  sphere  at  four  separations 
from  a  surface  of  index  n  =  O.S  ♦  5.01  illuminated  at  normal  incidence 

(a  =  180',  A  =  1).  * - *  d  =  X/16 .  O - O  d  =  X/8.  X - X  d  =  X/4. 

- - -  d  =  X/2. 

Fig.  5.  The  Mueller  matrix  elements  from  a  small  sphere  at  four  separations 
from  a  surface  of  index  n  =  0.5  ♦  S.Oi  illuminated  at  at  a  101*.  (A  =  1). 


d  =  X/16.  O - O  d  =  X/8.  X - X  d  =  X/4. 


d  =  X/2. 
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Stokes  vectors,  Mueller  matrices  and  polarized  scattered  light:  experimental 
applications  to  optical  surfaces  and  all  other  scatterers 


William  S.  Bickei  and  Gorden  Videen 

University  of  Arizona,  Department  of  Physics 
Tucson,  Arizona  85721 


ABSTRACT 


We  discuss  scattering  in  the  context  of  the  Stokes  vectors  and 
Mueller  matrices  that  characterize  the  interaction.  In  order  to  study  surface 
structures  using  light-scattering  techniques  it  is  useful  to  examine  the  nature 
of  light  scattered  from  perfect  and  perturbed  mirror  surfaces. 

L  INTRODUCTION 

A  highly  motivated,  systematic,  and  fundamental  approach  to  surface 
scattering  requires  that  the  initial  surfaces  be  fundamental  and  that  the  con¬ 
tamination  to  produce  the  surface  scattering  be  known  and  controlled.  With 
this  in  mind  any  rough  surface  can  be  considered  to  be  a  perturbed  perfect 
surface  that  has  reached  its  final  condition  through  some  continuum  contami¬ 
nation  process. 

The  powerful  Stokes  vector,  Mueller  matrix  light-scattering  tech¬ 
niques  can  be  applied  directly  to  study  surfaces.  Scattering  from  perfect  sur¬ 
faces  can  be  theoretically  predicted  and  experimentally  measured.  Experimen¬ 
talists  have  an  advantage  in  that  they  can  measure  what  theorists  cannot  cal¬ 
culate.  If  the  data  are  to  be  used  to  check  theory,  the  experimentalist  must 
relate  the  light  scattering  signal  to  the  exact  structure  and  orientation  of  the 
surface  that  scattered  the  light. 

Consider  how  a  perfect  optical  surface  (lens)  is  treated  by  taking 
into  account  its  geometrical  and  diffraction-limited  properties.  Geometrical 
optics  is  used  to  predict  where  a  point  on  the  ob)ect  will  be  focused  to  a 
point  on  the  image.  Diffraction  theory  will  predict  how  the  image  point  is 
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actually  a  diffraction  pattern;  i.  e.,  not  all  light  from  the  object  point  ends  up 
at  the  image  point.  Nevertheless,  the  diffraction-limited  image  is  exactly  pre¬ 
dicted  by  theory  and  itj  intensity  distribution  is  related  to  fundamental  con¬ 
stants.  The  problem  arises  when  the  perfect  diffraction-limited  optic  and  its 
image  is  perturbed  by  a  defect  in  or  on  the  optic.  This  defect  scatters  light 
into  all  directions  —  out  of  the  paths  so  well  defined  by  geometrical  and 
diffraction  optics. 

The  question  is  where  does  this  scattered  light  go  and  what  does  it 
do  to  the  image  (information).  There  will  be  a  loss  of  image  intensity,  but 
more  troublesome  is  the  light  initially  intended  for  one  part  of  the  image  that 
is  scattered  to  another  part  of  the  image  causing  loss  of  contrast  and 'image 
definition.  We  can  write  that  the  object  (point  A)  is  transformed  into  a  diffuse 
diffraction-limited  geometrical  image  (point  A‘)  by  a  transformation  matrix  CL, 
D,  SI  where  L  is  the  geometrical  function  of  the  lens.  D  is  the  diffraction 
function  and  S  is  the  scattering  function  which  contains  all  of  the  scattering 
parameters  of  the  defect. 

2,  POLARIZED  LIGHT  SCATTERING 

What  is  S?  If  the  defect  were  a  perfect  sphere  the  scatter  from  it 
is  nothing  more  than  the  diffraction  by  a  spherical  object  —  a  problem  solved 
by  Gustav  Mie  in  1908.'  As  the  sphere  becomes  more  irregular  in  shape  and  in¬ 
homogeneous  in  optical  constants,  exactly  solveable  sphere  diffraction  goes 
over  to  the  statistical  aspects  of  irregular  particle  scatter.1  A  scatterer  is  said 
to  be  characterized  by  its  scattering  properties  when  we  know  how  it  will 
rearrange  the  properties  of  light  incident  on  it.  For  example,  consider  perfect 
sphere  scattering.  A  sphere  with  radius  r,  refractive  index  n  and  aborption  4 
has  a  scattering  matrix  [S3.  When  it  is  illuminated  with  light  of  intensity  I., 
wavelength  X*  polarization  IL,  at  angle  8.,  it  creates  a  scattered  field  of  I.,  X., 
IT.,  and  8..  We  can  write  (I„  X.,  n.,  8.)  =  [SKI.,  X.,  IL,  8.)  where  [S3  is  exactly 
known  for  spheres.  The  matrix  CS3  can  be  calculated  only  for  certain  highly 
symmetric  cases.  It  can  be  approximated  for  slightly  distorted  systems,  but 
can  only  be  measured  for  truly  irregular  particles.  Therefore,  in  order  to  learn 
how  a  defect  on  or  in  a  surface  (an  irregular,  imperfect  surface)  will  affect  all 
aspects  of  the  incident  light,  it  is  better  to  measure  what  happens  exactly  in¬ 
stead  of  relying  on  approximation  and  theoretical  models.  Generally  two  things 
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are  required  of  scattering  experiments:  1.  produce  a  defect  and  measure  its 
scatter;  or  2.  measure  the  scatter  and  predict  the  defect  —  the  inversion  pro¬ 
cess.  Before  we  extract  information  from  light-scattering  data  it  is  important 
to  be  sure  that  we  have  all  the  information  available. 

T  POLARIZED  SCATTERING  NEPHELOMETER 

The  experimental  set-up  that  can  measure  the  polarized  intensities 
scattered  by  the  scatterer  CS]  into  the  angles  &  and  9  is  shown  in  Fig.  1.  The 
input  optics  can  be  selected  to  be  an  open  hole  Col,  or  horizontal  lirffear  po¬ 
larizer  Ch3.  or  ♦  45*  linear  polarizer  C*1  or  a  right-handed  circular  polarizer  Crl. 
The  exit  optics  choices  are  the  same  and  can  be  chosen  independently  of  the 
input  optics  and  can  swing  with  the  detector  through  the  scattering  angle  $ 
from  0*  to  180*.1 


Fig.  1.  Arrangement  of  entrance-exit  polarizers  to  measure  the  matrix  elements 
of  a  scatterer. 

The  following  example  shows  how  to  determine  v.  ..at  scattering  ma¬ 
trix  elements  Si  are  involved  when  a  particular  set  of  input-output  polarizers 
are  used  to  prepare  and  analyze  the  scattered  light.*  For  our  first  choice  we 
assume  that  the  arbitrary  scatterer  (defect)  is  illuminated  with  horizontally- 
polarized  light  (hi.  The  scattered  Stokes  vector  will  be  [V.]  =  CS]»th]  or 
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We  see  that  the  scatterer  CS]  mixes  the  initial  pure  polarization  state  Ch J  to 
produce  a  scattered  Stokes  vector  w-th  mixed  polarizations.  The  first  compo¬ 
nent  (S«  *  Sa)  is  the  total  intensity  and  is  a  sum  of  two  matrix  elements.  If 
this  scattered  light  is  now  passed  through  a  *45*  linear  polarizer  O]  we  get 
[V.]  =  MCV.J.  Specifically  we  have 
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The  first  component  of  the  final  Stokes  vector  is  now  a  mixture  of  four  ma¬ 
trix  elements.  The  element  sum  (S«  ♦  Sa  ♦  Sit  ♦  Sa)  is  the  total  intensity  that 
will  be  measured  by  the  detector.  We  put  the  result  of  all  such  calculations 
from  all  16  Stokes  vector  combinations  into  a  final  array  shown  in  Fig.  2.  The 
main  point  of  Fig.  2  is  that  there  are  *  X  4  *  16  measurements  to  be  made  to 
completely  characterize  the  polarized  light  scattered  from  the  defect.  These 
measurements  can  be  routinely  made  by  nephelometers  that  incorporate  the 
various  polarizers  in  their  entrance-exit  optics.  The  16  matrix  elements  are 
needed  to  completely  characterize  the  scattered  light  or  to  use  in  the  inversion 
process  to  determine  the  properties  of  the  scatterer. 
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Fig.  2.  Matrix  array  showing  matrix  element  combinations  measured  with  var¬ 
ious  arrangements  of  input  and  output  polarizers. 
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_L  SURFACE  ORIENTATION 

The  geometrical  orientation  of  the  rough  surface  with  respect  to  the 
entrance-exit  beams  must  be  addressed  early  because  there  are  an  infinite 
number  of  possible  orientations,  all  of  which  might  give  a  different  S,  for  the 
same  surface.  Figure  3  shows  some  of  the  geometrical  parameters  involved  in 
surface  scattering  studies.  The  laser  beam,  after  preparation  into  a  definite 
polarization  state,  strikes  the  surface  at  angle  a  and  is  scattered  into  all  4it. 
The  light  scattered  into  angle  $  is  analyzed  with  the  polarization  exit  optics 
and  detected  by  the  photomultiplier.  The  angle  t  is  a  surface  tilt  measured  in 
the  surface  plane.  Our  work  has  shown  that  geometry  is  important  and  that 
no  universal  scaling  factor  exists  and  that  no  best  orientation  exists  for  all 
cases. 


Fig.  3.  Optical  and  sample  arrangement  for  surface  scattering  measurements. 


S.  SURFACE  SCA3 
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Wa  now  show  some  general  results  of  scattering  from  some  sur¬ 
faces  to  show  how  light-scattering  nephelometry  and  BRDF  are  related.5*  Fi¬ 
gure  4  shows  four  matrix  elements  for  a  reflective  aluminum  surface.  The  4 
order  of  magnitude  angular  decrease  of  the  total  intensity  matrix  element  S« 
is  an  indication  of  the  quality  of  the  reflecting  surface.  Note  that  even  though 
the  total  scattered  intensity  is  down  by  over  S  orders  of  magnitude  near  S  - 
ISO*  (back-scatter)  the  polarization  of  the  scattered  light  is  significant  as  de¬ 
monstrated  by  matrix  Su  and  Sm. 
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Fig.  4.  Four  Mueller  matrix  elements  for  a  reflecting  aluminum  surface. 


Figure  S  shows  the  matrix  elements  for  the  same  surface  but  now 
roughened  to  saturation.  Further  roughening  will  not  change  the  surface  cha¬ 
racter.  All  hints  of  the  location  of  the  specular  peak  are  gone.  S«  as  well  as 
Sm  are  independent  of  surface  orientation;  and  polarizations,  as  indicated  by  all 
matrix  elements,  remain  large. 


Fig.  5.  Four  Mueller  matrix  elements  for  a  saturated  rough  aluminum  surface. 
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Now  we  show  how  dramatically  the  S.,  from  surface  defects  depend 
on  the  angle  at  which  the  surface  (and  its  defect)  is  illuminated.  Figure  6 
shows  the  matrix  element  Su  for  a  rectangular  line  (h  -  0.46pm,  w  *  1.10pm) 
on  a  smooth  aluminum  surface  illuminated  at  near-grazing  incidence  and  at 
normal  incidence.  Normal  incidence  not  only  restricts  the  angular  range  over 
which  the  data  is  received  but  it  also  wipes  out  all  the  phase  information 
needed  to  characterize  the  defect. 


Fig.  6.  Matrix  element  Su  for  a  rectangular  line  on  a  mirror  surface  measured 
at  two  different  angles  of  incidence. 

8.  CONCLUSIONS 

One  general  goal  of  light  scattering  is  to  develop  an  algorithm  to 
put  into  a  computer  which  will  predict  either  the  scattered  field  (from  particle 
properties)  or  particle  properties  (from  the  scattered  field).  The  data  (particle 
properties  or  scattered  field)  must  be  determined  experimentally.  We  must  de¬ 
termine  how  much  scattering  data  is  needed,  how  good  it  must  be  and  how 
well  it  describes  the  scatterer  in  a  practical  way.  We  also  must  determine  the 
amount  of  information  contained  in  the  various  matrix  elements  and  whether 
signal  changes  can  be  related  to  changes  in  the  optical  or  geometrical  proper¬ 
ties.  We  see  that  polar  nephelometry  gives  sixteen  matrix  element  signals.  The 
record  shows  that  polar  nephelometry,  eilipsometry,  BRDF,  and  other  optical 
techniques  can  all  complement  each  other  and  yield  important  information 
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about  surface  scattering  when  uoed  with  care,  keeping  in  mind  their  limitations 
and  range  of  validity. 
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The  light-scattering  Mueller  matrix  for  a  surface  contaminated  by 
a  single  particle  in  the  Rayleigh  limit 
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Abstract 

A  ray- tracing  model  was  used  to  derive  the  light-scattering  Mueller 
matrix  element  curves  for  a  dipole  near  a  perfect  surface  as  a  function  of 
angle  of  Incidence,  scattering  angle  and  surface  index  of  refraction.  This  sy¬ 
stem  represents  a  fundamental  system  composed  of  a  perfect  plane  surface 
and  the  perfect  (Rayleigh)  scatterer. 

Subject  terms i  dipole,  contamination,  light-scattering,  Mueller  matrix,  surface. 
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Mle  theory,'  which  can  be  used  to  predict  exactly  the  scattering 
from  spheres  of  all  types,  has  recently  been  the  subject  of  a  great  deal  of 
modification.  The  modification  is  brought  about  by  examining  the  simple 
sphere  when  it  comes  into  contact  with  the  perfect  surface.  The  impetus  has 
come  from  contamination  studies,  often  centered  around  wafer-producing  tech¬ 
nology  in  the  electronic  and  computer  sciences,  but  the  applications  are  nume¬ 
rous,  ranging  from  degradation  of  optical  modulation  transfer  functions  to 
energy  losses  in  laser  systems. 

The  modification  to  Mle  theory  investigates  how  the  scattering  sy¬ 
stem,  composed  of  a  single  sphere,  is  changed  when  the  system  Is  composed 
of  a  sphere  and  a  surface.  Since  each  system  taken  independently  has  been 
solved  exactly,  the  sphere  by  Mie  and  the  surface  by  Fresnel,  a  logical  first 
order  solution  would  be  to  combine  these  two  systems  into  something  similar 
to  what  Nahm  and  Wolfe*  refer  to  as  the  double-interaction  model.  This  model 
ignores  the  interaction  term  between  the  sphere  and  the  surface.  For  Rayleigh 
scattering,  the  scattering  intensity  is  proportional  to  (r/X)\  and  the  effect  of 
the  Interaction  term  when  compared  to  the  other  terms  is  negligible.  There¬ 
fore,  when  the  sphere  is  driven  to  the  Rayleigh  limit,  r  «  X,  this  model 
should  be  able  to  predict  the  scattering  exactly. 

THE.QRX 

Theoretically,  we  treat  the  system  as  a  dipole  illuminated  by  two 
plane  wavea,  one  incident  and  one  reflected  off  the  surface.  The  scattered  ra¬ 
diation  is  elan  in  two  parts,  one  directly  from  the  dipole  and  one  reflected  off 
the  surface.  Figure  I  shows  the  various  paths  by  which  light  can  be  scattered 
into  the  detector  and  the  convention  used  for  measuring  angles.  All  angles  are 
measured  from  the  normal  to  the  surface.  Figure  2  shows  the  four  paths  that 
a  beam  may  travel  for  scattering  on  the  incident  side  of  the  detector.  Path  I 
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is  the  direct  path  and  Is  what  comes  from  Mie  theory.  The  beam  whose  angle 

of  incidence  with  the  surface  normal  is  3.  excites  the  dipole,  which  scatters 

radiation  through  angle  d  to  the  detector.  Path  II  shows  the  beam  reflecting 

off  the  surface  before  striking  the  dipole.  Radiation  is  then  scattered  directly 

to  the  detector.  Path  Ill  shows  the  beam  striking  the  dipole  directly,  with  the 

scattered  radiation  striking  the  surface  before  reaching  the  detector.  Path  IV 

shows  the  beam  reflecting  off  the  surface  before  striking  the  dipole  whose 

scattered  radiation  also  reflects  off  the  mirror  before  reaching  the  detector. 

* 

Before  looking  at  the  polarizations,  we  will  complete  this  analysis 
by  looking  at  the  scattering  that  penetrates  the  surfaoe,  where  cos&  >  90*. 
Although  most  of  the  Interest  seems  to  be  in  the  case  where  the  surface  is 
highly  reflecting,  the  transmission  case  la  such  a  simple  extension  that  it 
should  not  be  neglected.  Figure  3  shows  the  two  paths  that  a  beam  may  travel 
in  reaching  a  detector.  Path  V  shows  the  beam  striking  the  dipole  directly 
with  the  scattered  radiation  being  refracted  by  the  surfaoe  to  the  detector. 
Path  VI  shows  the  beam  reflecting  off  the  surface  before  striking  the  dipole, 
with  the  scattered  radiation  being  refracted  by  the  surface  before  propagating 
to  the  detector.  For  a  high  Index  surface,  all  refracted  beams  lie  very  near  the 
normal.  To  facilitate  viewing  the  information,  the  transmitted  portion  of  the 
graphs  are  plotted  as  a  function  of  the  beam  angle  incident  on  the  surface. 
This  is  equivalent  to  ««*H"g  the  reflecting  medium  a  thin  slab  of  material 
with  air  on  either  side,  and  patting  an  anti-reflection  coating  on  the  exiting 
surface.  The  beam  angle  Incident  on  the  front  surface  is  equal  to  the  beam 
angle  leaving  the  rear  surface.  Since  we  are  going  to  calculate  the  Mueller 
matrix  elements,  we  need  to  consider  the  TE  and  TM  polarizations  indepen¬ 
dently. 

Tn  calculate  the  Mueller  matrix,  we  first  calculate  the  scattering 
amplitude  matrix  which  Is  explained  by  Bohren  and  Huffman4  and  others.  M  For 
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our  parpoMi  we  define  the  scattering  amplitude  matrix  by  the  following  equa¬ 
tion: 

e£*\  /s.  s.\  /  ES] 

e;“J  =  vs*  s.jVEnj.  in 

where  the  vectors  represent  the  incident  and  scattered  electric  fields  for  the 
TE  and  TM  polarizations.  If  we  place  our  detector  in  the  plane  of  incidence 
(defined  by  the  incident  beam  and  the  normal  to  the  surface),  it  can  be  shown 
that  the  matrix  elements  St  and  S.  are  Identically  zero.  We  can  now-examine 
the  TE  mode  and  calculate  the  element  Si. 


THE  TE  MODE 


The  vibration  of  the  dipole  for  the  TE  mode  is  perpendicular  to  the 
plane  of  incidence.  As  a  result,  the  radiation  from  a  dipole  in  this  mode  has 
no  angular  dependence  in  this  plane  and  is  therefore  slightly  easier  to  under¬ 
stand  than  the  TM  mode.  Since  all  the  beams  have  the  same  temporal  depen¬ 
dence.  we  can  normalize  the  beam  from  path  I,  and  represent  all  other  beams 

in  terms  of  this  beam’s  amplitude  and  phase.  A  beam  following  path  II  will  be 
oat  of  phase  with  the  normalized  beam  when  it  reaches  the  dipole  by  an 

amount  5(d.,d),  where 


5(»,d)  *  -  25-2dcosd, 


(2) 


d  is  the  between  the  dipole  and  the  mirror,  and  *  is  the  wavelength 

of  the  light.  Furthermore,  the  amplitude  and  the  phase  are  affected  by  the 
Fresnel  reflectance  factor1,  RuUW,  where 


fcrfd) 


Zjcosd  -  Z.[I  -  (n</ ml* sln*S  ]** 
Zacosft  +  Z.[l  -  (ai/m^sin^]1* 


(3) 


,  JilSl  .  (4) 

Z»  tom 

The  subscripts,  1  and  2,  correspond  to  different  sides  of  the  surface.  In  our 
case,  the  plane  wave  is  traveling  from  material  1  and  striking  material  2.  For 
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most  cases,  the  light  initially  travels  through  air  or  vacuum  withoi  %  1  and  ti> 
*  n*.  Path  QI  is  similar  to  path  I,  bat  the  phase  difference  and  reflectance 
factor  arn  a  function  of  x-d.  Finally,  for  path  IV  we  have  to  consider  the 
phases  and  reflectance  factors  of  both  the  previous  paths.  The  net  result  of 
the  four  beams  at  the  detector  is  the  following: 

Slid)  =  1  *  Rr«(d.)e,8(*-d>  +  Rrt(x-d)e 18(11 'd  d>  +  R»i(d.)Rn(x-d)e‘(S(***d)*S(l'**d)> .  (5) 
where  cosd  <  0. 

We  now  consider  the  transmission  of  the  TE  mode.  Path^V,  like 
path  I  undergoes  no  phase  lags.  However,  one  must  consider  the  effect  of  the 
Fresnel  transmission  factor  on  the  scattered  beam  which  is  given  by 


Tti(d) 


_ 2  ZaCQSd _ 

Zx»sd  +  Zi[l  -  (ni/na)*sin*d]w* 


(6) 


The  beam  following  path  VI  undergoea  the  same  phase  lag  on  the  incident 
beam  as  that  of  path  II.  The  scattering  amplitude  matrix  element  for  the 
transmitted  beam  can  be  expressed  as  follows: 

S,(d)  »  [  1  Rn(d»)e,,<**d)  ]Ti«(d),  (7) 

where  coed  >  0.  With  the  scattering  amptitnde  matrix  element  for  the  TE 
mode  solved,  we  now  examine  the  TM  mode. 


THE  TM  MODE 


The  terms  used  to  calculate  the  scattering  matrix  for  the  TM  mode 
are  almost  Identical  to  those  of  the  TE  mode.  The  Fresnel  coefficients  of 
reflection  and  transmission  of  the  surface  are  slightly  different  and  may  be 
written  as  follows: 
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TWd) 


_ 2  Zaco  s& _ 

Zicosd  ♦  Zj[I  -  (nt/m^sin^]1'* 


(9) 


The  major  difference  Is  the  angular-dependent  dipole  scattering  that  occurs  for 
this  mode.  Bickel  and  Bailey'  provide  an  inciteful  discussion  of  the  reason  why 
the  scattering  of  this  mode  is  proportional  to  the  cosine  of  the  angle  between 
the  incident  and  the  scattered  beam.  As  the  detector  scans  in  the  plane  of  in¬ 
cidence,  the  amount  of  dipole  vibration  that  the  detector  sees  corresponds  to 
the  projection  of  the  dipole  in  the  detector  plane.  This  causes  the  cosine  de 
pendence.  With  a  bit  of  trigonometry,  we  can  derive  the  following  equations 
for  the  scattering  amplitude  matrix  element  S*<$) : 

Si(&)  3  cos(d-9J  *  cos(x-0-0.)Ri»i<«Je,,(a-d>  ♦  cos<*-a-0»>R™<x-0)e,*(1'~*d) 

+  cos(0-0*)Rn.(d»)Ri».(x-$)e,ma*d)*,<*_*,d>>,  for  coed  <  0  and.  (10) 


S»<0>  a  [  cos(0-0.)  +  cosiit Rm(&.) e1* ]  Tr*(3),  for  coed  >  0.  (11) 

With  both  polarizations  taken  into  account,  we  now  direct  our  attention  to  the 
scattering  Mueller  matrix  elements. 

IHE  MUELLER  MATRIX 

From  the  relations  given  by  Bohren  and  Huffman,  we  can  use  the  scattering 
amplitude  matrix  elements  to  calculate  the  Mueller  scattering  matrix  elements. 
In  the  plane  of  incidence,  only  four  of  the  sixteen  elements  of  the  Mueller 
matrix  are  of  Interest.  These  are 

Sm  *  |[lS4*  ♦  ISS*  ]. 

S*  ■  |  [  ISd*  -  ISil*  ]/s«. 

So  3  Re(SiS»*)/Sn,  and 


Su  3  Im(  S»Si*)/S«. 


(12) 
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When  we  consider  only  the  scatter  in  the  plane  of  incidence,  the  other  scat¬ 
tering  element*  are  either  zero  or  simple  multiples  of  these  elements. 

Using  these  expressions,  the  four  Mueller  matrix  elements  were  cal¬ 
culated  for  a  dipole  with  no  surface.  The  surface  can  be  removed  by  setting 
the  surface  index  equal  to  i.  These  elements  are  shown  in  Figure  4.  It  is 
noted  that  the  scattering  angle  8  goes  from  0*  to  360*.  As  expected,  symme¬ 
tries  occur  about  6  *  180*.  This  is  a  fundamental  case  discussed  in  many  of 
the  references,  especially  Bickel  and  Bailey.  These  scattering  elements  are 
useful  references  and  will  be  compared  with  those  generated  when  a  surface  is 
placed  near  the  dipole. 

SCATTERING  AS  A  FUNCTION  OF  PARTICLE- SURFACE  SEPARATION 

Figure  S  shows  the  four  scattering  Moeller  matrix  elements  as  a 
function  of  dipole- surface  separation  for  a  beam  incident  at  normal  incidence 
and  a  surface  of  index  n  *  10.0  -  0.01.  Such  a  surface  coold  be  made  with  a 
series  of  1/4-wave  dielectric  layers  on  a  glass  snbstrate.  It  is  more  instinctive 
to  look  at  the  effects  of  a  high  index  dielectric  surface  intitially  (rather  than 
a  simple  glass  plate  or  a  metallic  surface)  because  a  high  index  interface  gives 
interfering  beams  that  are  of  nearly  equal  amplitude.  Also,  a  metal  interface 
creates  phase  shifts  upon  reflection.  Before  examining  Figure  S  in  detail,  it  is 
instinctive  to  examine  the  scattering  amplitude  matrix  elements.  To  provide 
some  Insight  for  the  TE  case,  we  rewrite  the  scattering  amplitude  matrix  ele¬ 
ment  (equation  S)  la  the  following  form: 

S«W  «  [l  ♦  Rn(0Je“w*  d>][l  +  Rt«(*-8>e,M*"*,d)].  (13) 

The  lint  term  which  is  independent  of  8  controls  the  amplitude  of  the  overall 
pattern.  As  the  dipole-surface  separation  is  Increased,  a  stationary  detector  at 
8  will  witness  the  intensity  element  Sn(8)  fluctuating  sinusoidally  with  distance 
due  to  the  first  term.  For  a  beam  incident  along  a  normal  to  a  dielectric 
surface  illuminating  a  particle  next  to  the  surface  (d  *  0),  the  intensity  mea- 
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sured  by  a  fixed  detector  will  be  in  a  trough  (minima)  since  there  is  a  phase 
change  at  the  surface.  As  the  particle  moves  away  from  the  surface,  the  signal 
increases  until  it  is  1/4  wave  from  the  surface,  at  which  point  the  signal  has 
cycled  into  a  peak  (maxima)  since  S(8o,d)  =  -it.  At  1/2  wave  from  the  surface, 
the  signal  has  returned  to  the  trough.  The  intensity  difference  between  a  peak 
and  a  trough  can  be  calculated  using 

M  =  log{|[l  ♦  R(0)] ]*}  -  log{j[l  -  R(0)]|’}  .  114) 

For  n  =  10.0  -  0.01,  this  yields  il  >  2  or  a  two  order  of  magnitude  change  in 
the  scattering  by  moving  the  particle  1/4  wave.  This  is  most  easily  seen  in  the 
transmission  region  (0  3  0*)  of  the  S«  curves  of  figure  S. 

The  second  term  in  the  S«(8)  equation,  controls  tig,  shape  of  the 
pattern.  If  the  first  order  approximation  assumes  R(8)  to  be  constant,  then  the 
pattern  behaves  sinsusoldaily  with  5(8,d).  At  0  *  180*,  we  not  only  see  the 
pattern  fluctuating  two  orders  of  magnitude  with  the  amplitude  envelope,  but 
also  varying  another  two  orders  of  magnitude  due  to  the  scattering  (second) 
term  of  equation  13,  for  a  total  fluctuation  of  four  orders  of  magnitude.  For 
other  scattering  angles,  the  fluctuation  of  the  amplitude  envelope  and  the 
fluctuation  at  a  particular  scattering  angle  will  be  out  of  phase  and  will  tend 
to  cancel  each  other.  For  a  dipole-surface  separation  of  3X/4  and  the  detector 
placed  at  8  *  180*.  5(8, d)  >  3k.  At  this  location,  there  is  a  peak  in  the  scat¬ 
tering  lnte natty.  A  trough  will  appear  at  any  scattering  angle  where  the  phase. 
S(8,d>  *  2mx.  where  j£  is  a  positive  Integer  not  greater  than  5(180*. d)/2x  (3/2 
in  this  easel.  Keeping  d  constant  with  m  *  1,  we  can  solve  our  phase  equation. 
5(8, d)  *  2u,  for  scattering  angle  8,  yielding  132*  and  228*.  At  about  these  lo¬ 
cations  in  Figure  S  we  see  troughs  whose  value  haan't  changed  much  from  the 
amplitude  envelope  of  the  half  wave  separation.  These  nodes  move  in  to  the 
center  of  the  pattern  (toward  180*)  as  the  dipole-surface  distance  increases. 
Another  set  of  troughs  (for  n  *  0)  la  Just  beginning  to  form  for  5(8, d)  3  0 
which  corresponds  to  8  >  90*  and  270*.  but  they  are  not  as  apparent  at  this 


point  because  the  Fresnel  reflection  term  dominates. 
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The  response  of  the  polarization  matrix  elements  (S*.  S»,  and  Sj.)  is 
not  as  dramatic  as  the  total  intensity  matrix  element,  S« .  Spikes  tend  to 
appear  on  these  elements  where  there  are  troughs  in  the  total  intensity  ele¬ 
ment.  These  spikes  seem  similar  to  those  that  appear  in  the  Rayletgh-Gans  li¬ 
mit  of  single  sphere  scattering.'  S«  has  an  interesting  feature:  for  very  small 
dipole- surface  separations  (d  ~  X/1000),  the  polarization  has  switched  sign. 
This  appears  to  be  a  function  of  the  refractive  index  of  the  surface  \shich  we 
now  consider. 

Figure  6  shows  the  S*  scattering  as  a  function  of  dipole  separation 

from  a  surface  of  index  n  *  28.0  -  95.01,  an  index  roughly  corresponding  to 

aluminum  at  10(un  wavelength.  First  note  that  the  transmission  portions  of  the 
polarization  elements  carry  very  little  information:  i.e.,  they  are  either  zero  or 
I00X  over  the  entire  range.  This  is  not  too  disturbing  since  the  skin  depth  of 
this  metal  would  effectively  prevent  most  transmission  measurements.  This 
information  loss  also  occurs  for  high  index  dielectric  surfaces,  but  is  more 
apparent  In  figure  6  because  the  reflectance  of  aluminum  Is  greater  than  the 
dielectric  of  figure  5.  An  interesting  feature  has  cropped  up  (n  matrix  element 
Sn(&)  at  normal  Incidence  when  d  *  X/2.  At  8  *  180*;  where  a  node  occurred 
for  the  dielectric  surface  of  figure  S,  a  small  maximum  exists.  Reflections 

from  a  metallic  surface  change  the  state  of  the  incident  beam  from  linear 

polarization  to  elliptical  polarization.  As  a  result  of  these  reflections,  the 
shape  and  position  of  the  nodes  are  different  than  the  nodes  for  a  dielectric 
surface.  The  nodes  and  spikes  for  aluminum  are  sharper  than  for  the 
dielectric,  huh  this  is  s  function  of  the  reflectance*,  the  aluminum  la  more 
reflecting,  so  the  spikes  are  sharper.  Other  than  the  small  displacements  of 
the  maxima/minima  and  the  loaa  of  information  in  matrix  element  Sm  for  large 
dipole-surface  separations,  it  is  hard  to  distinguish  the  scattering  from  the 
metai  surface  of  figure  6  from  the  high-index  dielectric  surface  of  figure  5. 
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Figure  7  shows  the  &*  scattering  as  a  function  of  dipole  distance 
from  a  surface  of  Index  n  *  1  S  -  0.01,  which  roughly  corresponds  with  glass 
In  the  visible  spectrum.  For  a  tow- index  surface,  the  matrix  element  curves  are 
much  smoother,  and  very  closely  resemble  the  curves  for  a  dipole  without  a 
surface  shown  in  Figure  4.  As  the  dipole-surface  distance  is  increased,  the 
elements  change,  while  retaining  this  general  shape. 

SCATTERING  AS  A  FUNCTION  OF  INCIDENT  ANGLE 

We  now  consider  the  final  remaining  parameter  for  dipole-surface 
scattering,  the  incident  angle  of  the  Incoming  radiation.  We  examine  the  spe¬ 
cific.  but  arbitrary,  case  for  the  particle  near  the  surface  at  d  3  X/10.  This 
case  could  represent  a  2nm  diameter  dust  particle  on  a  surface  Illuminated 
with  10.6  um  light.  Although  a  radius  of  X/10  borders  on  the  Rayleigh  limit,  it 
has  practical  applications. 

Figure  8  shows  the  S«  scattering  as  a  function  of  Incident  illumina¬ 
tion  angle  for  a  dipole  located  X/10  from  a  surface  of  index  n  *  10.0  -  0.01. 
First  we  note  the  lack  of  symmetry  in  the  elements.  Prior  to  this  study,  all 
the  elements  have  been  symmetric  about  8  *  180*.  As  8.  Is  increased  first  to 
30*  and  then  to  60*.  the  elements  vary  significantly  from  those  of  the  isolated 
dipole  of  figure  4.  Specular  peaks  occur  in  those  directions  where  the  beam 
takes  a  direct  path  and  does  not  Interact  with  the  particle.  These  occur  at  8  = 
8*  on  the  tranaimissioa  side  of  the  surface  and  8  *  x  -  8.  on  the  reflection 
side  of  the  surface.  It  may  seem  surprising  to  And  that  the  peak  of  the  scat¬ 
tered  radiation  is  not  necessarily  at  the  specular  beam  on  the  reflecting  side 
of  the  surface. 

Established  patterns  in  the  polarization  matrix  elements  for  scatter¬ 
ing  from  a  single.  Isolated  sphere  no  longer  hold  when  the  incident  beam  is 
no  longer  normal  to  the  surface.  Element  S»»  which,  for  an  isolated  sphere.  Is 
equal  to  zero  both  at  specular  and  at  180*  from  specular,  no  longer  obeys  this 
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rate.  For  a  high-index  surface  the  signal  doea  tend  to  -tOOX  for  scattering 
along  the  surface  ($3  *90*).  Element  S»,  which  for  an  isolated  sphere  is  equal 
to  100X  at  specnlar  and  -100X  in  the  backscatter  (180*  from  specular)  no  lon¬ 
ger  obeys  these  rules  either.  For  a  highly  reflecting  surface  it  does  tend  to 
zero  along  the  surface  (8=  *90*).  Element  Si*,  which  for  a  lone  sphere  is  equal 
to  zero  at  specular  and  180*  from  specular,  now  tends  to  zero  only  along  the 
surface  (8=  *90*).  These  tendencies  can  be  derived  from  the  scattering  ampli¬ 
tude  matrix  elements.  When  8  =  *90*.  R(8)  =  -1  and  $(8,d)  =  0,  the  elements 
reduce  to  zero.  For  surfaces  of  high  index  of  refraction,  the  Fresnel  reflection 
coefficient  Rn(8)  approaches  -I  more  rapidly  than  Rn»<8),  which  goes  to  zero  at 
Brewster's  angle.  As  a  result  of  this  dissimilarity,  the  beam  is  almost  entirely 
TM  polarized,  and  the  Mueller  matrix  element  S«  will  go  to  IOOX.  The  other 
polarization  elements  must  go  to  zero  at  this  point.  Since  Si  and  St  go  to  zero 
along  the  surface,  element  S«  goes  to  zero  along  the  surface. 

Figure  9  shows  the  scattering  as  a  function  of  incident  angle  for  a 
dipole  X/10  from  a  surface  of  index  n  *  28.0  -  95.01.  These  curves  are  similar 

to  the  curves  of  Figure  8  where  n  *  10.0  -  0.01.  The  only  differences  are  the 

magnitudes  in  some  locations  and  the  loss  of  information  in  the  reflected  por¬ 
tion  of  the  Sm  matrix  element.  The  change  to  a  higher,  complex  index  has  had 
very  Itttle  effect  on  the  elements. 

Figure  K)  shows  the  scattering  as  a  function  of  incident  angle  for  a 

dipole  X/10  from  a  surface  of  index  a  *  1.5  -  0.01.  The  matrix  elements  for  a 

low-index  surface  are  expected  to  retain  the  shape  of  the  elements  for  the 
isolated  dipole  shown  In  Figure  4.  The  peaks  in  the  elements  are  very  near  the 
same  angular  locations  that  they  would  be  with  no  surface  present.  The  maxi- 
ma/mMma  of  matrix  elements  S«  and  Sm  reveal  the  incident  angle  of  the 
Incident  radiation.  The  trough  located  90*  from  specular  in  the  S«  curve  is  still 
close  to  90*  from  specular.  The  maxima/ minima  of  the  Sm  curves  also  become 
shifted  toward  the  specular  and  backscatter  locations  as  the  incident  angle  is 
changed. 
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roMfLUSlON 

Determining  the  nature  of  systems  by  studying  the  light  scattered 
from  them  demands  an  understanding  of  how  basic  properties  of  the  systems 
affect  the  scattering.  Reducing  the  size  of  the  particle  to  the  Rayleigh  limit 
eliminates  the  effects  of  the  particle's  geometry  and  refractive  index  from  the 
problem.  This  effectively  reduces  the  number  of  parameters  to  three:  the 
complex  index  of  refraction  of  the  surface,  the  particle-surface  separation,  and 
the  incident  illumination  angle.  Understanding  the  trends  in  a  simply  system 
consisting  of  only  three  parameters  leads  to  a  better  understanding  of  what 
happens  in  more  complex  systems.  If  we  are  unable  to  comprehend  how  a 
simple  system  scatters  light,  we  cannot  hope  to  understand  the  scattering 
from  a  system  that  is  more  complex. 

Finally,  it  may  be  asked  If  any  of  our  conclusions  are  valid  for  a 
particle  of  finite  size.  Figure  II  compares  two  sets  of  curves.  One  Is  the  set 
of  Mueller  matrix  elements  for  a  dipole  separated  d  *  X/10  from  an  aluminum 
mirror  illuminated  at  an  angle  of  incidence  of  30*.  The  other  is  a  set  of  theo¬ 
retical  Mueller  matrix  elements  for  the  same  system  with  the  dipole  replaced 
by  a  X/10  radius  sphere  of  Index  i.5  resting  on  the  surface.  The  theoretical 
system  has  been  derived  by  solving  the  boundary  conditions  at  both  the  sphere 
and  the  surface.  Its  solution  is  the  subject  of  another  paper*  The  two  curves 
are  quite  similar,  with  differences  arising  from  the  geometry  and  refractive 
index  of  the  finite-size  sphere.  The  scattering  from  such  a  system  can  be 
measured.  Weber  end  Hirleman  made  scattering  measurements  from  a  single 
sphere  of  radius  ~X/2  using  a  He-Ne  laser  for  illumination*  A  higher  power 
laser  in  the  Infrared  region  would  facilitate  the  measurements.  On  a  more 
practice!  note,  this  model  could  be  used  to  measure  sub-wavelength  contami¬ 
nant  sizes  and/or  densities  by  piecing  detectors  at  appropriate  scattering 
locations  and  monitoring  signal  changes  in  the  matrix  elements  as  the  surface 
is  scanned. 
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Figure  Caption* 

Fig.  1.  Th#  geometry  of  a  dipole- surface  scattering  system:  the  dipole,  repre¬ 
sented  by  a  circle,  rest*  a  distance,  d.  from  a  surface  of  index  n.  The 
incident  angle  0.  and  the  scattering  angle  $  are  measured  from  the  sur¬ 
face  normal. 

Fig.  2.  The  path*  that  an  Incident  beam  may  follow  to  the  dipole  and  after 
interaction  before  being  detected  on  the  incident  side  of  the  surface. 

Fig.  3.  The  paths  that  a  beam  may  follow  before  being  detected  on  the  tran¬ 
smitted  side  of  the  surface. 

Fig.  4.  Th*  Mueller  matrix  element*  for  an  isolated  dipole. 

Fig.  5.  The  Moeller  matrix  elements  as  a  function  of  dipole  distance  from  a 

surface  of  index  n  *  10.0  -  0.01  illuminated  at  normal  incidence.  » - • 

Distance  *  X/1000.  0—0  Distance  *  X/4.  X - X  Distance  =  X/2. 

+ - ►  Distance  »  3X/4. 

Fig.  6.  The  Moeller  matrix  elements  as  a  function  of  dipole  distance  from  a 

surface  of  index  n  *  28.0  -  95.01  illuminated  at  normal  incidence.  « - • 

Dietance  *  X/1000.  O— — O  Distance  *  X/4.  X - X  Distance  *  X/2. 

» . »  Distance  *  3X/4. 

Fig.  7.  The  Mueller  matrix  elements  as  a  function  of  dipole  distance  from  a 

surface  of  Index  n  »  li  -  0.01  illuminated  at  normal  incidence.  • - • 

Distance  •  X/1000.  O - O  Distance  ■  X/4.  X - X  Distance  *  X/2. 

- - ►  Distance  *  3X/4. 

Fig.  8.  The  Moeller  matrix  elements  as  a  function  of  incident  angle  of  illumi¬ 
nation  for  a  dipole  a  distance  X/10  from  a  surface  of  index  n  *  10.0  - 
0.01.  « - •  0.  *  0*.  0—0  0.  *  30*.  X - X  ■  60*. 
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Fig.  9.  The  Mneller  matrix  elements  as  a  function  of  incident  angle  of  Illumi¬ 
nation  for  a  dipole  a  distance  X/10  from  a  surface  of  index  n  =  28.0  - 
95.01.  • - •  8.  s  0*.  O 0  8.  =  30*.  X - X  8.  *  60*. 

Fig.  10.  The  Mueller  matrix  elements  as  a  function  of  incident  angle  of  illumi¬ 

nation  for  a  dipole  a  distance  X/10  from  a  surface  of  index  n  =  1.S  - 
0.01.  • - •  8.  *  0* .  O - O  8.  =  30*.  X - X  8.  =  60*  . 

Fig.  11.  The  Mueller  matrix  elements  as  a  function  of  incident  angle  of  illumi¬ 
nation  for  a  particle  X/10  from  a  surface  of  index  n  *  28.0  -  95.01  and 

Incident  angle  of  30*.  • - •  A  dipole.  O—O  A  sphere  of  index  n  = 

1.55  -  0.01  and  radios  X/10. 
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Abstract 

Although  the  general  case  of  electromagnetic  waves  scattered  from 
a  sphere  when  illuminated  by  a  plane  wave  was  formulated  by  Mie  more  than 
half  a  century  ago,  almost  all  measurements  and  calculations  of  the  total 
scatter  have  concentrated  on  the  far-field  region.  We  examine  the  electric  and 
magnetic  field  topographies  as  well  as  the  scattering  Mueller  matrices  in  the 
near-field  region  of  a  Mie  sphere. 
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L  Introduction 

In  the  nineteenth  and  early  twentieth  centuries,  one  of  the  most 
important  problems  in  the  application  of  Maxwell's  theory  to  phenomenon  tak¬ 
ing  place  in  nature  was  the  interaction  of  electromagnetic  waves  with  matter. 
An  important  example  is  the  work  of  Rayleigh'  in  explaining  the  blue  color  of 
the  sky  in  terms  of  scattering  from  particles  small  in  comparison  with  the 
optical  wavelengths.  An  analysis  of  electromagnetic  waves  scattered  from  a 
sphere  was  published  by  Clebsch  in  1863,’  and  like  Fresnel's  reflection  equa¬ 
tions,  his  solutions  appeared  before  the  development  of  Maxwell's  equations. 
Although  the  radiation  scattered  by  a  sphere  was  solved  by  Mie  in  1908,*  and 
refined  by  Debye  in  1909 4  when  he  investigated  the  incident  electromagnetic 
wave  pressure  on  particles,  the  result  was  also  obtained  by  Lorenz  in  1898.* 
The  geometry  of  spherical  systems  is  especially  appealing  because  spheres 
possess  high  symmetry  and  they  occur  naturally,  being  a  condition  of  minimum 
potential  energy.  Extensive  applications  to  atmospheric  sciences  occur  where 
rainbows,  planetary  atmospheres  and  other  systems  are  explored  in  books  by 
van  de  Hulst,*  Kerker/  and  Bohren  and  Huffman.' 

The  relations  of  Mie's  general  theory  are  relatively  simple  when  the 
distance  kR  to  the  sphere  lay  in  the  electromagnetic  far-field  region  where 
only  scattered  fields  decaying  as  1/kR  are  significant.  Fifty  years  after  Mie 
presented  his  theory,  King  and  Wu*  lamented  the  fact  that  "sufficient  data  do 
not  seem  to  be  available  to  construct  a  complete  representation  about  a  con¬ 
ducting  sphere"  as  that  shown  for  an  infinitely-long,  conducting  cylinder  (in 
referring  to  figures  in  the  text  showing  the  near- field  region  of  the  cylinder). 
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In  the  seme  text,  they  devoted  an  appendix  to  plots  of  surface  currents  on 
perfectly-conducting  spheres  of  increasing  size  illuminated  by  incident  plane 
waves  which  were  obtained  from  a  large,  state-of-the-art  military  computer 
program  directed  by  Nelson  Logan  for  USAF  Cambridge  Research  Center. 

Most  theoretical  work  and  scattering  applications  consider  only  the 
far  field,  optically  represented  as  the  Fraunhofer  region.  For  instance,  all  three 
scattering  programs  in  the  book  by  Bo  hem  and  Huffman  (ref.  8)  calculate  the 
far  fields.  One  notable  exception  is  the  work  of  Aden9  who  ingeneously  mea¬ 
sured  the  near-zone  fields  of  metal  and  water  spheres  along  the  back-scatter 
direction  in  the  microwave  region  of  the  spectrum.  With  the  advent  of  better 
computers,  the  entire  electromagnetic  fields  can  now  be  calculated  easily.  It  is 
instructive  to  see  what  happens  as  the  detectors  move  into  the  near- field 
region,  optically  represented  as  the  Fresnel  region.  We  want  to  know  the  error 
involved  in  assuming  a  far-field  solution  when  we  know  our  detector  is  not  at 
infinity.  We  also  want  to  know  just  how  the  near-field  region  differs  from  the 
far-field  region. 

In  experiments  where  scattering  systems  are  illuminated  with  visible 
wavelengths  (X  *  0.6328pm,  for  instance),  a  detector  mechanism  only  one 
„  meter  distant,  would  certainly  seem  to  be  operating  in  the  far-field  region. 
At  this  distance,  changing  the  scatterer- detector  distance  by  a  factor  of  two 
would  not  change  the  relative  shape  of  a  scattering  signal,  therefore  the  de- 
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tector  is  considered  to  be  in  the  far  field  for  all  practical  considerations. 
Indeed,  Bell's'1  measurements  of  the  light-scattering  Mueller  matrix  for  a 
fiber,  seemed  to  verify  this  (although  his  optical  system  was  typical,  it  cer¬ 
tainly  was  not  universal). 

In  this  paper  we  explore  the  near-field  region  of  a  Mie  scattered 
First  we  will  calculate  the  entire  field.  We  will  then  calculate  the  scatter  that 
would  be  measured  in  the  far-field  region,  and  compare  the  results  with  what 
would  be  measured  in  the  near-field  region. 

!L  Mie  Theory 

Figure  1  shows  the  arrangement  of  the  scattering  system.  A  sphere 
of  radius  a,  and  complex  refractive  index  n'  is  centered  on  the  origin.  A  plane 
wave,  electrically  polarized  in  the  x  direction,  propagates  in  the  positive  z 
direction.  We  detect  the  electric  and  magnetic  fields  at  position  R.  $,  <p  from 
the  center  of  the  sphere  (R  is  the  sphere-detector  distance). 

We  report  only  the  primary  results  of  Mie  theory  since  formal  deri¬ 
vations  exist  in  many  of  the  cited  texts.  The  electric  and  magnetic  components 
of  the  fields  incident  upon  the  sphere  are 
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where  k  *  2 it/\  is  the  wavenumber,  o>  is  the  angular  frequency,  and  p  is  the 
permeability  of  the  incident  medium  at  X.  The  vectors  1 U'1  and  N—1'1  are  given 
by 
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The  superscript  identifies  the  Bessel  functions  used  to  represent  the  fields. 
For  example,  j  *  1  corresponds  with  the  use  of  the  functions  z.‘"(p)  *  j.(p),  j  = 
2  corresponds  with  the  use  of  the  functions  z»'”(p)  *  y.(p),  and  j  *  3  corre¬ 
sponds  with  the  use  of  the  functions  z.lll(p>  *  h."’(p)  »  j.(p)  ♦  ly.(p).  The  func¬ 
tions  Pr(cosd)  are  the  normalized  associated  Legendre  polynomials  defined  by 
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For  an  incident  electromagnetic  plane  wave  traveling  along  the  positive  z  axis 
and  polarized  along  the  x  axis,  the  coefficients  are 
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where  a  pj-(p),  and  5»(p)  =ph."'(p).  These  coefficients  can  be  reduced  to 
much  simpler  expressions  when  the  refractive  index  of  the  sphere  approaches 
large  values,  in  which  case 
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We  will  detect  the  fields  in  two  different  planes,  one  in  the  y-z 
plane  (9  »  90*)  which  lies  perpendicular  to  the  incident  electric  field  (the  TE 
plane,  or  H-plane),  and  the  other  in  the  x-z  plane  (9  *  0*)  which  lies  per¬ 
pendicular  to  the  incident  magnetic  field  (the  TM  plane,  or  E-plane).  We  can 
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represent  the  fields  in  the  form  of  matrices: 
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S*  =  ^2n(n^l)^~^  P.'(cosd)f. 
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These  expressions  hold  for  the  fields  everywhere.  To  get  to  the  far-field  re¬ 
gion  we  let  kR  •*  co,  where  these  expressions  are  further  reduced  by  taking 
advantage  of  the  asymptotic  behavior  of  the  spherical  Hankel  function: 
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»P 


kR  I-  sind 


kR  *-  sind 


9 
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As  the  sphere-detector  distance  is  increased,  matrix  elements  S»  and  Si’  (which 
vary  as  1/R*)  approach  zero  more  rapidly  than  do  the  other  matrix  elements 
(which  vary  as  1/R).  We  sec  here  another  significance  of  the  far  field.  The 
E  and  H  field  vectors  are  not  only  orthogonal  to  one  another  and  to  propaga¬ 
tion  direction  I,  but  are  also  in  time  and  spatial  phase  so  as  to  compose  a 
wave  front  moving  at  the  speed  of  light  in  the  wave  medium  (equations  12c,d). 

1IL  Hfi2L  fit  il  IicZ 

Equations  9-12  can  be  used  to  determine  whether  the  radiation  from 
a  scattering  system  is  in  the  near  or  far  field.  As  an  example  we  start  by 
examining  the  scattered  fields  from  a  dielectric  (quartz)  sphere  of  refractive 
index  n'  *  1.S5  ♦  O.Oi  in  the  optical  region  of  the  spectrum  (X  •  0.6328pm 
HeNe  laser  radiation).  We  arbitrarily  choose  the  radius  of  the  sphere  to  be 
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a  2  O.S25(im,  a  size  large  enough  to  provide  some  oscillatory  structure  on  the 
scattering  curves. 

A*  The  Electric  Fields 

Figure  2  shows  the  intensity  of  the  electric  field*  as  a  function  of 
scattering  angle  at  four  different  observation  distances.  The  TE  curves  corre¬ 
spond  to  measurements  when  the  incident  electric  field  is  transverse  to  the 
scattering  plan*  (a  measurement  of  Si"Si).  The  TM  curves  correspond  to  mea¬ 
surements  of  the  non-radial  component  of  the  electric  field  when  the  incident 
magnetic  field  is  transverse  to  the  scattering  plan*  (a  measurement  of  &*&). 
The  RD  curves  correspond  to  measurements  of  the  radial  component  of  the 
electric  field  when  the  magnetic  field  is  transverse  to  the  scattering  plane  (a 
measure  of  S*a&).  Four  sets  of  curves  are  shown  respectively  for  sphere-de¬ 
tector  distances  of  R  »  IX,  R  «  3X,  R  *  9X,  and  R  *  oo.  To  facilitate  compari¬ 
son,  the  elements  have  all  been  normalized  for  sphere-detector  distance  [mul¬ 
tiplied  by  (IcR)1  for  TE  and  TM  and  (kR)‘  for  RD],  Not*  that  at  the  endpoints 
(8  *  O’  and  8  »  180*)  the  TE  elements  are  equal  to  the  TM  elements  and  the 
RD  elements  are  zero. 

These  curves  show  that  at  optical  wavelengths,  the  far  field  is  not 
very  far  away  from  the  scatterers.  This  is  shown  by  the  fact  that  three  curves 
(R  *  3X,  R  »  9X,  and  R  >  oo)  in  each  set  of  four  have  very  similar  shape  and 
almost  coincide,  while  the  fourth  curve  (R  *  IX)  is  quite  different.  When  the 
sphere-detector  distance  is  only  R  ■  3X,  the  matrix  elements  are  already  simi- 
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lar  to  those  for  R  >  oo.  When  R  =  9X  (less  than  6pm  at  optical  wavelengths!) 
the  elements  are  almost  indistinguishable  from  far-field  elements.  This  shows 
that  near-field  phenomena  occurs  only  very  close  to  the  particle,  and  that  it  is 
very  difficult  to  measure  near-field  signals  in  scattering  experiments  at  optical 
wavelengths.  An  infinitely-small,  non-interfering  detector  would  have  to  probe 
the  field  at  R  ~  X. 

B.  The  Magnetic  Fields 

We  next  look  at  the  magnetic  fields  and  compare  them  with  the 
electric  fields  for  the  same  system.  Figure  3  shows  the  intensities  of  the  mag¬ 
netic  fields  as  a  function  of  scattering  angle  at  four  different  observation  dis¬ 
tances  (note  that  TE  corresponds  to  S*mS»,  TM  to  S»*S»\  and  RD  to  S»'"S$'). 
Four  sets  of  curves  are  shown  respectively  for  the  same  sphere-detector  dis¬ 
tances  of  R  *  IX,  R  »  3X,  R  ■  9X,  and  R  *  ®  used  to  calculated  the  electric 
field  intensities  of  Fig.  2.  To  facilitate  comparison,  the  elements  have  all  been 
normalized  for  sphere-detector  distance  and  for  field  strength  (multiplied  by 
up/k).  As  with  the  electric  field  curves,  three  curves  (R  *  3X,  R  *  9X,  and  R  = 
oo)  in  each  set  of  four  have  very  similar  shape  and  almost  coincide,  while  the 
fourth  curve  (R  *  IX)  is  quite  different.  Only  in  the  far  field  when  R  1  ®  does 
the  magnetic  field  behavior  approach  that  of  the  electric  field,  the  ratio  bet¬ 
ween  them  being  up/k.  The  TE  and  TM  electric  and  magnetic  field  intensity 
curves  are  identical  in  the  far  field.  Since  the  radial  component  approaches 
zero  much  more  rapidly  (as  1/R4)  than  the  other  components  (which  decay  as 
1/R*),  it  does  not  contribute  to  the  intensity  in  the  far  fields. 


c.  Jh£  Matrix  Elements 
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Finally,  we  generate  the  Mueller  matrix  elements'3  from  the  above 
expressions  for  the  fields  and  examine  their  behavior  as  a  function  of  sphere- 
detector  distance  R.  The  Mueller  matrix,  which  is  defined  in  terms  of  the 
electric  field  components  Si  and  Sa,  is  especially  attractive  experimentally 
because  it  completely  quantifies  the  polarization  state  of  the  system.  The 
matrix  does  not,  however,  include  the  radial  component  of  the  system.  This 
component  which  decays  as  1/R1  contributes  only  to  the  total  field  and  not  to 
radiation  from  the  system  since  the  radiation  term  decays  as  1/R.  An  ideal 
(infinitely-small,  perfectly-impedence-matched)  detector  placed  in  the  near¬ 
field  region  cannot  distinguish  between  the  radiation  terms  and  the  total 
electric  field.  Only  four  elements  of  the  16-element  matrix  are  non-zero  and 
unique.  Only  these  are  necessary  to  characterize  a  symmetric  scatterer  like  a 
sphere.  These  four  Mueller  matrix  elements  for  the  a  *  0.525pm  radius  quartz 
sphere  calculated  from  the  data  shown  in  Figs.  2  and  3  are  presented  in  Figure 
4. 


The  Mueller  matrix  is  an  especially  useful  diagnostic  tool  for  expe¬ 
rimentalists  because  it  is  self-calibrating.  The  polarizations  are  independent  of 
-*  intensity.  It  is  important  to  emphasize  the  value  of  the  matrix  elements  as 
discriminators  of  small  differences.  The  electric  fields  measured  at  R  ■  IX  in 
Fig.  2  can  easily  be  distinguished  from  the  others  measured  at  R  ■  3X,  R  «  9X, 
and  R>  a.  However  they  are  proportional  to  the  source  intensity  and  detector 
efficiency.  If  these  parameters  are  unknown,  the  absolute  intensities  are  un- 
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known  and  the  similarity  in  shapes  will  prevent  discrimination  of  these  four 
curves.  The  situation  is  entirely  different  with  the  Mueller  matrix  elements  of 
Fig.  4.  Although  the  curves  are  similar,  they  are  distinguishable.  A  detector 
measuring  element  Sh  near  13S*  can  easily  distinguish  all  the  elements  at  the 
four  different  sphere-detector  distances,  and  the  signal  is  completely  indepen¬ 
dent  of  the  source  intensity. 

#■ 

IV.  What  about  conductivity? 

Earlier  in  the  paper  we  simplified  the  equations  representing  the 
scattering  amplitude  coefficients  by  letting  the  sphere  refractive  index 
approach  infinity.  Making  the  sphere  a  perfect  conductor  provides  an  interest¬ 
ing  example  because  it  provides  one  limit  to  which  the  refractive  index  of  the 
sphere  may  approach.  Another  interesting  example,  known  as  the  Rayleigh-Gans 
region,  occurs  when  the  refractive  index  of  the  sphere  approaches  the  refrac¬ 
tive  index  of  the  encompassing  medium.  In  this  limit,  where  the  refractive 
index  of  the  sphere  is  identical  to  the  refractive  index  of  the  encompassing 
medium,  the  sphere  disappears.  It  is  optically  invisible  and  no  scatter  can 
occur.  The  intensities  of  the  electric  fields  as  a  function  of  scattering  angle  at 
different  sphere-detector  distances  for  a  perfectly-conducting,  a  3  0.525pm  ra¬ 
dius  sphere  Illuminated  at  X  »  0.6328pm  are  shown  in  Figure  5  and  the  corre¬ 
sponding  intensities  of  the  magnetic  fields  are  shown  in  Figure  6.  The  curves 
for  the  conducting  sphere  are  much  smoother  than  those  for  the  dielectric 
sphere  (shown  in  figures  2  and  3)  because  the  fields  are  unable  to  penetrate 
the  conducting  sphere.  The  scatter  from  a  perfectly-conducting  sphere  and  a 
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dielectric  sphere  can  be  compared  to  the  reflection  of  a  light  source  from  a 
metal  mirror  (with  one  boundary)  and  a  sheet  of  glass  having  two  boundaries. 
The  metal  mirror  provides  a  well-defined  reflection,  but  the  intensity  of  the 
light  reflecting  from  the  sheet  of  glass  will  oscillate  through  maxima/minima 
due  to  interference  occuring  between  the  fields  which  reflect  off  the  front  and 
rear  boundaries  of  the  glass.  These  oscillations  are  perhaps  more  pronounced 
in  the  Mueller  scattering  matrix  shown  in  Fig.  7. 

Finally  we  should  see  if  the  approximation  of  large  refractive  index 
has  any  physical  significance.  To  do  this  we  examine  the  scatter  from  an  a  » 
O.S2Spm  radius  partially-conducting  sphere  of  aluminum  (n  3  1.1  ♦  6.0i)  illumi¬ 
nated  at  X  3  0.6328pm  as  a  function  of  scattering  angle  at  four  different 
sphere-detector  distances  R.  The  matrix  elements  for  this  system  are  shown  in 
Figure  8.  They  are  very  similar  to  the  elements  of  Figure  7  for  perfectly  con¬ 
ducting  spheres.  Therefore,  even  though  the  refractive  index  of  aluminum  is 
finite,  its  scattering  properties  can  be  estimated  by  considering  scattering  in 
the  limiting  case  of  a  perfect  conductor. 

V.  The  Total  Field 

The  near  field  can  be  distinguished  from  the  far  field  by  placing  a 
“perfect'*  detector  within  a  few  wavelengths  of  a  sphere  whose  size  is  compar¬ 
able  with  the  wavelength  of  the  illuminating  radiation.  Unfortunately  the  de¬ 
tector  will  measure  a  mixture  of  the  scattered  fields  and  the  incident  fields 
since  there  is  no  way  to  filter  out  the  incident  fields  from  the  scattered 


fields.  Therefore  the  incident  fields  must  he  included  in  measurements  taken 
in  the  near  field. 

This  does  not  occur  experimentally  for  most  scattering  measure¬ 
ments  since  they  are  made  in  the  far-field  region.  The  laser  illumination, 
though  not  an  infinite  plane  wave,  is  very  nearly  constant  over  the  scatter- 
ing  system  (sphere).  Also,  because  the  laser  beam  has  a  small  diameter  (~ 
1000a),  the  angular  subtense  of  the  incident  field  is  very  narrow,  occuring 
mostly  in  the  forward  scatter  (&  ~  0*).  Therefore  it  is  not  collected  by  the  de¬ 
tector  system  at  large  angles.  Even  though  such  systems  are  not  illuminated 
by  an  infinite  plane  wave  (the  laser  beam  intensity  cross-section  is  approx¬ 
imately  Gaussian),  the  resulting  scatter  is  very  nearly  what  is  predicted  with 
plane-wave  illumination  (ref.  13). 

In  order  to  study  the  interference  of  the  scattered  radiation  with 
the  incident  plane  wave  we  investigate  a  system  with  a  simple  scattered  radia¬ 
tion  pattern:  an  a  *  0.2SX  radius,  perfectly-conducting  sphere.  Now  we  can 
distinguish  between  the  effects  due  to  the  fields  scattering  from  the  sphere 
and  the  effects  due  to  the  fields  interfering  with  the  plane  wave.  Figure  9 
shows  the  intensities  of  the  electric  fields  for  this  system  as  a  function  of 
scattering  angle  at  four  different  sphere-detector  distances.  Reducing  the  size 
of  the  sphere  relative  to  the  incident  wavelength  reduces  the  oscillatory  fre¬ 
quencies  in  the  scattering  as  a  function  of  scattering  angle  and  the  total 
amount  of  scattered  radiation.  The  interference  effects  are  best  seen  when  the 
scattered  intensities  are  comparable  to  the  incident  plane  wave  intensities 
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where  the  measurements  are  made.  We  keep  in  mind  that  the  plane  wave  in¬ 
tensity  does  not  depend  on  detector  position;  whereas,  the  total  scattered  ra¬ 
diation  decays  as  1/R*.  The  curves  of  Figure  9  are  for  a  sphere  whose  scat¬ 
tered  fields  have  only  a  small  angular  dependence. 

Figure  10  shows  the  interference  of  the  scattered  electric  fields  with 
the  incident  plane  wave  from  an  a  3  0.2SX  radius,  perfectly -conducting  sphere 
as  a  function  of  scattering  angle  at  four  different  sphere-detector  distances. 
The  intensity  on  these  graphs  is  not  normalized  for  distance.  At  a  sphere-de¬ 
tector  distance  of  infinity,  the  fields  give  the  well-known  “Rayleigh  curves" 
and  are  due  entirely  to  the  plane  wave  (It*  3  1,  Im  3  cos’d,  and  ho  3  sin’d).  As 
the  sphere-detector  distance  is  reduced,  oscillationa  about  these  values  appear 
in  the  interference  patterns, 

XL  Topographical  Plots 

It  is  apparent  that  the  near- field  interference  is  very  complex  and 
rich  with  structure.  In  order  to  view  this  complexity  from  a  different  vantage 
point,  we  generated  the  topography  of  the  electric  and  magnetic  field  intensi¬ 
ties  in  figures  11  and  12,  respectively.  These  plots  show  the  intricacies  of  the 
near-field  and  elucidate  certain  features  which  are  not  apparent  on  the  other 
figures.  Intensity  is  shown  in  the  vertical  axis.  The  sphere  lies  on  the  z  axis 
and  the  scatter  is  shown  only  to  the  left  of  the  z  axis  (in  the  x-z  or  y-z 
plane). 
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There  are  several  points  to  make  from  these  plots.  First,  the  inter¬ 
nal  fields  within  the  perfectly-conducting  sphere  are  zero.  Second,  boundary 
conditions  force  the  electric  fields  parallel  to  the  surface  (TE  and  TM  of 
Figure  tl)  and  the  magnetic  fields  perpendicular  to  the  surface  (RD  of  Figure 
12)  to  zero.  In  these  three  plots,  the  fields  outside  the  sphere  approach  zero 
as  the  sphere-detector  distance  approaches  the  sphere  radius  (R  -*■  a).  In  the 
other  three  plots  (RD  of  Figure  11  and  TE  and  TM  of  Figure  12),  there  are 
discontinuities  in  the  field  intensities  at  the  boundary  of  the  sphere.  To  satisfy 
the  boundary  condition  at  the  sphere  surface,  a  surface  current  is  induced  on 
the  sphere  equivalent  to  the  magnetic  field  at  the  surface  of  the  sphere.  Third, 
all  radial  fields  along  the  z  axis  are  zero.  Fourth,  the  TE  fields  are  identical 
to  the  TM  fields  along  the  z  axis.  Fifth,  the  TE  and  TM  field  intensities  decay 
as  1  +  l/R2.  And  sixth,  standing  waves  appear  in  the  intensity  distribution. 

The  standing  waves  that  appear  in  figures  11  and  12  are  the  result 
of  the  scattered  spherical  waves  interfering  with  the  incident  plane  waves.  To 
gain  insight  into  the  origin  and  implication  of  these  standing  waves,  we  exa¬ 
mine  the  scalar  components  along  the  z  axis.  The  amplitude  of  the  unit-nor¬ 
malized  incident  plane  wave  is  of  the  form  t"  The  scattered  spherical  waves 
are  of  the  form  Ae'**  for  negative  values  of  z  and  Be**  for  positive  values  of  z 
where  A  and  B  contain  additional  phase  and  amplitude  information  (the  spheri¬ 
cal  amplitudes  must  decay  as  1/R,  for  instance).  In  the  negative  z  direction, 
the  sum  of  these  exponentials  (e**  ♦  Ae  **)  indicates  a  decaying  standing  wave 
and  a  traveling  wave,  while  in  the  positive  z  direction  the  sum  (e**  ♦  Be**)  in¬ 
dicate  only  a  traveling  wave.  Figure  13  shows  the  electric  and  magnetic  fields 
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along  the  z  axis  for  a  perfectly-conducting,  a  =  0.25X  radius  sphere.  The 
standing  waves  appear  in  the  back-scatter  region  (negative  z).  Note  that  the 
electric  field  standing  waves,  which  have  a  node  (I  =  0)  at  the  sphere  surface, 
are  180*  out  of  phase  with  the  magnetic  field  standing  waves  which  have  an 
antinode  (I  ~  3.7S)  at  the  sphere  surface.  The  spatial  frequency  of  the  standing 
waves  is  twice  the  frequency  of  the  incident  radiation. 

We  now  consider  the  shapes  of  the  nodal  and  antinodal  regions  in 
the  scattering  plane.  These  regions  are  produced  by  the  interference  effects  of 
the  incident  plane  wave  and  the  scattered  spherical  waves.  Antinodal  regions 
occur  in  regions  where  these  two  sets  of  waves  meet  in  phase,  and  nodal  re¬ 
gions  occur  in  regions  where  they  meet  out  of  phase.  We  are  looking,  in 
effect,  at  contours  in  space  where  the  distance  from  some  reference  plane 
wave  is  equal  to  the  distance  from  the  center  of  the  spherical  plane  waves. 
This  is  precisely  the  definition  of  a  paraboloid.  Figure  14  shows  contour  maps 
of  the  TE  electric  field  intensities  and  the  TM  magnetic  field  intensities  for  a 
perfectly-conducting,  a  *  0.25X  radius  sphere  interfering  with  the  incident 
plane  wave  is  shown  in  Figure  14.  The  parabolically-shaped  standing-wave 
patterns  are  clearly  seen. 

VII. 

We  have  examined  electromagnetic  waves  scattered  from  spheres  in 
the  near-field  and  found  that  the  field  distributions  quickly  approach  those  of 
the  far-field  (within  about  10  wavelengths  for  spheres  whose  radius  is  approx- 
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imately  the  wavelength  of  the  incident  radiation).  We  examine  the  near-field 
region  where  the  scattered  fields  cannot  be  separated  from  the  incident  fields. 
The  most  striking  feature  of  the  near-field  region  is  the  parabolicaliy- shaped 
standing  waves  that  form  in  the  backscatter  region. 

Acknowledgements 

We  thank  Steve  Robinson  and  Jim  Dugan  for  many  enlightening  dis¬ 
cussions.  This  research  was  supported  in  part  by  the  Air  Force  Office  of 
Scientific  Research  and  by  an  Itek  fellowship. 


20 


Reference* 

1.  Lord  Rayleigh,  "On  the  Scattering  of  Light  by  Small  Particles,"  Phil.  Mag.  41. 

447-454  (1871). 

2.  Clebsch.  A.,  "Ueber  die  Reflexion  an  einer  Kugeiflache,”  J.  reine  angew. 
Math.  61,  19S-262  (1863). 

3.  Mie,  G..  "Beitrage  zer  Optik  tr'uber  Meiden  speziell  kolloidaler  Metallosun- 
gen,"  Ann.  Physik,  2S,  377-445  (1908). 

4.  Debye,  P.,  "Der  Lichtdruck  auf  Kugein  von  beliebigem  Material,”  Ann.  Physik, 
30,  57-136  (1909). 

5.  Lorenz,  L.,  Oeuvres  Scientiflaues  (Copenhagen,  1898). 

6.  van  de  Hulst.  H.  C.,  Light  Scattering  by  Small  Particles  (Dover,  New  York. 
1981). 

7.  Kerker,  M.,  The  Scattering  o£  Light  and  Other  Electromagnetic  Radiation. 
(Academic  Press,  New  York,  1969). 

8.  Bohren,  C.  and  Huffman,  D.,  Absorption  and  Scattering  of  Light  by  Small 
Particles  (Wiley,  New  York,  1983). 

9.  King,  R.  W.  P.,  Scattering  an<j  Diffraction  s£  Wave*  (Harvard  University 
Press,  19S9). 

10.  Boyer,  J.  M.  and  Coleman,  R.,  "Investigation*  on  the  future  development 
and  verification  of  the  E.W.F.  theory”  NASA  Technical  Report  NAS-8-11138, 
GMSFC,  Huntsville,  Alabama  (1964). 

11.  Boyer,  J.  M.,  ”A  radio  wave  mechanism  to  account  for  the  known  distribu¬ 
tion  van  Allen  belts  about  the  earth,”  Nature  190.  S97-S99  (1961). 

12.  Aden,  A.  L.  "Electromagnetic  scattering  from  sphere*  with  sizes  compar¬ 


able  to  the  wavelength.”  J.  Appl.  Phys.,  22>  601-605  (1951). 


13.  Bell,  B.  W.  and  Bickei  W.  S.,  "Single  fiber  light  scattering  matrix:  an  expe¬ 
rimental  determination,"*  Appi.  Opt.  20,  3874-3879  (1981). 

14.  Bickei,  W.  S.,  and  Bailey,  W.  M.  ,  ’"Stokes  vectors,  Mueller  matrices,  and 
polarized  light,'*  Am.  J.  Phys.  S3,  468-478  (1985). 


Figure  Captions 


Fig.  i.  The  scattering  geometry  showing  the  sphere  of  radius  a  centered  on  the 
origin. 

Fig.  2.  The  normalized  scattering  electric  field  intensities  for  an  a  3  0.525pm 
radius  quartz  sphere  illuminated  by  X  =  0.6328pm  light  for  sphere^detector 
distances  of  R  3  IX  (•—*),  R  -  3X  (O- ~ O),  R  *  9X  (X— X),  and  R  3  oo  (♦*—►). 

Fig.  3.  The  normalized  scattering  magnetic  field  intensities  for  an  a  3  0.525pm 
radius  quartz  sphere  illuminated  by  X  3  0.6328pm  light  for  sphere-detector 
distances  of  R  3  IX  (•-—*),  R  3  3X  (0—0),  R  3  9X  (X—X),  and  R  3  oo 

Fig.  4.  The  scattering  Mueller  matrix  for  an  a  3  0.52Spm  radius  quartz  sphere 
illuminated  by  X  3  0.6328pm  light  for  sphere-detector  distances  of  R  3  IX 
(*—*>,  R  3  3X  (O—O),  R  3  9X  (X—X),  and  R  3  oo  (— -). 

Fig.  5.  The  normalized  scattering  electric  field  intensities  for  an  a  3  0.52Spm 
radius,  perfectly-conducting  sphere  illuminated  by  X  3  0.6328pm  light  for 
sphere-detector  distances  of  R  3  IX  (*--■),  R  3  3X  (O—O),  R  3  9X  (X — X). 
and  R  3  ao  (*—-»). 

Fig.  6.  The  normalized  scattering  magnetic  field  intensities  for  an  a  3  0.525pm 
radius,  perfectly-conducting  sphere  illuminated  by  X  3  0.6328pm  light  for 
sphere-detector  distances  of  R  3  IX  (•—*),  R  3  3X  (O—O),  R  3  9X  (X — X). 
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Fig.  7.  The  scattering  Mueller  matrix  for  an  a  =  0.525pm  radius,  perfectly-con¬ 
ducting  sphere  illuminated  by  X  =  0.6328pm  light  for  sphere-detector  distances 
of  R  =  IX  <« — *),  R  =  3X  (O — 0),  R  =  9X  (X— X),  and  R  =  as  <« «■). 

Fig.  8.  The  scattering  Mueller  matrix  for  an  a  =  0.525pm  radius  aluminum 
sphere  illuminated  by  X  =  0.6328pm  light  for  sphere-detector  distances  of  R  = 
IX  <• — *),  R  =  3X  (O—O),  R*9X  (X— X),  and  R  *  ao  (< - ►>. 

Fig.  9.  The  normalized  scattering  electric  field  intensities  for  an  a  *  0.2SX 
radius,  perfectly-conducting  sphere  at  sphere-detector  distances  of  R  =  O.SX 
(*— •>,  R  =  0.7SX  (O—O),  R  =  1.0X  (X— X).  and  R  *  oo  («— ). 

Rg-  10.  The  normalized  scattering  electric  field  intensities  for  an  a  *  0.25X 
radius,  perfectly-conducting  sphere  interfering  with  the  incident  plane  wave  at 
sphere-detector  distances  of  R  »  O.SX  (•—*),  R  *  0.75X  (O— 1 0),  R  =  1.0X 
(X — X).  and  R  »  oo  (< - ►). 

Rf.  11.  The  normalized  topography  of  near-zone  electric  field  intensities  for  an 
a  a  0.25X  radius,  perfectly-conducting  sphere  interfering  with  the  incident 
plane  wave. 

Rf*  12.  The  normalized  topography  of  near-zone  magnetic  field  intensities  for 
an  a  «  0.25X  radius,  perfectly-conducting  sphere  interfering  with  the  incident 
plane  wave. 
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Fig.  13.  The  normalized  scattering  electric  and  magnetic  field  intensities  for  an 
a  =  0.2SX  radius,  perfectly-conducting  sphere  interfering  with  the  incident 
plane  wave  measured  on  the  z  axis. 

Fig.  14.  Contour  maps  of  the  normalized  scattering  electric  (TE)  and  magnetic 
(TM)  field  intensities  for  an  a  =  0.2SX  radius,  perfectly-conducting  sphere 
interfering  with  the  incident  plane  wave. 
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